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Abstract 


This report serves as a user's guide for the ACCESS-3 
computer program. ACCESS-3 is a research oriented program 
which combines dual methods and a collection of approximation 
concepts to achieve excellent efficiency in structural synthesis. 
The finite element method is used for structural analysis and 
dual algorithms of mathematical programming are applied in the 
design optimization procedure. 

The ACCESS-3 program retains all of the ACCESS-2 capa- 
bilities and the data preparation formats are fully compatible. 
The following new features have been added in the program: 

. four distinct optimizer options: 

. interior point penalty function method (NEWSUMT) 

. second order primal projection method (PRIMAL2) 

. second order Newton- type dual method (DUAL2 ) 

. first order gradient projection-type dual method (DUAL1) 

. pure discrete and mixed continuous-discrete design 
variable capability 

. zero order approximation of the stress constraints. 
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Approximation Concepts Code for 
Efficient Structural Synthesis 

User's Guide 


1. INTRODUCTION 

The ACCESS computer programs have been developed to 
demonstrate the effectiveness of an automated structural synthe- 
sis capability formed by combining finite element analysis 
techniques and mathematical programming algorithms using 
approximation concepts. Structures with prescribed configu- 
ration and given material properties are optimized so that 
their structural weight is minimized by modifying the sizing 
of finite elements , i . e . cross-sectional areas or thicknesses. 

The ACCESS-1 program (see Refs. 1 and 2) was designed 
to test the effectiveness of the coordinated use of approxi- 
mation concepts on problems of relatively small scale, subject 
to simple static constraints. As reported in Ref. 1, effi- 
ciency, in terms of the number of finite element structural 
analyses needed to obtain near optimal designs, was improved 
significantly over previously reported capabilities having 
comparable generality. However, many practical design pro- 
blems were beyond the capacity of ACCESS-1 and consideration 
of more complicated constraints than those treated in ACCESS-1 
was desirable. 

The ACCESS-2 computer program was developed in response 
to these needs and to build a body of experience that can be 
used to set effective guidelines for future developments of 
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large scale industrial application problems (see Refs. 3 and 
4) . Through the use of dynamic array allocation techniques 
and data transfer by sequential data files, ACCESS-2 is 
capable of treating larger problem sizes than its prede- 
cessor ACCESS-1. A thermal load analysis capability was 
added, providing experience with problems involving load 
vectors which depend on design variables. Frequency con- 
straints were also installed. In addition to the usual 
Taylor series expansion with respect to the reciprocals of 
linked design variables, additional options for represent- 
ing natural frequency constraints as first or second order 
Taylor series expansion with respect to regular linked var- 
iables were implemented. Finally the element library was 
extended. A constant strain triangular element with arbit- 
rary orthotropic material properties (CSTOR) was included 
to model laminated fiber composite material membrane struct- 
ures. A thermal shear panel element (TSP) was introduced 
to take uniform soak temperature effects into account, with 
emphasis on midplane symmetric wings. 

The ACCESS-3 program continues along these lines and 
provide a further improved structural synthesis capability 
by combining dual methods and approximation concepts (see 
Ref. 5) . The detailed analytical development of the methods 
implemented by the ACCESS-3 computer program as well as 
numerical results representing a substantial body of com- 
putational experience will be found in Ref. 6. 
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Approximation concepts are used to convert the general struct- 
ural synthesis problem into a sequence of explicit problems 
of separable algebraic form. The dual method formulation 
exploits the separable form of each approximate problem to 
construct a sequence of explicit dual functions. These dual 
functions are maximized subject to nonnegativity constraints 
on the dual variables, which are the Lagrangian multipliers 
associated with the linearized behavior constraints. 

The main advantages of the dual methods lies in their 
high level of computational efficiency, which is due to the 
fact that the dimensionality of the dual problem is rela- 
tively low for many structural optimization problems of 
practical interest. Another important advantage of the dual 
formulation is its ability to accommodate discrete design 
variables, e.g. available cross-sectional areas of bars, 
available gage sizes of sheet metal, the numbers of plies in 
a laminated composite skin, etc... The ACCESS-3 program 
provides two dual optimization algorithm options: (a) DUAL2 

(2nd order Newton type of algorithm) , which is restricted 
to pure continuous design variable problems; and (b) DUAL1 
(1st order gradient projection type of algorithm) , which can 
solve pure discrete and mixed continuous-discrete design 
variable problems. 

In addition, a 2nd order primal projection algorithm 
called PRIMAL2 has been introduced, which operates, like 
NEWSUMT, on each explicit approximate problem expressed in 
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terms of the reciprocal design variables. Hence a collec- 
tion of four distinct optimizer options are available in 
the ACCESS-3 computer program: NEWSUMT, PRIMAL2 , DUAL2 

and DUAL1. 

Another new feature of ACCESS-3 is that the stress 
constraints can be replaced with zero order explicit ap- 
proximations instead of first order ones. The zero order 
approximations are obtained using classical stress ratio 
formulas. They can be expressed as simple side constraints, 
which is especially beneficial when dual methods are em- 
ployed. A selection criterion permits automatic subdivi- 
sion of the stress constraints in two categories: those 

requiring first order approximation (full linear Taylor 
series) and those for which zero order approximation (side 
constraint) is accurate enough (see Ref. 7) . 

In summary, the main feature of ACCESS-3 lies in the 
joining together of approximation concepts and dual methods. 
This solution scheme can be interpreted as a generalized 
optimality criteria method. Another new capability is the 
zero order approximation of the stress constraints based 
on the conventional "Fully Stressed Design" optimality 
criterion. Therefore the ACCESS-3 program can be regarded 
as an advanced research tool where mathematical programming 
and optimality criteria approaches coalesce to provide an 
efficient and reliable structural weight minimization method. 
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2. THE ACCESS-3 COMPUTER PROGRAM 
2 . 1 Program organization 

The fundamental structure of the ACCESS-3 program is 
outlined in Fig. 1. Upon activation the preprocessor reads 
and prints out the input data in a readable format. The pre- 
processor then computes all the ancillary data that are inde- 
pendent of changes in the design variables and it stores the 
results in appropriate arrays as well as in temporary external 
files (see Table 1) . When preprocessing is completed success- 
fully, the design process control (DPC) block is activated 
and it initializes the design iteration process. At the out- 
set the design given in the input data is transferred to the 
approximate problem generator (APG) , and this design is ana- 
lyzed by the finite element method. Constraint functions are 
evaluated using the response quantities obtained from the 
finite element analysis and then the initial set of critical 
and potentially critical constraints is identified and tagged. 
Explicit approximate expressions for these tagged constraints 
are computed using the Taylor series expansion with respect 
to appropriate intermediate design variables. Reciprocals 
of independent design variables are used as intermediate 
variables throughout the program, except for an optional use 
of the independent design variables themselves when expanding 
frequency constraints. In ACCESS-3, the objective function is 
the structural weight and it may be expressed exactly and expli- 
citly in terms of the independent design variables or their 
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reciprocals. Thus, the APG block can generate an approximate 
problem statement of the form: 

Minimize W(X) 

X = (x^,X2/ • • 

Subject to 

H (X) ^ 0 q = 1,2, . . .Q 

where W(X) and all H (X) are explicit analytic functions of X. 
Note that the number of constraints Q for this approximate 
problem is much smaller than that of the original structural 
design problem, because only the tagged constraints are included 
and all other constraints are temporarily ignored during a 
particular design stage. 

The data which define the approximate problem are sent 
back to DPC and subsequently given to the optimization algorithm 
block (OA) . The primary function of OA is to carry out a 
numerical search process which will improve the design by 
operation on the current approximate problem statement. Since 
OA deals with problems that are stated in algebraically expli- 
cit form, it is not even aware that these problems are related 
to structural design. Therefore, any established algorithm 
for inequality constrained minimization of a function of many 
variables may be used. However the main feature of ACCESS-3 
is that special purpose, highly efficient OA have been selec- 
ted, which capitalize upon the special mathematical structure 
of the explicit problems generated by APG. Unlike its prede- 
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cessor ACCESS-2, in which only one general purpose OA was 
implemented (NEWSUMT) , the new ACCESS-3 program offers a collec- 
tion of four available options (NEWSUMT, PRIMAL2 , DUAL2 and 
DUALl) . The user can select an OA from this collection, tak- 
ing into account the size of the problem, the nature of the 
constraints and the definition of the design variables (con- 
tinuous or discrete). Section 3.8 gives a description of the 
available OA. 

After carrying out a numerical search with the approxi- 
mate problem, the optimization algorithm (OA) block proposes 
an improved design x' to DPC. This step completes one stage 
of the design iteration procedure. 

In summary, one stage of iteration includes one finite 
element structural analysis, one constraint deletion process, 
one sensitivity evaluation for retained constraints, and one 
optimization of an approximate problem. Since the final 
design is subjected to a detailed finite element analysis, 
the total number of finite element analyses equal the number 
of iteration stages plus one, which will be typically around 
10. The iterative design process is terminated when one of 
the specified convergence criteria is satisfied. 

2 . 2 Program Implementation 

The ACCESS-3 computer program may be executed as a stand 
alone program. It consists of approximately 15000 FORTRAN 
statements. Two operational versions of the program exist, 
one for IBM 360/370 systems and the other for CDC 6600/7600 
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systems. Since it contains no machine dependent statements, 
it can be made operational on various computers, provided 
enough main memory capacity and auxiliary data storage support 
are available. 

Auxiliary storage files are required as shown in Table 1. 
Files 10, 11, 12, 13, 14 and 15 are required for all problems. 
File 16 is required only when type 4 elements (TSP) are used 
in the structural model. Files 18, 19, 20, 21 and 2.2 are 
required only when second order expansions of frequency con- 
straints are specified. File 17 is required only when zero 
order approximation of stress constraints is specified. 

The required size of blank common is very problem depen- 
dent: i.e. it depends on the structural analysis model (number 
of nodes, elements and load conditions) , the number of inde- 
pendent design variables, and the constraint types included. 

For certain problems, it also depends on the initial design. 
Hence, it is rather difficult to give explicit formulas which 
estimate the size of blank common requirements. Table 2 gives 
actual blank common array size requirements for several example 
problems . 

Overlay or segmentation of the program can be designed 
easily by referring to Fig. 2. The simple 3 level overlay is 
adequate to solve most of the meaningful problems. If an 
operating system allows more flexible overlay structure, it is 
possible to decrease the core requirement further. 
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However, the net gain acquired by the elaborate overlay may 
not be significant, since most of the core is used for data 
and not for instructions. 

All routines are written in standard FORTRAN IV language 
and they have been tested on: 

(a) the IBM 360/91 using the FORTRAN-H compiler at UCLA 

(b) the CDC 6600 at the NASA Langley Research Center. 
Implementation on other computers will be straightforward 
provided those computers have the required main memory capacity 
Except for the blank common arrays, 330^ Q k and 95-^k bytes are 
required on IBM 360/91 without and with program overlay, res- 
pectively. On a CDC 6600, the corresponding basic memory 
requirement is 100gk words with overlay. 

2 . 3 Restrictions and Limitations 

The amount of main memory storage required for solution 
of a particular problem depends upon many factors, including 
the number of nodes, the number of elements, the number of 
design variables, the element types used, the kinds of con- 
straints imposed, and even the initial design employed, etc. 

For static problems, it is necessary to retain two system 
stiffness matrices and the load vectors in core. For dynamic 
problems, three system matrices must be retained in core. If 
a problem involves dynamic constraints and thermal shear panel 
elements, four system matrices must be in core simultaneously. 
Also a complete approximate problem statement (all retained 
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constraint values and all the corresponding derivative com- 
ponents) must be in core for the OA block. It is difficult to 
estimate the array size required for a system stiffness matrix 
in advance. Only the nonzero skyline of an upper half matrix 
is stored, hence the memory requirement depends on the node 
numbering scheme. For medium size problems (300-600 DOF) , 
the density of nonzero elements in the matrix is usually 20-50% 
and a first approximation can be made by estimating the density 
based on observation of the finite element model. The main 
memory storage required for the integer portion of the blank 
COMMON is usually less than 10,000 words, but the real variable 
portion is very dependent on the nature of the problem. For 
problems in which the number of constraints retained tends to 
be larger and in which there are many independent design varia- 
bles (e.g. structures involving laminated fiber composite skins) 
the constraint derivative array size [i.e. (Number of Design 
Variables) x (Number of Constraints Retained) ] may limit the 
problem size, since this large array must be in core in addition 
to the instructions and local variables. Furthermore, a large 
number of available discrete values may also limit the problem 
size, since a separate set of discrete values must be associ- 
ated to each independent linked design variable, in view of 
the design variable normalization process used in ACCESS-3. 

When first order approximations are used, frequency con- 
straints can be imposed on any subset of frequencies within 
the lowest NFREQ frequencies. If second order approximations 
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are employed, all frequencies in the lowest NFREQ frequencies 
must be bounded. 

Capabilities for aeroelastic constraints are not available 
in this version, therefore NMODE must be zero and the flight 
condition specification flags must all be zero. 

There are 4 available optimization algorithms in ACCESS-3. 
However only NEWSUMT is a general purpose optimizer, which can 
be used even when the explicit approximate constraints are not 
linear in the reciprocal design variables. As a consequence, 
the possible combinations of frequency constraint approxima- 
tions and optimization algorithms are not all acceptable; 

Table 3 shows the available combinations that may be used. 

It is also important to recognize that only the DUAL1 
optimization algorithm is applicable to problems involving dis- 
crete design variables. For the case of pure continuous design 
variable problems, all four optimizers are applicable. The 
algorithm options available for various kinds of problems are 
summarized in Table 4. 
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3. INPUT DATA PREPARATION 


It is assumed that the reader is familiar with elastic 
structural analysis via the finite element displacement method, 
as well as with associated structural modelling techniques 
and typical data preparation procedures. Sufficient informa- 
tion for preparing the input data card images is given in 
Section 4, therefore, the explanations given in this section 
are limited to topics which require somewhat detailed technical 
discussion in order to avoid possible misunderstandings. 

3 . 1 System of Units 

Input data of the ACCESS-3 computer program may be pre- 
pared in any system of units as long as they are consistent. 

For example, if it is decided that the units for length and 
force are to be centimeters and Newtons, respectively, then 

the corresponding units for pressure load or allowable stress 
o 

must be N/cm . Note that the material constant specification 
calls for the specific weight of the material, not its mass 
density. To be consistent lumped nodal mass should be given 
using weight rather than mass units. Example problems given 
in Appendix B, as well as the corresponding computer input data, 
are presented using numerical values associated with the US 
units, simply because all the examples were originally presented 
in the literature using US units. 
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3 . 2 Node Numbering Scheme 


The system stiffness and mass matrices are stored in 
a vector form within the skyline of the nonzero elements, i.e. 
there are no operations or no storage allocations with elements 
that remain zero during the solution. This scheme allows 
somewhat more flexible node numbering arrangement than the 
ordinary band equation solver. It is better, however, to follow 
the same guidelines in preparing data as for a banded matrix 
solution scheme; i.e., differences among node numbers associ- 
ated with an element must be kept as small as possible for all 
elements . 

3 . 3 Symmetric Wing Model 

If the webs of a midplane symmetric wing are modelled 
with SSP elements, only the upper (or lower) half of the wing 
is modelled. Assuming that the X-Y reference plane is the 
plane of symmetry, the X and Y displacement components and 
loading components are then anti-symmetric. Displacements and 
loadings in the Z direction are identical for both sides of 
the X-Y plane. For example, if a cantilever beam such as that 
shown in Fig. 3(a) is to be modelled using two SSP elements, 
then the simplified model should be that shown in Fig. 3(b). 

Note that only half of the load P needs to be applied to the 
node 3, since the other half is implicitly applied to the con- 
jugate node 3' (which does not exist explicitly in the model). 
The SSP elements are always perpendicular to the X-Y plane of 
symmetry. 
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The assumed displacement function for SSP elements cannot 
accommodate uniform thermal expansion of each SSP element. If 
specified midplane symmetric temperature changes are specified 
for a midplane symmetric structural model, in which the verti- 
cal webs are represented by SSP elements, ACCESS 3 branches 
and makes a separate calculation which adds in the midplane 
symmetric temperature change effects. This is accomplished by 
assembling equilibrium equations for the midplane symmetric 
structure with all of the SSP elements replaced by TSP elements 
while only considering midplane symmetric temperature change 
loading. These equilibrium equations are solved for displace- 
ments S due only to midplane symmetric temperature changes . 
These thermally induced midplane symmetric displacements are 
superimposed on the previously computed midplane antisymmetric 
displacement state due to mechanical loads only. Treating the 
symmetric and antisymmetric contributions separately reduces 
the number of displacement degrees of freedom that need to be 
considered in each of the two analyses and for thin wings it 
also tends to improve the accuracy of the analysis by avoiding 
the poor conditioning often associated with simultaneous treat- 
ment of bending and membrane response. The strain state is 
computed based on the total displacement, and the stress state 
is computed by transforming the strain state using the stress- 
strain relationships. 
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3 . 4 Design Variable Linking and Stress Constraint Regionali - 
zation 

The general concept of design variable linking is dis- 
cussed in Sec. 2.3.1 of Ref. 1. In the ACCESS-3 computer pro- 
gram, if the sizes of some group of finite elements of the same 
type are controlled by a single design variable, these elements 
are said to belong to the same design variable linking group. 

The sizes of elements in a design variable linking group are 
modified in proportion to the initial sizes given in the input 
data. 

Design variable linking groups are also used to define 
"regions" for the regionalization of stress constraints. The 
general idea of regionalization is described in Sec. 2.4.1 of 
Ref. 1. Elements which belong to the same design variable 
linking group form a region and only one stress constraint per 
load condition (the most critical) is considered for each group 
in any stage of the iterative design procedure. Selection of 
the critical stress constraints within a region is not rigidly 
fixed, but dynamically updated at the beginning of each stage. 

If the location of the critical stress constraints shifts fre- 
quently within a region from stage to stage the iteration pro- 
cess may be unstable, although this type of instability was 
not observed in solving any of the problems given in Ref. 1. 
However, if the user desires to remove the regionalization of 
stress constraints, it is only necessary to specify IGLINK = -200. 
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3.5 Failure Criteria 


The CSTOR element is implemented to model structures made 
with orthotropic materials including multi-layered fiber com- 
posite laminates. While strength failure criteria for iso- 
tropic metal alloy materials are imposed using the von_Mises 
combined effective stress, strength failure criteria for CSTOR 
elements are selected from 3 available options. 

They are: 

A. Maximum strain criteria 

S L * e L - a L AT S 5 L 
^ e T - a T AT £ e£ 

Y LT^ “ Y LT 


B. Stress interaction formulas 
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C. Tsai-Azzi Criterion 
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: allowable shear stress 


longitudinal tensile stress 
longitudinal compressive 

transverse tensile stress 
transverse compressive 


: Poisson's ratio relating to contraction in 

the longitudinal direction due to extension 
in the in-plane transverse direction 

: Poisson's ratio relating to contraction in 

the in-plane transverse direction due to 
extension in the longitudinal direction 


Among the three alternative strength criteria, the maximum 
strain criterion is the most conservative while the stress inter 
action formulas are usually the least conservative. 
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3 . 6 Computation of Constraints 


All constraints, except the side constraints, are normal- 
ized so that potentially critical constraint functions in the 
feasible region assume values between 0.0 and 1.0. Constraint 
functions are defined as follows: 

5i2Ei52§IB22£_Q2I}2l:£2i2£§ 
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where d^ is the projection of the distance between the two 
points on a plane normal to the displacement components 6^ 
and 6 2 ’ 

Stress (Strain) Constraints 


a< U > - a 
„(U) 


a - a 


(L) 


(L) 


^ 0 


k 0 


18 




For interaction formulas and Tsai-Azzi failure 


criteria, see Section 3.5, B and C, 
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3 . 7 Zero Order Approximation of the Stress Constraints 

It is well known that in a structural synthesis problem, 
the stress constraints can often be efficiently treated using 
the classical "Fully Stressed Design" (FSD) concept. In this 
approach a stress ratio formula is employed to transform the 
stress constraints into simple side constraints, which can be 
interpreted as zero order explicit approximations (see Ref. 6, 
page 39) . 

The zero order approximation of stress constraints leads 
to a significant reduction in the number of behavior con- 
straints retained in each explicit approximate problem. This 
feature is particularly beneficial when dual methods are 
employed, because the dimensionality of the dual problem 
corresponds to the number of constraints represented by first 
order approximations. On the other hand, the FSD procedure 
does not always converge to the true optimum and is sometimes 
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the source of instability or divergence of the optimization 
process. However, in practical structures, it is observed 
that many of the stress constraints can be approximated with 
sufficient accuracy by zero order explicit approximations 
using FSD, while others require a more accurate approximation, 
using first order Taylor series expansion with respect to the 
reciprocal design variables. 

The ACCESS-3 program provides the capability of selecting 
automatically the stress constraints for which a zero order 
approximation by FSD is sufficiently accurate. For each 
retained potentially critical stress constraint, the following 
test is accomplished (see Ref. 6) ; 


STR - LRDV X GRD 
STR 


<1 EPS 


STR denotes the reference value describing the stress state, 
i.e., the tensile or compressive stress in a TRUSS element, 
the Von Mises combined effective stress in a CSTIS element, 
the longitudinal, transverse or shear strain in a CSTOR 
element, etc.... (see Section 3.5). LRDV represents the linked 
reciprocal design variable describing the element in which 
the current STR is evaluated 7 *, GRD stands for the first par- 
tial derivative of STR with respect to LRDV. 

All the stress constraints for which the test above is 
satisfied will be replaced with zero order explicit approxima- 

^Note that LRDV = 1 in view of the normalization process 
used in ACCESS-3. 
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tions using stress ratio formulas (side constraints) while the 
others continue to be transformed into first order explicit 
approximations using Taylor series expansion (full linear 
constraints) • Of course this selection of zero/first order 
approximations for the stress constraints must be repeated at 
each design stage, exactly like the well known truncation 
procedure for selecting the potentially critical constraints. 

The severity of the test depends on the value adopted 
for the tolerance EPS, which is provided by the user. If 
EPS is taken close to 1 , a small number of stress constraints 
will be first order approximated. The smaller is the value of 
EPS, the larger will be the number of first order approximated 
stress constraints. 

3 . 8 Optimization Algorithms 

The ACCESS-3 program includes four distinct optimization 
algorithms (OA) that the user can select depending upon the 
nature of the constraints , the expected number of strictly 
critical first order approximated constraints, the number of 
design variables, and their continuous or discrete character. 

NEWSUMT_Optimizer 

The NEWSUMT optimization algorithm is the same as in the 
ACCESS-2 program, where it was the only available option. 

This optimizer implements a sequence of unconstrained minimi- 
zations techniques using a modified Newton's method and a 
quadratic extended penalty function feature to facilitate the 
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unconstrained minimizations. One virtue of this interior 
penalty function type of OA is that it can usually be controlled 
so as to provide an improved design that is also feasible with 
respect to all of the constraints at each stage in the design 
process. NEWSUMT is thus particularly interesting when the 
constraints of the primary problem are highly nonlinear in the 
reciprocal variables, in which case each approximate problem 
must be solved only partially to preserve the design feasi- 
bility. 

Another advantage of the NEWSUMT optimizer lies in its 
generality. Unlike the other optimization algorithms avail- 
able in ACCESS 3, NEWSUMT can indeed accommodate explicit con- 
straints which are not linear in the reciprocal design varia- 
bles. As a result, NEWSUMT must be employed when second order 
Taylor series expansions are generated to represent the fre- 
quency constraints . 

From the point of view of the computational cost, how- 
ever, the NEWSUMT optimizer is far less efficient than the 
other available options. It should be selected only in the 
special cases previously indicated. 

PRIMAL2_0ptimizer 

PRIMAL2 is a second order projection algorithm for pro- 
blems with separable objective function and linear constraints. 

It uses a weighted projection operator to generate a sequence 
of Newton's search directions that are constrained to reside 
in the subspace defined by the set of active constraint hyper- 
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planes (see Ref. 6, page 80 ) . Like NEWSUMT, the PRIMAL2 opti- 
mizer generates a sequence of improved feasible designs with 
respect to the explicit approximate problem. Hence PRIMAL2 
can be adequately used for seeking a partial solution to each 
approximate problem, in such a way that the constraints of 
the primary problem remain almost satisfied. This is accom- 
plished by prescribing an upper limit on the number of one- 
dimensional minimizations performed before updating the approx 
imate problem statement. Of course, PRIMAL2 is also well 
suited to solve exactly each approximate problem, in which 
case it produces the same iteration history as the dual met- 
hods, with comparable efficiency. 

PRIMAL2 is computationally more economical than NEWSUMT. 
However, in the current version of the ACCESS-3 program, the 
PRIMAL2 algorithm has not been tested extensively enough and 
it should be used with circumspection. 

DUAL2 Optimizer 

The dual method formulation, which exploits the separable 
form of the approximate problem, consists in maximizing the 
explicit dual function subject to nonnegativity constraints on 
the dual variables. The efficiency of this approach is due 
to the fact that the dimensionality of the dual space, which 
is primarily dependent on the number of critical behavior 
constraints, is relatively low for many structural optimiza- 
tion problems of practical interest. In contrast to the 
primal optimizers NEWSUMT and PRIMAL2 , which usually seek a 
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partial solution to each approximate problem, reducing the 
weight while remaining feasible, the dual methods efficiently 
find the "exact" minimum weight solution of each separable 
approximate problem (see H2.4.1, Ref. 6). Therefore, at the 
end of any stage, the design may not be strictly feasible, 
in which case scale up is needed to obtain a feasible design. 
It should be noted, however, that the design infeasibility, 
if any, is usually small and decreases stage by stage. 

DUAL2 is a dual method which employs a second order 
Newton type of algorithm to find the maximum of the dual 
function when all the design variables are continuous . Since 
the DUAL2 optimizer has been found to be very efficient in 
practice, it is the recommended option for pure continuous 
varialbe problems. 

DUALl_Optimizer 

DUAL1 is a dual method which employs a specially de- 
vised first order gradient projection type of algorithm to 
find the maximum of the dual function when the design vari- 
ables are all discrete or mixed continuous^discrete . The 
DUAL algorithm incorporates special features for handling 
the dual function gradient discontinuities that arise from 
the primal discrete variables. These discontinuities occur 
on specific hyperplanes in the dual space. The DUAL algo- 
rithm determines usable search directions by projecting the 
dual function gradient on the intersection of the successively 
encountered first order discontinuity planes. When a maximum 


24 



has been obtained, the algorithm is restarted releasing all 
of the previously accumulated discontinuity planes. The 
whole maximization process is terminated when two successive 
restarts yield the same dual point. When all the design 
variables are continuous, the DUAL1 algorithm reduces to a 
special form of the conjugate gradient method; however it 
is generally less efficient than the DUAL2 optimizer for 
pure continuous variable problems. 

Con trol_over_ Convergence 

Strictly speaking each new design generated by the 
optimization algorithm is an improved design only with re- 
spect to the approximate problem statement. When this new 
design is analysed by means of the finite element method, 
it may turn out that some behavior constraints are violated. 
This situation may occur when the design changes in one stage 
exceed the applicable range of the approximate problem state- 
ment. 

The NEWSUMT optimizer is capable of locating feasible 
designs starting from an infeasible design; however vio- 
lation of constraints usually has a deleterious effect on the 
convergence characteristics. The PRIMAL 2 optimizer contains 
a built-in scaling procedure which readily finds a feasible 
critical design starting from any design (n.b. this scaling 
process is only applicable to static constraint violations) . 
It is worthwhile noticing that constraint violation can 
usually be controlled or eliminated by appropriate use of 
the maximum step size parameter STEPMX (move limit) and its 
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dynamic modification feature via the parameters STEPMX- 
multiplier and STEPMX-lower limit (see Section 4; block of 
data XXIV) . Furthermore, constraint violation is less 
likely to occur if the optimization algorithms are used to 
seek a design improvement that corresponds to only a partial 
solution of each approximate problem. When the NEWSUMT 
option is selected, this can be done by reducing the number 
of response surfaces and/or increasing the response factor 
cut ratio (see Ref. 8) . When the PRIMAL2 optimizer is used, 
the same effect is obtained by reducing the maximum allow- 
able number of one-dimensional minimizations. 

On the other hand, the dual optimizers DUALl and DUAL2 
generate a sequence of not necessarily feasible designs and 
their performance is not adversely affected when design in- 
feasibility is encountered. However, if for some reason the 
user wants to control constraint violation, this can still 
be done through the STEPMX parameters. 

3. 9 Printout Control Parameters 

There are two parameters used to control the line printer 
output quantity, namely IPRINT and JPRINT. The greater the 
integer numbers assigned to these parameters, the more de- 
tailed output will be printed. IPRINT controls printout 
from all programs except the optimizers. Brief summaries 
of the output items are given in Table 5. Standard output 
will be obtained from the optimizers (see Tables 6,7,8 and 
9) . The standard value is JPRINT = 0 for all optimizers. 
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4. INPUT DATA DESCRIPTION 


All input data are read in with fixed format, hence 
column positions of the punched data are of critical import- 
ance. Especially note that all blank columns are regarded 
as zeroes for numerical inputs. 

I. Job description and heading cards (II, 79A1) 

The first column is used as follows 

0 or blank: ordinary heading cards, whose contents 

in columns 2-80 will be printed on the 
page of the output. 

1 : indicates tl^at this is the last heading 

card and input data cards follow. 

2 : request for immediate normal termination 

of this job. 

Any number of cards may be used to describe or to comment 
the job. Note that the last heading card must have "1" 
punched in the first column. Without this, all of the 
data may be regarded as heading cards. 

II. Primary control cards 

Card 1 (715) 

IOPT : 1 = Input data check only 

2 = Structural analysis only 

3 = Structural analysis and constraint 

function evaluation 
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see section 3.8 


IPRINT 


IGLINK 


IANALY(l) : 


IANALY ( 2 ) : 


IANALY ( 3 ) : 


IANALY ( 4 ) : 


! 4 = optimization by the NEWSUMT optimizer 

5 = optimization by the PRIMAL2 optimizer 

6 = optimization by the DUAL1 optimizer 

7 = optimization by the DUAL2 optimizer 

Printout control parameter except for 
output from each optimizer 
Standard output: IPRINT=2 (see Table 5) . 

0 = standard execution 

-200 = removal of stress constraint region- 
alization 

-300 = removal of stress constraint region- 
alization for fixed size elements 
only 

1 = Compute displacement 

0 = Skip displacement calculation 

1 = Compute stress/strain for all elements 

0 = Skip stress/strain calculation 

1 = Compute eigenanalysis 

0 = Skip eigenanalysis 

0 always 
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Card 2 (1215) 


IN 

IBN 

INL 

IMATIS 

IMATOR 

INITVG 

ILOWBG 

IUPPBG 

ITHLDG 

IPRLDG 

IDISVG 

NVAL 


Total number of nodes 

Number of boundary nodes 

Number of load conditions 

Number of isotropic materials 

Number of orthotropic materials 

Number of initial value groups for design 

variables 

Number of minimum size groups for design 
variables 

Number of maximum size groups for design 
variables 

Number of thermal load groups 

Number of pressure load groups 

Number of discrete value groups for design 

variables 

Maximum number of available discrete values 
in each IDISVG discrete value group 
Default_option ; if NVAL=0, the program 
adopts NVAL=20 


Card 3 (1015) 

IETP(i) , i = 1,2, . . .10 


Number of elements in the i-th element type 
i = 1 TRUSS 

2 CSTIS (CST isotropic) 

3 CSTOR (CST orthotropic) 
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4 SSP ( symmetric shear panel) 

5 PSP (pure shear panel) 

6 TSP (thermal shear panel) 


III. Node Coordinates 


(15 , 5X, 3E10.4) 


IN cards are required to specify the node coordinates of 
node numbers 1 through IN. The order of the cards may 
be random. 

n. : Number of the i-th node 


IV. 


l 


X . 
ni 

: X 

coordinate 

of 

the 

node 

n . 
l 

Y , 
ni 

: Y 

coordinate 

of 

the 

node 

n . 

l 

Z 

ni 

: Z 

coordinate 

of 

the 

node 

n . 

l 


Boundary Conditions (415) 

If all 3 degrees of freedom associated with a node are 
free, the node is not a boundary node. Otherwise it is 
a boundary node and for each boundary node, a card is 
required. 

bn^ : i-th boundary node number 


V. 


IBX 


IBY. 


bn . 

l 

bn . 


IBZ 


bn . 

l 


Element Data 


\ constraint code: 0 

) 1 


free 

fixed 


If IETP(i) = 0 for the i-th element type, no data is 
required. For dach element type with IEPT(i)^ 0, 
IETP(i)+l cards are required. 

Card 1 } element type number (15) 
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Card 2-IETP ( i) +1 


: element information (1215) 


NP 

NQ 

NR 

NS 

LGN 

IGN 

LBGN 

UBGN 

MTLGN 


SCC 


DVGN 


element number 

node number corresponding to the internal 

node number P 

node number Q 

node number R 

node number S 

linking group number, = 0 for the fixed size 
elements 

initial value group number 
lower bound group number 
upper bound group number 
material group number 

> 0 for isotropic materials: 1,2,... 

< 0 for orthotropic materials: -1,-2,... 

side constraint code 

-1 : element size restricted by the lower 

bound only 

0 : non-negativity constraint only 

1 : element size restricted by the upper 

bound only 

2 : element size restricted by both lower 

and upper bounds 

discrete value group number, = 0 for the 
continuous design variables. DVGN^O is 
valid only with DUALl optimizer (I0PT=6) . 
§E22i2i_2Eti22 : if DVGN = 999 , the program 
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adopts the following set of available dis- 
crete values for element M^:" 

DMIN, 2xDMIN , 3 X DMIN, . . . ,NVAL xDMIN 
where DMIN denotes the minimum gauge for 
element (depending upon the LBGN value 
defined on the same data card) . The special 
option DVGN=999 is useful for fiber compo- 
site material problems, where DMIN repre- 
sents the smallest change in lamina thick- 
ness. For example, if the laminates are 
required to be symmetric, DMIN will be equal 
to the thickness of two plies. 

Comments 

1. Elements must be numbered starting from 1 through 
IETP(i) for each element type. For example, if a 
structure is modeled using 100 TRUSS elements and 
300 CST elements, TRUSS element numbers are 

1,2, 3... 100 and CST element numbers are 1 , 2 , 3 , . . . 300 . 
Within an element type, order of element data cards 
may be random, 

2. NR and/or NS are not required for element types with 
only 2 or 3 nodes per element. 

3. LGN, linked group number , starts from 1 for each 
element type. For example, if a structure is 
modeled with 100 TRUSS and 300 CST elements, with 
10 and 30 design variables allocated to TRUSS and 
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CST, respectively, then the linked group number 
for TRUSS runs from 1 through 10 and that for CST 
ranges from 1 through 30. 


VI. Initial Values (7E10.4) 

INITVG real numbers must be given. If INITVG > 7, two 
or more cards are required. The first value of the 
first card indicates the initial value for the group 
number 1, and so on. 

VII. Lower Bound Values (7E10.4) 

Minimum gauge values. ILOWBG real numbers must be given. 
If ILOWBG > 7, two or more cards are required. If 
ILOWBG = 0, no card is required. 

VIII. Upper Bound Values (7E10.4) 

Maximum gauge values. IUPPBG real numbers must be given. 
If IUPPBG > 7, two or more cards are required. If 
IUPPBG = 0, no card is required. 

IX. Available Discrete Values (7E10 . 4/7E10 . 4 , f . ) 

IDISVGxNVAL real numbers must be given, i.e.: for each 

of the IDISVG discrete value groups, NVAL discrete values 
must be specified in ascending order. If, for a given 
group, the number of available discrete values, say 
NVALG , is less than NVAL, the remaining ( NVAL-NVALG ) 
positions must be left as blank. 

If IDISVGxNVAL > 7, two or more cards are required. If 
IDISVG = 0, no card is required. 
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X. Isotropic Material Data (6E10.4) 

IMATIS cards are required and on each card the following 
6 real numbers must be given: 

Elastic modulus 
Poisson's ratio 
Specific weight 
Thermal expansion coefficient 
Allowable compressive stress 
Allowable tensile stress 

XI. Orthotropic Material Data ( 7E10 . 4/7E10 . 4/6E10 . 4 ) 

IMATOR x 3 cards are required, i.e. for each material 
group, 3 cards are required, containing the following 
data. 

Card 1 


E 

V 

Y 
a 



LT 


v 

Y 

a l 

a„ 


LT 


Card 2 


m. 


n T 


-t 

e L 


: Longitudinal elastic modulus 

: Transverse elastic modulus 

: Shear modulus 

: Longitudinal Poisson's ratio 

: Specific weight 

: Longitudinal thermal expansion coefficient 

: Transverse thermal expansion coefficient 

: \ Direction cosines of the longitudinal 
: > axis with respect to system reference 
: / coordinates 

: Tensile allowable longitudinal strain 
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Compressive allowable longitudinal strain 
Tensile allowable transverse strain 
Compressive allowable transverse strain 


Allowable shear strain 
Tensile allowable longitudinal stress 
Compressive allowable longitudinal stress 
Tensile allowable transverse stress 
Compressive allowable transverse stress 
Shear allowable stress 


1. The transverse Poisson's ratio v TL is internally 
computed using the relation v mT E T = v Tm E_ 

1 JLj Li Li 1 

2. Depending upon the failure criteria applied to the 
specific material, either strain allowables or 
stress allowables are left unspecified. Failure 
criteria options will be specified later in the 
category XXI. 

XII. Lumped Nodal Loads 

Two card groups are required to specify lumped nodal loads 
applied to the structure. 

Card Group 1 (1415) 

Number of nodes subject to lumped nodal loads for each 
load conditions. INL integer numbers must be given. 
Card Group 2 (15, 5X, 3E10.4) 

For each load condition, the specified number (by the 



Card 3 

^LT 



Comments 
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group 1 cards) of cards must be given to identify 
the node numbers and associated load components in 
the reference coordinate system. 

Supply one blank card if there is no lumped nodal load. 
XIII. Pressure Load Data 

No card is required if IPRLDG =0. If IPRLDG > 0, the 
following 5 groups of cards must be given. 

Card Group 1 (1015) 

Number of elements subject to pressure load for each 
element type. (Presently, only CSTIS and CSTOR 
elements can be subject to pressure loads) . 

Card Group 2 (1415) 

Pressure load ON-OFF for each load condition, 
k 

ONOFF =0 No pressure load for load condition k 
= 1 Pressure load should be considered 
for the k-th load condition. 

Card Group 3 (1415) 

Element numbers subject to pressure loads for all 
member types mentioned in card group 1. For each 
element type, the first element number subject to 
pressure load must be punched in columns 1-5; 
namely the group 3 cards should be subgrouped for 
different element types. 

Card Group 4 (1415) 

For each load condition corresponding to a load con- 
• • • k 

dition with ONOFF =1, an identical amount of data 
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similar to that specified in the card group 3 must be 
given. Those numbers designate the pressure magni- 
tude group numbers, which are the pointers to the 
pressure magnitude applied to the corresponding ele- 
ment type and element number. This set of cards 
should be given for all load conditions with 
ONOFF k = 1. 

Card Group 5 (7E10.4) 

Pressure load magnitude for each pressure load group 
must be given. IPRLDG real numbers are required. 

Comments ; 

1. The direction of the pressure force is determined 
by the node numbering scheme of the triangular 
element and also by the sign of the pressure load 
magnitude specified in the card group 5. When 
the P, Q and R nodes of the triangle are in counter 
clockwise order and the corresponding pressure 
magnitude has a positive sign, positive pressure 
is applied to the surface of the triangular region. 

P P 




Pressure Load Sign Convention 
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2. Pressure applied on a single triangular surface 

must be uniform; no variation of pressure over an 
element surface can be represented. 

XIV. Inertia Load Data 

Self-weight in a gravitational field or uniform trans- 
lational acceleration will be accounted for by specifying 
this set of data. Note that rotational inertia loads 
cannot be considered. Two groups of cards are required. 
Card Group 1 (1415) 

INERTL k : k = 1 , 2 , . . . INL 

Inertia load ON-OFF for each load condi- 
tion. 

1 ; Inertia load exists 

2 : No inertia load for the load condi- 

tion 

Card Group 2 (4E10.4) 

For each load condition with INERTL 0, one card 
will be required. 

ACC : Magnitude of acceleration in units of the 
standard earth gravitational field 
(e.g. 4g) 

X : 

Y : 

Z : 

Supply one blank card if there is no inertia load. 


Direction cosine components of the 
acceleration vector in the reference 
coordinate system. 
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XV. 


Thermal Load Data 


No card is required if ITHLDG =0. If ITHLDG > 0, the 
following 5 groups of cards must be given. 

Card Group 1 (1015) 

Number of elements subject to thermal load for each 
element type. 

Card Group 2 (1415) 

Thermal load ON-OFF for each load condition. 

ON-OFF =0 No thermal load for load condition k 
1 Thermal load should be considered for 
the k-th load condition. 

Card Group 3 (1415) 

Element numbers subject to thermal loads for all 
member types mentioned in card group 1. For each 
element type, the first element number subject to 
thermal load must be punched in columns 1-5; namely 
the group 3 cards should be subgrouped for different 
element types. 

Card Group 4 (1415) 

For each load condition corresponding to the load 

k 

condition with ON-OFF = 1, an identical amount of 
data similar to that specified in the card group 3 
must be given. Those numbers designate the tempera- 
ture magnitude applied to the corresponding element 
type and element number. This set of cards should 
be given for all load conditions with ONOFF = 1. 
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Card Group 5 


(7E10.4) 


Temperature change for each thermal load group must 
be given. ITHLDG real numbers are required. 

Comments : 

1. Each element is considered to have uniform tempera- 
ture. 

2. Temperature change should be computed with respect 
to an appropriate uniform reference temperature. 

Note that if all elements are made of the same 
material and assume the same temperature, then 
thermal stresses are not induced. 

XVI . Flight Condition Data 

This block of data will be reserved for future development 
of the ACCESS-3 program which may include aeroelastic 
constraints. Supply one blank card. 

XVII. Lumped Nodal Mass Data 

Card 1 (15) 

NMASS : Number of lumped nodal masses 

Card 2 - (NMASS + 1 ) (15, 5X, E10.4) 

N^ : Node number to which the mass is 

attached 

w N : Weight of the mass 

l 

Note that the magnitude w N must be 

l 

given in weight units, not in mass 
units . 

Supply one blank card if there is no lumped nodal mass. 
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XVIII. 


Stress Constraint Approximation Data (3E10.4) 

The automatic selection between zero and first order 
approximations for the stress constraints proceeds as 
follows. A retained potentially critical stress constraint 
will be replaced by its zero order approximation (using 
stress ratio) , rather than by its first order Taylor 
series expansion (linearization) , if the following test 
is satisfied; 


STR-GRD 

STR 


£ EPS 


0 £ EPS £ 1 


where STR is the constrained quantity and GRD is the 
relevant gradient component (see Section 3.7). Two 
limiting cases are instructive; 

if EPS = 0 ; all retained stress constraints are replaced 
by full linear constraints (first order 
approximation) . 

if EPS = 1 : all retained stress constraints are replaced 
by simple side constraints (zero order 
approximation) . 

Initially the tolerance EPS is set to be EPS-initial 
and, at the end of each design stage, EPS is updated by 


EPS = EPS x (EPS multiplier) 

Since EPS-multiplier is chosen to be less than 1, EPS 
is decreased stage by stage, which means that more and 
more stress constraints are linearized as the design 
proceeds (but simultaneously more and more stress 
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constraints are truncated; see card group XXI) . 

Three real numbers must be specified: 

EPS-initial : initial tolerance for zero/first 

order approximation 
EPS-min : lower limit of EPS 

EPS-multiplier : EPS modification multiplier. 

Supply one blank card if this capability is not used. 
XIX-XXII Constraint Control Data 

There are 4 types of constraints which can be specified. 
Each constraint type may have different truncation control, 
although the method used is identical for all types of behavior 
constraints. The truncation strategy is similar to the one 
used in ACCESS-1 (Refs. 1 and 2), but the sign convention 
defining the feasible region is reversed. 

If a q constraint function at a design a is evaluated 
as h (a) (see Section 3.6), h (a) is compared with a truncation 

q q 

boundary value (TBV) which is determined by: 

TBV = + {Min [h (a) - C] }x TRF + C 

q q 

where Min is applied to all q's in the constraint type. 

q 

Initially, TRF is set to be TRF-initial and at the end of 
each design stage, TRF is updated by 

TRF = TRF x (TRF multiplier) 

Since TRF-multiplier is chosen to be greater than 1, TRF is 
increased stage by stage, i.e. TBV is decreased, which means 
that more and more constraints are truncated as the design 
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proceeds (see Fig. 4) . 

It should be noted that the side constraints do not appear 
in this block of data, because there is no need to truncate 
them. The side constraint codes are specified in the element 
data . 


XIX. Displacement Constraint Control Data 
Card 1 (15) 

NDPC : Number of constrained displacement degrees 

of freedom 


Card 2 (5E10.4) 

TRF-initial 
TRF-max 
C-cutof f 
TRF-multiplier 
Min. Norm Ftr. 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
Minimum constraint normalization 
factor. Constraints are usually 
normalized by the absolute values 
of the limiting values. 


Card 3 - (NDPC+2 ) (315, 5X, 2E10.4) 

Node . : Node number associated with the 

i 

i-th displacement constraint 
Ixyz : Direction identifier 

0 = not used 

1 = X direction 

2 = Y direction 

3 = Z direction 
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Code 


: -1 = Lower bound only 

0 = No constraint 


1 = Upper bound only 

2 = Both 

Lower Bound : Lower bound of the displacement 

component 

Upper Bound : Upper bound of the displacement 

component 

XX. Slope/Relative Displacement Constraint Control Data 

This constraint type is restricted to place bounds on 
relative displacement components of two arbitrary nodes. In 
other words, the difference between Y-displacement components 

j_ r_ 4- 

of the L 1 " and U nodes may be bounded. But the difference 

th 

between the Z-displacement of L^ and X-displacement component 

4- H 

of U c node cannot be bounded. 

Card 1 (15) 


NSLC : Number of slope/relative-displacement 

constraints 


Card 2 (5E10.4) 

TRF-initial 
TRF-Max 
C-cutof f 
TRF -multiplier 
Min. Norm. Ftr 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
Minimum constraint normaliza- 
tion factor 
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Card 3- (NSLC+2) (315, E10.4) 


Node 


(L) 


Node ( “> 


xyz 


Upper Bound 


Note : 


Node number of the L node 

. th 


associated with the i 
constraint 


slope 


Node number of the u node 

*th 

associated with the 1 slope 

constraint 

Direction and code 

0 : not used 

1 : X direction \ 

2 : Y direction 

3 ; Z direction 

4 : X direction \ 

5 : Y direction 

6 : Z direction 

Upper bound of the slope/ 
rel. displ. 


i 




relative 

displacement 

slope 


1. If I =1, for example, the constraint function is 

XV Z 

, , rT Node ^ T7 Node 

1 - (U - U ) /Upper Bound i 0 

X X 

2. If I VV17 = 4, for example, constraint function is 

AX4I 


U 


Node 


(U) 


1 - 


X 


- u 


Node 


(L) 


X 


D 


/Upper-bound ^ 0 


YZ 


45 



3 . 


If lower bound is to be specified, node^ and 

node ^ should be exchanged to transform it to an 

upper bound constraint. 

XXI . Stress/Strain Constraint Data 

Card 1 (1015) 

MTYP 

Code : Stress/Strain constraint code (see 

Section 3.5) 

Except for element type 3 

-1 = read stress constraint code element by element 

0 = no stress constraint 

1 = all elements in this element type are con- 

strained by lower bounds on compression stress 

2 = all elements in this element type are con- 

strained by upper bounds on tensile stress or 
Von Mises combined stress (Element Type 1 or 
Types 2, 4, 5, 6) 

3 = effectively this implies that both codes 1 

and 2 are applied simultaneously 
For element type 3 

-1 = read strain constraint code element by element 

0 = no strain constraint imposed 

1 = maximum strain envelope criteria imposed on 

all elements 

2 = stress interaction criteria imposed on all 

elements 

3 = Tsai-Azzi criteria imposed on all elements 
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Card 2 (7E10.4) 


TRF-initial 

TRF-max 

C-cutof f 

TRF-multiplier 

Min. Stress Norm. Ftr. 


Min. Strain Norm. Ftr. 


TEBCF 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
Minimum stress constraint 
normalization factor 
Minimum strain constraint 
normalization factor 
Truss Euler buckling control factor 
If TEBCF > 0, TEBCF stands for the specified mean 
radius r of the truss element assuming tubular cross 
section. Stress constraint is 

2 Er 2 /2Z 2 } 

If TEBCF < 0, it stands for the thickness (t) to mean 
radius (r) ratio of the truss element assuming 
cylindrical cross section. Stress constraint is 

° 2 Maxto allowable , tt 2 EA/[4S. 2 - (^) ] > 


0 > Max{0 allowable , it 


If TEBCP = 0, no Euler buckling constraints are con- 


sidered . 

Card 3 (1415) 

Stress/strain constraint specification for element type 

MTYP 

code, Code MTYP <0. If all Code are positive, no 

cards are required. 

For each element type with Code MTYP = -1, stress/strain 


47 



code must be given to all elements sequentially start- 
ing from element number 1. 
Element_stress/strain_constraint_code : 

Stress code 

-1 : only compression stress is bounded 

0 : no constraint 

+1 : only tensile (truss only) or Von Mises 

combined stress is bounded 
+2 : both compressive and tensile stress 

are bounded. 

Strain Code 

MTYP 

same as Code specification 

XXII. Natural Frequency Constraint Data 
Card 1 (215) 

NFREQ : number of lowest frequencies to be bounded 

NSPACE : frequency constraint approximation scheme 

0 = first order Taylor series expansion with 

respect to linked reciprocal variables 
(linear in the optimization design space) . 

1 = first order Taylor series expansion with 

respect to linked direct variables 
(with NEWSUMT optimizer only: I0PT=4) 

2 = second order Taylor series expansion 

with respect to linked direct variables 
(with NEWSUMT optimizer only: I0PT=4) 
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Card 2 


(7E10.6) 


: Initial truncation factor 

: Upper limit of TRF 

: Cutoff base value 

: TRF modification multiplier 

: Minimum constraint normalization 

factor 

: Eigenvalue analysis convergence 

criteria (see note below) 

: Acceleration of gravity 

If 0.0, American standard unit 

is assumed and replaced by 
2 

386.0 in/sec . 

Note: Subspace iteration algorithm is used to obtain 

eigenvalues and eigenvectors. Iteration is judged 
to be converged if the relative differences of all 
eigenvalues are less than Eig. Conv. 

Card 3 (15, 2E10.4) 

Code^ ■*": constraint code 

-1 = lower bound only 

0 = not bounded 

1 = upper bound only 



2 = 

lower 

and upper 

bounds 



Lower 

Bound : 

lower 

bound on 

the 

. th 
i 

eigenvalue 

(«?) 

Upper 

Bound : 

upper 

bound on 

the 

. th 

l 

eigenvalue 

Uk) 


TRF-initial 
TRF-max 
C-cutof f 
TRF-multiplier 
Min. Norm. Ftr. 

Eig. Conv. 

Acc. Gravity 
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XXIII. Flutter Constraint Data 


XXIV. 


Supply one blank card. 

Optimizer Control Parameters 

Four distinct options are available for solution of the 
explicit approximate problem generated at each design stage 
(see Section 3.8). The following block of data takes on 
different meaning depending upon the value selected for 


IOPT (Card II) 


If IOPT = 4 : NEWSUMT Optimizer Control Cards 


Card 1: (515) 


JPRINT 


MAXSTG 

MAXRSF 


MAXODM 


JSIGNG 


Optimizer printout control (see Table 6) 
standard output = 0 
Maximum allowable number of stages 
Maximum number of response surfaces 
per stage; i.e. response factor is 
reduced MAXRSF times before the approxi- 
mate problem is updated 
Maximum allowable number of one dimen- 
sional minimizations per response surface 
sign of feasible region 
1 

-1 


feasible region is gg(a)>0 


feasible region is g^(a)iO 


Card 2: (8E10.4) 


EPSSTG : Stage convergence criterion. 

Overall iteration is judged to be con- 
verged if both of the following conditions 
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I W - W , I /W £ EPSSTG 
' P p-1 1 ' p 

' W P-1 " W p-2' /W p-l * EPSSTG 

EPSODM : Unconstrained minimization convergence 

criterion. Convergence is obtained if 
the relative values of total function at 
the ends of 3 successive one dimensional 
minimizations are not different by EPSODM. 
RACUT : Response factor decrease ratio 

RAMIN : Minimum response factor 

STEPMX : Maximum step size at each stage (move 

, limit) . All design variable components 

are constrained by 

STEPMX ^ ^i ^ STEPMX i - 1, . . .B. 

ITP : Initial transition point for the extended 

penalty function 
Power Fr : specify = 0.5 

Coefficient : specify =1.0 

Card 3 (2E10.4) 

STEPMX-mul : STEPMX modification multiplier 

STEPMX-£.£. : Lower limit on the STEPMX 

If IOPT = 5 : PRIMAL2 Optimizer Control Cards 
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Card 1 (315) 


JPRINT 


MAXSTG 

MAXODM 


Optimizer printout control (see Table 7) 
Standard output = 0 
Maximum allowable number of stages 
Maximum allowable number of one dimen- 
sional minimizations per stage; i.e. 
MAXODM search directions are computed 
before the approximate problem is updated. 

if MAXODM = 0, complete 
solution of the approximate problem is 
performed at each stage. 


Card 2 (5E10.4) 


EPSSTG 

DUMMY 

DUMMY 

DUMMY 

STEPMX 


: Stage convergence criterion. 

Same as NEWSUMT 
: not used 

: not used 

: not used 

: Maximum step size at each stage (move 

limit ) Same as NEWSUMT 


Card 3 (2E10.4) 

STEPMX-mul : STEPMX modification multiplier 

STEPMX- £ . H . : Lower limit on the STEPMX 


If IOPT = 6: DUAL1 Optimizer Control Cards 
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Card 1 (515) 


JPRINT 

MAXSTG 

MAXRES 

MAXODM 


I COMB 
Card 2 
EPSSTG 

EPSODM 

TAUMAX 


: Optimizer printout control (see Table 8) . 

Standard output = 0. 

: Maximum allowable number of stages 

: Maximum allowable number of restarts in 

the solution of the dual problem (dis- 
crete and mixed discrete-continuous cases) . 

: Maximum allowable number of one dimensional 

maximizations per restart. MAXODM should 
be at least equal to the number of con- 
straints retained in the first stage. 

: Specify = 1 

(5E10.4) 

: Stage convergence criterion 

Same as NEWSUMT 

: Dual maximization convergence criterion 

: Maximum step size in dual space. 

Standard_option : if TAUMAX = 0, the 
program automatically estimates an 
appropriate value for TAUMAX and simul- 
taneously determines a good starting 
point for dual maximization in the first 
stage. TAUMAX f 0 can be used in order 
to reduce the storage requirement for 
computation of the intercept-distances 
to discontinuity planes 
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DUMMY 


not used 


STEPMX . Maximum step size at each stage (move 

limit in the primal space) . 

Card 3 (2E10.4) 

STEPMX-mul : STEPMX modification multiplier 

STEPMX-il.il . : Lower limit on the STEPMX 


If IOPT = 7 : Dual 2 Optimizer Control Cards 


Card 1 (215) 
JPRINT 

MAXSTG 

Card 2 (5E10.4) 
EPSSTG 

EPSODM 


DUMMY 

DUMMY 

STEPMX 

Card 3 (2E10.4) 
STEPMX-mul 
STEPMX- A . &. 


Optimizer printout control (see Table 9) 
Standard output = 0. 

Maximum allowable number of stages 

Stage convergence criterion 
Same as NEWSUMT 

Dual maximization convergence criterion. 
Convergence is achieved when the norm 
of the vector made up of the primal 
constraint values is less than EPSODM. 
not used 
not used 

Maximum step size at each stage (move 
limit in the primal space) . 

STEPMX multiplier 
Lower limit on STEPMX. 
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Figure 1. Basic Organization of ACCESS 3 
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Figure 3. SSP Element Model Example 
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Table 1. Temporary Files 


File Name 


Contents 


10 


11 

12 

13 

14 

15 

16 

17 


18 

19 


20 


21 


22 


Stiffness matrix components associated with 
unit values of independent design variables 
and load vector components which are indepen- 
dent of design variables 

Mass matrix components associated with unit 
values of independent design variables 

Load vector components due to thermal loads 
and dependent on independent design variables 

Load vector components due to inertia loads 
and dependent on independent design variables 

Constraint gradients 

Input data and a part of preprocessor output 

Thermal shear panel stiffness matrix compo- 
nents. Used only when IETP(6)^0 

Constraint gradients (except stress constraints) 
Required only when zero order approximation 
of the stress constraints is considered. 

Eigenvector sensitivity vectors, if computed 

Mass matrix post-multiplied by eigenvectors. 
Required only when second order expansion of 
frequency constraints is used. 

Original system stiffness and mass matrices. 
Required only when second order expansion of 
frequency constraints is used. 

Modified [K-XiM] in the eigenvector sensitivity 
computation. Stored in decomposed form. 

[H b - x i ' <<i=l, NEIG), b=l, B) 

Required only when second order expansion of 
frequency constraints is used. 
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Table 2 . Required Blank Common Size 


Problems 

Elements 

Total No. 

Free Displ. 

No. of 

Total No . of 

Required 



of Elements 

d.o.f .s. 

Design Variables 

Constraints 

Real Array 

Integer Array 

Wing Carry-Through 








Truss Model 

TRUSS 

63 

42 

63 

256 

4291 

1606 

(Static) 








Swept Wing (Metal ) 

CST 

60 

120 

18 

268 

5238 

2305 

(Static) 

SSP 

70 






Delta Wing-Composite 
(Static & Dynamic) 

CSTOR 

252 

70 

105 

60 

2661 

13090 

9433 

Delta Wing- Composite 

CSTOR 

252 






(Static, Thermal 

SSP 

70 

105 

60 

2661 

18652 

9103 

and Dynamic) 

TSP 

70 






Delta Wing-Composite 

CSTOR 

252 






( Static , Thermal , 

SSP 

70 

105 

g 0 ^12 continuous 
'"48 discrete 

2661 

22327 

8567 

Dynamic , discrete 

TSP 

70 





variables) 












Table 3. Available Options for Frequency Constraints 


A=o) 2 

Approx. 

Algorithm 

1st order 
Reciprocal 
DV's 

1st order 
Direct 
DV's 

2nd order 
Direct 
DV's 

NEWSUMT 

* 

* 

* 

PRIMAL2 

* 


- 

DUAL1 

* 

- 

: 

DUAL2 

* 

- 

- 


*available combination in ACCESS-3 program 


Table 4. Algorithm Options for Various Kinds of Problems 


Kinds 
of DV's 

Algorithm 

Pure 

Continuous 

Pure 

Discrete 

Mixed 

Continuous - 
Discrete 

NEWSUMT 

* 

- 


PRIMAL2 

* 

- 

- 

DUAL1 

* 

* 

* 

DUAL2 

* 

- 

- 


* available for application in ACCESS- 3 program 
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Table 5 Analysis Printout Control - IPRINT 


All messages above the horizontal line corresponding to each 
value of IPRINT are printed ' 


IPRINT 

Information Printed 


Constraint identification code 
Posture table at each stage 
Time statistics of the job 

o 

Messages prior to any error termination 


Input data in readable format 
Initial and final nodal displacements 


Initial and final eigenanalysis results 


Element sizes and weight information 
Scaling factor and scaled weight 
New list of linearized constraints (after 
zero order stress approximation) 
Independent linked design variable values 


(including lower and upper limits) 


Modified truncation factors 
Initial and final element stresses 

9 * 

Initial and final values of all constraints 


Element stresses and constraint values at 
each stage 

Lower limits, upper limits and allowable 


discrete values after design variable 
normalization 

Q 

Interface data (with optimizer) 

0 

Listing of data cards 
Element geometry data 
Load vectors 

! A 

Gradients of retained constraints 

4 

Element stiffness/mass matrices 
Master stiffness/mass matrices 
Integer and real pointer vectors for 
dynamic array allocation 
Debugging of integer and real arrays 


* Standard Values 



Table 6 NEWSUMT Optimizer Printout Control - JPRINT 


All messages above the horizontal line corresponding to each 
value of JPRINT are printed 


JPRINT 

Information Printed 


Control and system parameters 
Initial design analysis results 
Response surface penalty multipliers 
ODM's results summary 
Final results of optimization 

n* 

Time and counting statistics 

1 

Results for each response surface 


Direction finding data 

4> 

ODM's results at each design point 

O 

Golden section search data 


Penalty function detailed data 


* Standard value 
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Table 7 PRIMAL2 Optimizer Printout Control - JPRINT 


All messages above the horizontal line corresponding to each 
value of JPRINT are printed 


JPRINT 

Information Printed 


Initial and final results summary 

u 

Final values of all linear constraints 
and associated dual variables 
Identification of each constraint added 


to or dropped from the active set. 

_L 

Final design variables 
ODM's results summary 


Active set strategy data 

Z 

Design variables and search directions 


* Standard value 
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Table 8 DUALl Optimizer Printout Control - JPRINT 


All messages above the horizontal line corresponding to each 
value of JPRINT are printed 


JPRINT 

Information Printed 

0* 

Control and system parameters 
Final results summary 

Final dual variable values 
Summary of dual solution analysis 
(discrete case) 

Results in brief summary form (each 




restart) 

1 

Initial dual variable values (each 
restart) 

Final primal variable and constraint 


values 


Analysis of upper and lower bound 


solutions (discrete case) 


Results in brief summary form (each ODM) 


Initial primal variable and constraint 
values (each restart) 

List of best discrete solutions 

3 

Primal and dual variable values (each 


ODM) 

Search direction and sensitivity of 


primal variables in ODM 


Newton iteration results 

4 

Detailed ODM's data 

5 

Detailed Newton iterations data 


* Standard value 
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Table 9 DUAL2 Optimizer Printout Control - JPRINT 

All messages above the horizontal line corresponding to each 
value of JPRINT are printed 


JPRINT 

Information Printed 

0* 

Final results summary 

1 

Identification of the starting point 
(first stage only) 

Final values of all dual variables and 
associated primal constraints 
Identification of the dual subspace 
at each iteration 

X 

9 

Final primal variable values 

Summary of results after each iteration 

9 

Data on search of a suitable dual start- 
ing point (first stage only) 

Search direction and dual variables at 
each iteration 
ODM's results 


Primal variable values after each 
iteration 


* Standard value 
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Fig. A-l Space Truss Element 


Strain-Displacement Relation (local coordinate) 


‘ " \l 



Stress Strain Relation (local coordinate) 

a = E £ - EaAT ( A— 2 ) 

where a : thermal expansion coefficient 

AT : average temperature change 


Force Displacement Relation (local coordinate) 


EA 

L 


1 

-1 


u ) 

-1 

1 


- 

U 0 I 




Q/ 


EaAT A 




= 0 


(A- 3) 


where F , F 

p q 

nodes , 


are externally 
respectively. 


applied force at P and Q 


Force Displacement Relation (reference coordinates) 
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Consistent Mass Matrix (reference coordinates) 



(A-5) 


where p : density 
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Stress-Strain Relation (local coordinate) 



Stress-Displacement Relation 



(local coordinate) 



EaAT 

1-v 



(A- 8) 


Local-Reference Displacement Relation 
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where X 


x 


unit vector parallel to the x-axis 


X : unit vector parallel to the y-axis 

Up , U q , U p : displacement vectors of P,Q,R nodes. 


[T] = 


Jl m n 
XXX 


Z m n 

y y y 


Z m n 

XXX 


Z m n 

y y y 


Z m n 

XXX 

l y m y "y 


{A— 10) 


Stiffness Matrix (local coordinate system) 
K = K + K 


n 


(A-ll) 


where 


K = 
n 


Et 


4A(l-v ) 


(s-b) -v(s-b)h -(s-b)s v(s-b)h (s-b)b 


vhs 

2 


Symm. 


-h" 

-vhs 

,2 


-vhb 

-bs 

vbh 
, 2 


0 

0 

0 

0 

0 
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where p : density 
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hs 


0 -bh 


-(s-b)s 0 (s-b)b 

-hs 0 bh 

s 2 0 -bs 

0 0 
b 2 

✓ 

1 coordinates) 

= f (A-12) 


0 0 
1 0 
0 1 
0 0 
1 0 
0 1 
0 0 
2 0 


(A-13) 


2 





Stress-Strain Relation (material axis) 



1-v v 
LT TL 


V LT E T 

1-V LT V TL 


v 

r=v 


tl e l 

LT V TL 


1-V LT V TL 



E L (a L +V TL a T } 

1-V LT V TL 


^ E T( a T +V LT a L) 
1-V LT V TL 



Strain Transformation Law (material-local) 



Lx 

Tx 


'Ly 

( 2 

'Ty 


^Lx^Ly 


^Tx^Ty 


2Z T £_ 2l r Z m l- Z m + Z m Z T 
Lx Tx Ly Tx Ly Tx Ty Lx 


[T] 



where 


Note : 


Stress 



z r 

-T 

= e_ • 

— *• 

e = 

COS0 , 1- 

-T -* 

= e_ * e = -sm 

Lx 

L 

X 

Tx 

T x 

&_ 

—T 

= e T * 

e = 

sin0 , 

-*T 

= e_ • e = cos0 

Ly 

L 

y 

Ty 

T y 

the 

direction of 

IT., is chosen 

so that 

( ®1 

X ef 2 ) 


x S'y) > 0. 



-Displacement Relation 
| = [D] [T] [B] u - AT • £ 


(A-16) 


(A-17) 
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Local-Reference Displacement Relation 


same as type 2 

Stiffness Matrix (local coordinate) 

K = K n + K s (A-18) 


(^(s-b) -C 2 (s-b)h -C 1 (s-b)s C 2 (s-b)h C 1 (s-b)b 0 



C 2 h 2 C 2 hs 

-Cjh 2 

-c 2 bh 

0 

K n 2bh(l-v LT v TL > 

C x s 2 

-C 2 hs 

-C^bs 

0 


Symm 

°2 h2 

c 2 bh 

0 




C 1 b 2 

0 





0 

where 

S 


(A-19 ) 

y 

c i = p e l + 

2 p 2 u 2 v lt e t + u 4 e t 




2 2 

C 2 = P ^ E L 

* (p 4 + M 4 )v lt E t + P 2 P 2 E t 





C 3 = y 4 E L + 2p 2 p 2 v LT E T + p 4 E T 


v lt e t v tl e l 


p = sin0 
U = cos 9 
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hfs-bjD.,^ 

h 2 D 1 -2h(s-b)D 2 

+ [h 2 - (s-b) 2 ] D 2 -h (s-b) d 3 

+ (s-b) 2 D 3 

-s (s-b) 

-hsD^ [ (s-b) s-h 2 ] D 2 

-(2s-b)hD 2 -h 2 D 3 

+(s-b)h D 3 

-h (s-b) 

-h 2 D 1 +h(2s-b)D 2 

+ [ (s-b) s-h 2 ] D 2 +sh D 3 

-s (s-b) D 3 

b(s-b)D^+bh D 2 

bhD. -b (s-b) D„ 

' 1 2 


b(s-b)D 2 -bh D 3 -bhD 2 +b (s-b) D 3 

2 2 

where = 4p p 

D 2 = 2pp (p 2 -y 2 ) 















Equilibrium Equation (local coordinate) 


K u + 


h = 1 


h = t AT <{ 


-(b-s) (p 2 h 1 +u 2 h 2 ) + 2hpM(h 1 -h 2 ) 
-hfM^ + p 2 h 2 ) + 2(b-s)py(h 1 -h 2 ) 
-sfp 2 ^ + y 2 h 2 ) -2hpy (h 1 ~h 2 ) 
h(y 2 h 1 + P 2 h 2 ) + 2spy(h 1 -h 2 ) 
b(p 2 h 1 + u 2 h 2 ) 


> 


-2bpy (h 1 ~h 2 ) 


(A-21) 


where 


„ E L (a L + 'WV 

1 t-Wn, 


h 2 


E T (C1 T + Wl, 1 
1 ~ V LT U TL 


Consistent Mass Matrix (reference coordinate) 


same as type 2 
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4. Type 4 - SSP: Symmetric shear panel element with uniform 

thickness and isotropic material 

This is a special element used to model relatively thin 
symmetric structures such as idealized supersonic lifting 
surfaces. Theoretical discussion is given in Ref. 1. It 
is assumed that this element models the upper (or lower but 
not both) half of the symmetric structure and the element 
plane of symmetry coincides with the X-Y plane. It is 
further assumed that all SSP elements are placed vertically 
with respect to the X-Y reference coordinate plane. 
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3 


. If the heights PP" and QQ' are different, the average 
(PP* + QQ')/ 2 is considered as the height of the 
element, i.e. b. 

4. No thermal load can be considered in this element 


Strain Displacement Relation (local coordinate) 



Stress-Strain Relation (local coordinate) 



(A-22) 


(A-23) 


Stress Displacement Relation (local coordinate) 


= E 


xy ' 


' _ 2n 
a 

0 

2 _n 

a 

v 

0 

| 

5P ) 

0 

0 

0 

0 

< 

l*P ( 

1 

1 

1 

1 


*Q 1 

,2 (1+v) b 

2 (1+v) a 

2 (1+v) b 

"2 (1+v) a 


/ 


(A-24) 
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Local to Reference Displacement Transformation 

U I 
V I 
W I 
0. 




l 

X 

m 

X 

0 

V 





0 


V 


0 

0 

±1* 


K 1 



0 

J l m 

X X 

0 

Vfy 




0 0 

±1* 







*(- 

) sign 

if Z P 

< 0 and Zq < 0 



(A-25) 


V, 


w. 


where SL „ and m are components of a unit vector e" along 
X X X 

the local x axis. 

Stiffness Matrix (local coordinate) 

F+3a -3 -F+3a 3 

3/a -3 -3/a 

F+3a 3 


K = 


Et 

12 (1+v) 


Symm. 


3/a 


(A— 26 ) 


where a = 


a 

b 


F = 


2 (1+v) 


a 
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Consistent Mass Matrix (local coordinate) 


M = 

w 2 


J +G 


12 

1 

3 


H 


Symm. 


0 3J“ G 


l 

9 


12 


H 


j + s 


IT H 


1 

6 


12 


1 

3 


18 


H =■ 


1 

9 


(A-27) 


where p = density 


2 2 
r = + JL + -V 

30 18 30a 2 


H = 


a + i 
a 


It may look strange that the mass matrix depends upon 
Poisson's ratio v through G. This is due to the fact that 
the assumed displacement field is derived based on assumed 
stress field. 
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5 . 


Type 5 - PSP: Pure symmetric shear panel element with uniform 
thickness and isotropic material 


This element is identical to a type 4 (SSP) element, 

except for a minor change in the assumed displacement state 

so that the stress state of the element is pure shear: i.e. 

a = a =0. This implies that e = e so. 
x y e x y 


Strain Displacement Relation (local coordinate) 

,T 


'xy 


[i - ii il (u 

_ b , a, b, a J \ 


p, v p , u Q , v 


Stress-Strain Relation (local coordinate) 


xy 


2 (1+v) Y 


xy 


Stress-Displacement Relation 
E 


1 1 1 ll ~ ~ . ~ •> 

b' - a' b' al { U P' V P' U Q' V Q } 


xy 2 ( 1+v) 

Local to Reference Displacement Transformation 

same as type 4 . 

Stiffness Matrix (local coordinate) 
same as type 4 except F e 0. 

Mass Matrix 


(A-28) 


(A-29 ) 


(A-30) 


Assumed to be the same as type 4, 
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Strain-Displacement Relation (local coordinates) 




Stress-Displacement Relation 



same as type 4 


\ 

i 
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Stiffness Matrix (local coordinate 



where a = §■ 
b 

Force Displacement Relation (local coordinate) 



(A-34) 


(A-35) 


Consistent Mass Matrix 

Assumed to be the same as type 4 

Notes As shown in the stress-displacement relation, stress 
distribution is linear with respect to x. In order 
to simplify the problem, an approximate stress dis- 
placement relation is used in computing stress and 
stress sensitivity. 
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(A— 36) 
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6 . 


Type 6 - TSP: Thermal symmetric shear panel element with 


uniform thickness and isotropic material 

Since SSP and PSP cannot be used for problems involving 
thermal loads, this special element is added to the ACCESS - 3 
element library. The TSP element is designed to be used under 
steady thermal soak load conditions such that the temperature 
change in each TSP element is uniform and therefore symmetric 
with respect to the X-Y plane. 

If the -structure is subject to both mechanical and thermal 
loads, two structural models must be created and analyzed separ- 
ately. One model is to use SSP elements to model shear panels 
and it is subject to only mechanical loads. The other model uses 
TSP elements to model the shear panels and it is subject to only 
thermal soak loads. These two models are created automatically, 
if the user specifies both SSP and TSP elements. Displacement 
and stress states of the structure subject to both thermal and 
mechanical loads are generated by superimposing the results ob- 
tained from the two separate models. 

Theoretically, it is also possible to consider the PSP - TSP 
element combination, but this is not implemented in the current 
version of the program. 

Note that the TSP option requires a significant amount of 
core memory and CPU time, since two system stiffness matrices 
are stored and decomposed. Sensitivity analyses of the responses 
must be carried out separately and superimposed afterwards. 
Therefore, analysis effort is nearly doubled when thermal effects 
need to be considered. 
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Appendix B 


Examples 

Two simple examples are given to illustrate input data preparation 
and program output for various features of the ACCESS 3 code . To help 
understand these examples , Figs . B1 and B2 represent the geometrical layout 
of the structures and the following indications are provided: 

(1) 10-bar cantilever truss 

. static constraints only; 

. automatic selection of stress constraints requiring first order 
approximation ; 

. equality constraints on displacement; 

. DUAL 2 optimizer 

Note that displacements at nodes 4 and 5 in the Y direction are 
required to be equal to -5.08 cm (-2.0 in) and -2.54 cm (-1.0 in), 
respectively. 

(2) 10 element delta wing 

. static constraints and frequency constraint; 

. mechanical and thermal loads; 

. mixed continuous-discrete problem; 

. titanium webs (continuous) and composite skin (discrete); 

. DUAL 1 optimizer. 

Note that the available discrete thicknesses for the CSTOR elements 
are {0.0254, 0.0508, 0.0762, 0 .1016 , . . . ,5 . 08} (cm) or {o.Ol, 0.02, 
0.03, 0.04 , . . . .2 .0} (in). 
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Fig. B1 10 bar cantilever truss 



Fig. B2 10 element delta wing 
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(1) 10 bar cantilever truss 
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COMPLETE OPTIMZaTICN BY DUAL2 OPT IM IZER 


ANALYSIS PRINT OUT CONTROL = 2 


DESIGN 


IN LINKED SIZING VARIABLE SPACE 


NUMBER 

NUMBER 

NUMBER 

number 

NUMBER 

NUMBER 

NUMBER 


OF NODES 

OF TOTAL ELEMENTS 
OF LINKED VARIABLES 
CF LC AC CONDITIONS 
OF BOUNDARY NODES 
CF ISOTROPIC MATERIALS 
CF ORTHCTROPIC materials 


6 

10 

10 

1 

6 

1 

0 


TRUSS 


CST 

ISOTROPIC 


CST 

ORTHOTROPIC 


SSP 


PSP 


TSP TUD TBD TBD TBD 


ELEMENTS 10 

LINKED VARIABLES 10 


NODE NUMBER 

X 

Y 

z 

1 

0.0 

0 . 0 

0.0 

2 

360.0000 

0.0 

0.0 

3 

720.0000 

0.0 

0.0 

4 

720.0000 

360. 0000 

0.0 

5 

360 .0000 

360 . OCOO 

0.0 

6 

0.0 

360.0000 

0.0 


VO 

UJ 


DISPLACEMENT dOUNUARY CONDITIONS 


BOUNDARY CODES* PRESCRIBED DISPLACEMENT 


NODE NO. 

X 

Y Z 

X 

Y 

z 

1 

i 

1 1 

0.0 

0 .0 

0.0 

. ... 2 _ 

0 

0 1 

0.0 

0.0 

0.0 

3 

0 

0 l 

0.0 

0.0 

0.0 

4 

Q 

0 1 

0.0 

0 .0 

o.a 

s 

a 

0 1 

0.0 

0.0 

0.0 

0 

1 

1 1 

0.0 

0.0 

Q.Q . 


* — 1=PR ESCR I BED « 0=FREE. 1=F IXED 



element 


NCUE NUMBERS 

LINKED 


ELEMENT SIZE 


MATERIAL 

SIDE CONSTRAINT 


NO. 

N1 

N2 N3 N4 

GROUP 

INITIAL 

LOW.dD 

UPP-.BL) 

GROUP 

cooe* 



TRUSS 

elements 









l 

1 

2 

1 

2 0. C 00000 

0. 100000 

0.100000 


— 1 

0 


2 

2 

3 

2 

20.000000 

0. 100000 

0.100000 


— 1 

0 


3 

3 

- _4 

2 

20.000000 

.. 0.100000 . 

Q.IOOQOQ 



0 


4 

4 

5 

4 

20.000000 

o . loooao 

0. 100000 

1 

«-t 

0 


s 

5 

6 

5 

20.000000 

o. looaoo 

0.100000 

1 

—1 

0 


6 

2 

5 

6 

20.000000 

0 .10000 0 

0.100000 

I 

— 1 

0 


z 

l_ 

.5 

7 

20.000000 

. . 0 * 1 00 00 Q 

g.taoooo 

1 

• 1 

0 


8 

2 

6 

a 

20.000000 

0 . 100000 

0.100000 

1 

• 1 

0 


0 

2 

4 

S 

20.000000 

0.100000 

0.100000 

1 

— 1 

0 


10 

3 

5 

1 0 

20. 000 JOO 

0. 100000 

0.100000 

1 

— l 

0 


* - 

2=F I XEQ AT INITIAL VALUE 

-1 =LOM£R 

BOUNDS CNLY 

.... - - 



— 

- 


— 


0 = NCN NEGATIVITY ONLY 1=UPPER BOUNDS CNLY 

2 = d CT H UPPER AND LOAER BOUNDS 

■&EC ON Q-hU-MJjER_=. -DISCRETE- VALUE fiRQUP NUMBER 


MATERIAL CONSTANTS 


isotropic materials 



tLhMt Nr 
■MO. 

N1 

NClJE NUMJLPS 
N2 NJ 

N4 

L INKED 
GROUP 

INITIAL 

element 

low . a a 

SIZE 

UPP.tJD 

MATERI AL 
GROUP 

SIDE CONSTRAINT 
CODE* 

r H JS5 

1 

hLLVtNTS 

l 2 


1 

20 • CQOOOU 

a . looooo 

0.100000 

1 


0 

2 

2 

3 


2 

20.000000 

0 • 100000 

0. 100000 

1 

— | 

0 

3 

3 . 

4 . 


3 

20.000000 

o . looaoo 

0.100000 

1 


0 

4 

4 

5 


4 

20.000000 

0 . 100000 

0 . 100000 

1 

— | 

0 

5 

5 

6 


5 

20.000000 

0 . 100000 

0. 100000 

1 

— 1 

0 

6 

2 

5 


6 

20.000000 

0.100000 

0.100000 

1 

— i 

0 

7 

.1 

.. 5 . 


7 

20.000000 

0. 100000 

0 .100000 


— | 

0 

a 

2 

6 


6 

20.000000 

0 . 100000 

o. looooo 

1 

— 1 

0 

9 

2 

4 


9 

20.000000 

0.100000 

0.100000 

I 

— 1 

0 

10 

3 

5 


i o 

20.000000 

0 . 1 JOOOO 

0.100000 

i 

“ I 

0 

* 

-2=FIXEQ AT INITIAL 

VALUE 

-1 =LOWER 

BOUNDS CNLY 


.... ■ • — • 



■ — ~ — 


O = NCN NEGATIVITY ONLY 1 =UPPER BOUNDS CNLY 

2=dCTH UPPER AND LOWER dCUNCS 

SECOND N UMdER = DISCRETE VALUE GROUP NUMBER 


MATERIAL CONSTANTS - ISOTROPIC MATERIALS 


GROUP YOUNG'S POISSON'S 

NO. MOOULUS RATIO 

SPECIFIC 

WEIGHT 

THERMAL 

EXPANSION. 

COMPRESSIVE 
A. STRESS 

TENSILE 
A. STRESS 

1 IOCCOCOO.O 0.3000 

0.100000 

1 

1? 
1 o 

—25000.0 

25000.0 


LOAD CONDITIONS 


LUMPED LOAD AT NODES 


NODE NUMBER X Y 2 


LOAD CONDITION 1 _ ; VO 

"2 0.0 -lOOOOC.OCO ' 0.0 

3 0.0 -lOOOOC.OOO 0.0 


PRESSURE LOAD 

NC PRESSURE LOAD SPECIFIED 

GRAVITY LOAD . . 

LOAD CONDITION NO. MAGNITUDE(G) DIRECTION COSINES 

— ■ 0 ' NO GRAVITY LOAD 


thermal load 

NO THERMAL LOAD SPECIFIED 


< 










CONSTRAINT data 


SIDE CONSTRAINTS 


SIDE CONSTRAINT SPEC I F IC AT IONS ARE GIVEN IN ThE ELEMENT DATA 


DISPLACEMENT CONSTRAINTS 


INITIAL T PONCA T I ON FACTOR 
MAXIMUM TRUNCATION FACTOR 
ft A SI S CUTOFF FACTOR 
MULTIPLIER FCR TRF UPDATING 
MINIMUM NORMALIZATION FACTOR 


NODE NO. 

UI SECT IDN** 

CODE * REGION NO. 

LOWER BOUND 

UPPER BOUND 

1 4 

2 

2 

-2. 000000 

-2.000000 

2 5 

2 

2 

-1 *000000 

-l.COOOOO 


♦ — 1 “LOW ER BOUND ONLY. Q=NO CONSTRAINT 

1 =UF?£R ECUND ONLY. 2=L0WER AND UPPER BOUNDS 
** IF Of SPHERICAL DISPLACEMENT BOUND 


0*1000 
0*8000 
1*0000 
1*2000 
0.0 100 


SLOPE/RELATIVE DISPLACEMENT CONSTRAINTS 


NO SLOPE CONSTRAINTS 


STRESS/STRAIN constraints 


INITIAL TRUNCATION FACTOR 
MAXIMUM TRUNCATION FACTOR 
BASIS CUTCFF FACTOR 
MULTIPLIER FOR TRF UPDATING 
MINIMUM NORMALIZATION FACTOR 
STRESS 

STRAIN 0. 

NO EULER BUCKLING CONSTRAINTS IMPOSED 

SELECTION FACTOR FOR . QaXQflQ 

ZERO ORDER APPROXIMATION 

cLEMtNT TYPE 1 

ALL ELEMENTS ARE CONSTRAINED EY BOTH LOWER AND UPPER BOUNDS 


Q«1 000 
0*8000 
I *0000 
1 *2000 


1000 *0000 
i COOOOE+OO 


-Sp- 

un 


OPTIMIZER CONTROL PARAMETERS 


CUAL 2 

PRINT CUT CONTFCL 
MAX. NC. CF STAGES 

DIMINISHING RETURN CRITERION AMCNG STAGES 
DIMINISHING RETURN CRITERION AMCNG 0*D*M* 
MAX. STEP SIZE ALLOWED IN A SINGLE STAGE 
STEP SIZE MODIFICATION FACTOR 
5 TE 3 SIZE MINIMUM ALLOWABLE 


I 

1 0 

0. 10G0E-03 
0*10 QOE — 0 3 
0* 10 OOE tO 4 
0*1 OOOEfO 1 
0*1 20OE+O 1 


CONSTRAINT IDENTIFICATION CCOfcS 


CONSTRAINT TYPE 
— 10 00 1 

l 

-10002 

io constraints in 

— 1 00 0 J -13004 

Tl-IS 

TYPE 

10005 

— 1 0 0 CO 

-10007 

-10003 

-I00Q9 

-IOQIQ 

CONSTRAINT type 
- l 001 20004 

2 

100120004 

4 constraints in 

-100220005 100220005 

This 

TYPE 






C CN STRAIN T TY^E 

3 

o constraints in 

This 

type 





- 

CONSTRAINT TYPE 
- 1 QQ 10 03 1 
-10010006 

I 001 00 J 1 
1 0010006 

20 constraints in 
-10010002 10010002 

-l 001 000 7 1 001000 7 

THIS TYPE 
-100 10003 
-10010008 

1 001000J 
1001 0006 

-I 00 10004 
-I 0010009 

10010004 

10010009 

-10010005 

-10010010 

1 0010005 
1 0010010 

CONSTRAINT TYPE 

s 

o constraints in 

This 

type 






constraint type 

6 

o constraints in 

This 

type 








CONSTRAINT IDENTIFICATION CCDtS 


CONSTRAINT TYPE 
... _ -IQ QO 1 

1 

-IQQQ2 

10 CONSTRAINTS IN 
-10003 -10004 

Tt- IS 

type 
1 000 5 

CONSTRAINT TYPE 
-l 001 20004 

2 

1 00 120004 

4 CCNSTRAINTS IN 
-100220005 100220005 

THIS 

type 

CONSTRAINT TY^tf 

3 

0 CCNSTRAINTS 

IN 

TFIS 

type 

CONSTRAINT TYPE 

... -ig.fi.lQ.QsU 

-10010006 

** 

J 00 1 000 l 
1 0010006 

20 CCNSTRAINTS IN 
-10010002 1 0C10002 

-10010007 10010007 

TUS TYPE 
-10010003 
-10010000 

CONSTRAINT TYPE 

5 

0 CCNSTRAINTS 

IN 

THIS 

type 

CONSTRAINT TYPE 

6 

0 CCNSTRAINTS 

IN 

This 

TYPE 

HUMi3ER . Of . JSEtlAV IOOR CONSTRAINTS 24 




NUMBER OF DESIGN 

VARIABLES 

1 0 





NUMd£9 3£ 0 iSCR£ r £ VAR I Afl.U E S 0 


-1O0C6 

-10007 

-1 0008 

-10009 

-laoio 

10010003 

lootoooa 

-10010004 
-l 0010009 

10010004 

10010009 

-10010005 

-10010010 

1 001000S 
1 0010010 


END OVERLAY PREPCO 


ENTER OVERLAY DESIGN 


VO 

Ov- ... 





STAGE NU 


l 


APPROXIMATE PROdLEM GENERATOR 


CURRENT VEMdER SIZE 


ME^dER TYPE NL V QEA L 

O. 20006*02 0420006+02 0420006+02 0420006+02 0420006+02 042000E+02 0.2000E+02 0.2000E+02 0.2000E+02 0.2000E+02 


CURRENT WEIGHT OATA 

MEMBER TYPE NUMBER 1 WtIGHT = 0 .B3il 292 E + 04 


VARIABLE STRUCTURAL WEIGHT 0.8392926+04 

FIXED 'TRLCTURAL WEIGHT 0.0 

TGTAL STRUCTURAL WEIGHT O . 83 9292 E+ 04 

NtHTSTSUCTURAL WEIQHTHT 0.0 . . .. _ 

TOTAL WEIGHT 0.839292E+04 

fit-ECK STAGE NC.= l 0. U9IE + 27 O.IOOOE + Ol MLST BE LESS THAN 0.100000Er<>3 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0.200000E+31 0.100000E+31 0.839292E+04 


ENTER OVERLAY ANALYS 



NOCAL DISPLACEMENTS 


NODE . _ . Jt ...... Y i NODE. X ... X.. ..... Z 


LOAD CCNDITICN l 

l 0.0 0.0 0.0 2 -C.36834E+00 -0.90I05E+00 0.0 

3 —0.4/61 IE+OC — 0.19698E + 01 0.0 4 0.42388E + 00 -O.I8975E+OI 0.0 

S 0.351 65E + 0 0 -0.83/1 76+00 0.0 6 0.0 0.0 0.0 


NQCAL DISPLACEMENTS 


NQQE. . X 


Y 


node x r 


z 


LOAD CCNDITICN l 


l 0.0 0.0 0.0 
3 -0. 4761 IEEOO -0.19698E+01 0.0 
5 0.35165EF00 -0.83M7E*00 0.0 


2 

4 

6 


-0.36834E400 
o.423aae*oo 
o. o 


-0.90I03E400 0.0 

—0.1 Q975E40 1 0.0 

0.0 0.0 


m — as33ssssZ3a333Xn»3X ixxxiixaxxxixxxsixixxxx 


4 . „ DISPLAC EMENT C ONSTRAINTS _ 1 TO 4 MOST CRITICAL CONSTRAINT* -0.162831 1 E*00 

0.5123E-OI -0.5123E-01 0.1628EF00 -0. I 62 8E-E0 0 


MTYP 

M 

LC _ 

S-CCMBINED 

SX 


SY 

1 

1 

1 

-0. 10231 6E+05 

-0 . 1023EE05 

c.o 

0.0 

1 

1 

1 

-0. 102316E+05 

-0.1023E+05 

0.0 

0.0 

. _1_ 

-2_ 

_J 

-0.29937^E+Q4 

-0.2994ET04 

0.0 

0. 0 

1 

2 

1 

-0.299372EE04 

-0 .2994E+04 

0.0 

0.0 

1 

3 

1 

0.20081 9E 404 

0.2006EF04 

0*0 

0.0 

1 

3 

l 

0.200619EE04 

0 .2006EF04 

0. ◦ 

0.0 

.. 1_ 

4... . 

1 . 

Q .200620E+04 

0.2006EE04 

0.0 

0.0 

1 

4 

1 

0. 20C62 0E+-0 4 

0.2006EE04 

0.0 

0.0 

1 

5 

1 

0. 9768 16E+04 

0.9768E+04 

0.0 

0.0 

1 

5 

1 

0.976Q 16E + 04 

0 •9768E4-04 

0*0 

0*0 

1 

6. 

_J _ 

Q._l 7 2446EE04 

0.1774E+04 

0. 0 

0. 0 

1 

6 

1 

0.177446E+04 

0. 1774EE04 

0.0 

0. 0 

1 

7 

1 

-0.674326E+04 

-0 .6743E + 04 

0.0 

0.0 

1 

7 

1 

“0. 674326E+04 

—0 .6743EF04 

0. 0 

0. 0 

X- 

.. _ a __ 


0.739875E+04 

0 .7399E + 04 

0.0 

0. 0 

1 

a 

1 

0. 73987 5fc *04 

0.7399EE04 

0.0 

0.0 

1 

9 

1 

-0.28272 IE+04 

-0 . 2337EE04 

0.0 

0. 0 

l 

9 

1 

-0 .283721 E + 04 

-0 .2837EE04 

0.0 

0.0 

1_ 

10 

1 .... 

0.423376EE04 

0 . 42J4EE04 

0.0 

0.0 

1 

1 0 

1 

0.423376EE04 

0 .4234E+04 

0.0 

0.0 


SX-THERM SY-TMERM. SXY-THERM. 


UO 

co 


20 STRESS/STRAIN constraints 

0. 59 C 7EtOO 0.1409E«-01 0.8803ETOO 

0.107 lE-fO t 0.92J0EE00 0.730JE+00 


5 TC 24 

0.11 20EE 0 1 0. 1080E + 01 

0 • 1 2 7 OEFO l 0.1 296EE01 


MOST CRITICAL CONSTRAINTS O. 5907344EE00 
0.9198E+00 0.1080E»01 0.9198E400 0.1391E+01 0.6093E+00 
0 . 7 0 40 E *0 0 0.8865E + 00 0.1113E+01 0.1169E+01 0.8306E«-00 


4 

CONSTRAINTS 

OUT 

OF 

4 

CUTOFF 

POINTS 

0. 883 71 7E100 






— 

_L .. 

2 



3 

4 




... 




0 

CONSTRAINTS 

OUT 

OF 

0 

CUTOFF 

POINT = 

0. 88371 7EE00 






10 

_. CONSTRAINTS 

OUT 

OF 

20 

CUTCFF 

PCINT= 

0.959073E+00 







s 

7 



10 

1 2 

14 

16 

17 

20 

21 

24 

10 

CONSTRAINTS 

OUT 

OF 

10 

RETAINED 

DUE TO 

VAPI ABLE LINKING 






— 

__ _5 _ . 

7 



to 

12 

1 4 

16 

17 

20 

21 

24 

0 

CONSTRAINTS 

OUT 

OF 

0 

CUTCFF 

POINTS 

0.959073E+00 






0 

. . CONSTRAINTS 

OUT 

OF 

0 

CUTOFF 

POINTS 

0.959073E+00 







AVAILABLE INTEGER ARRAY* 


314 


2 SO C 


OVERLAY ANALYS REQUIREMENTS 


POSTURfc TABLE 


retained 

TOTAL 

TYPE 

MEMBER 

NODE DIRECTION 

L«C * 

MODE 

CONSTRAINT VALUES 

D ISPLACEMENT 

CONSTRAINTS 

MOST 

CRITICAL = 

-0. 1626J1E+00 




I 

1 



2 4 

1 

-1 

0 • 5 1 2276 E— 01 

2 

2 



2 4 

1 

1 

-0.512276E-01 

3 

3 



2 £ 

1 

-l 

0* 162831E400 

4 

4 



2 £ 

1 

1 

— 0* 162831 E 4-00 

STRESS/STRAIN 

constraints 

MOST 

CPITICAL = 

0.590734E+ OC 




. . .. _ B 

5 

1 

1 


1 

-i 

0.590734E49Q 

6 

7 

1 

2 


1 

- 1 

0.680251E400 

7 

10 

1 

3 


1 

1 

0.919752E400 

8 

12 

1 

4 


1 

i 

0.919752E+00 


14 

i 

s 


i 

i 

0.609273E+00 

10 

16 

l 

6 


1 

i 

0 ,92902 1 E 400 

1 1 

17 

l 

7 


i 

-i 

0.730270E+00 

12 

20 

1 

a 


i 

i 

0.704050E400 

13 

21 

i 

9 


t 

-i 

o,ea6S|2E*oo 

14 

24 

1 

1 0 


1 

1 

0 . 830650E 4 00 


MODE STANDS FOP THE FCLLC»ING 

f 1 EG A-LVE=UCV.EF BOUND FC S I T I V6 = UPP EP BCUNC 

FOR STRESS CONSTRAINT, (CCDEF1) 

1 = VCN MISES EQUIVALENT STRESS 

2 = LONGITUDINAL STRAIN 

... 3 .= TRANSVERSE strain 

A = SHEAR STRAIN 

5 = FIRST EOLATION OF STRESS INTERACTION 

6 = SECOND EQUATION OF STRESS INTERACTION 
_7 = TSAI-A2ZI CRITERION 

FOR FRECUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 


AVAILABLE R EAL ARR AY * _7500 OVERLAY ANALYS REQUIREMENTS 304 


SELECTION OF 10 1ST ORDER APPROXIMATED STRESS CONSTRAINTS 


NEW NUMBER OF LINEARIZED CONSTRAINTS - NTCE 


14 


ANALYSIS TIME DATA 

assemble mass/stiffness matrix 

E.MdLE LOAD . VECTORS 

DECOMPOSE STIFFNESS MATRIX 
SOLUTION OF DISPLACEMENTS 
FREQUENCY ANALYSIS 
..F.LliJTEiL ANALYSIS 
CONSTRAINT EVALUATION 
POSTURE TABLE SET 
SELECTIVE GRADIENT EVALUATICN 


GRAND TOTAL CPU TIME 


0.279083E-01 
0.320435E-03 
0. 732422E-03 
0.732422E-03 
0.0 

0. 1C6B12E-03 
0 . 3353880-01 
0,3041 0 BE— 0 I 
0, 16S8J3E+00 

0. 2d625SE+00 


VO 

VO 


scaling factor 


0,1 162B3E+01 


END OVERLAY ANALYS 
SCALED WEIGHT O.S75S54E404 


NEW L I S T OF LINEARIZED CONSTRAINTS 


12 3 4 

5 7 10 12 14 

16 17 

20 

21 24 







SIDE 

CONSTRAINTS 




RELATIVE MOVE LIMIT 

VASLAdLE 

NUMBER 

0*1000 E— 0 2 

LOWER ACTUAL 

BOUND SIZE 

UPPER 

BOUND 


VARIABLE 

NUMBER 

LOWER 

BOUND 

ACTUAL 

SIZE 

UPPER 

BOUNO 


1 

J. 1 JOOt ♦ JO 

C*2000fct02 

0- 2000E+QS 

2 

0 

• 

i oooe*oo 

0 • 200 OE 40 2 

0 *2 00 OE + O 5 

3 

0. 10 0 OE + 00 

0*200 0 1 ♦- C 2 

0. 2000E+05 

4 

0 

• 

1 OOOEOO 

0*2000E-f02 

0*2 00 OE + O 5 

5 

J. 10 0 Qt£ *-0 0 

0 • 2 00 OE*- 02 

0. 2 OOOE ♦O 5 

6 

0 

• 

1 CO OE + OO 

0 • 2000E *02 

0*2 OOOE+OS 

7 

0. 1 0 0 OE *-00 

0 , 2 0 0 OE f 02 

J *2000E+05 

a 

0 

• 

1 00 0 E+00 

0* 2000E+02 

0-2Q00E + Q 5 

9 

0* lQOOEfOO 

0 • 2 00 OE 4-02 

0*2 030E+0 5 

to 

0 

• 

l OOCE+OO 

0* 2000E >02 

0*200 OE 4*05 


MOST VIOLATED SIDE CONSTRAINT 


UESIGN VARIABLE 


10 


CONSTRAINT VALUE 


0.9950E+00 



NEW LIS 


■JE LINEACI2ED CONSTRAINTS 


1 2 3 


7 i 0 


1 2 


1 4 


16 17 20 21 24 


Slot constraints 


RELATIVE MOVE LIMIT 


0. 1 COOE-02 


VARI ABLE 

NUMBER 


LOWER ACTUAL UPPER 

BOUND SIZE BOUND 


VARIABLE 

LOWER 

ACTUAL 

UPPER 

NUMBER 

BOUND 

SIZE 

BOUND 


I 

3 

S 

7 

..9 

J. I 000fc4 JO 
0. I000E4J0 
0. 1000E400 
0. 1 D00E400 

Q*. LOOSE t-M 

C.2000E402 

0.2000E4C2 

0.20006402 

0.2000E402 

0.2000E402 

0 • 20 JOEa-OS 
0. 2000E405 
0. 2000E405 
0 .2000E405 
0.2000E405 

2 0.1000E400 

4 0.1000E430 

6 0*1 CO 0E4OO 

a o.ioooe+oo 

10 O.IOOOE40Q 

0.2000E402 

0.2000E402 

0.2000E402 

0.2000E402 

0.2000E402 

0 *2000640 5 
0 *2 000640 5 
0.2000E405 
0.2000E405 
0.2000E40S 

MOST VIOLATED SIDE 

CONSTRAINT 

OESIGN VARIABLE 10 

CONSTRAINT VALUE 0.9950E400 



ENTER OVERLAY PREDL2 


AVAILABLE REAL ARRAY = 7500 OVERLAY PPEUU2 REQUIREMENTS 3B9 

AVAILABLE INTEGER AR 3 AY s 250C OVERLAY PRE0U2 REQUIREMENTS 327 


( 


loo 



MAXIMIZATION OF LUAL FUNCTION 


8 Y NEWTON 


METHOD 


iG PCIN T 

CONSTRAINT 1 

..XIEHA.TI.XN l. - . 

NORM OF PROJECTED GRADIENT 

TisT *OF“ ACT I V e constraints - 

ITERATION 2 

NCRM OF PROJECT EU GRAOIENT 
LIST OF ACT IVE CONSTRAINTS 
ITERATION 3 


DUAL VARIABLE = 


0.318613E+03 


NORM OF PRCJECTEO GRADIENT 
LIST OF ACTIVE CONSTRAINTS 


ITERATION 4 

NORM— QE. -BSOJELT ED GRADIENT 
LIST OF ACTIVE CONSTRAINTS 
ITERATION 5 

NCRM OF PROJECTED GRADIENT 


LIST OF ACTIVE CONSTRAINTS 
ITERATION b 


NORM OF PROJECTED GRADIENT 


ITERATION 


L 


NORM OF PROJECTED GRADIENT 
LIST OF ACTIVE CONSTRAINTS 
ITERATION 8 

LIST OF ACTIVE CONSTRAINTS 
ITERATION 9 

NORM OF PROJECTED GRADIENT 


LIST OF ACTIVE CONSTRAINTS 
ITERATION 10 


NORM OF PROJECTED GRADIENT 
.SX-DF.. ACJIVL£_CCNSTR&IMTSl 
ITERATION tl 


NORM OF PRCJECTEO GRADIENT 
LIST OF ACTIVE CONSTRAINTS 


0 • 81 6 56 4 E— 05 

1 7 

0.1301 77E+02 
1 7 

0.31 0776E+01 
l 7 

0.404 0 Q2EF01 
1 7 

0.1892Q3E+01 
1 7 

0.733 5 7 6E+ 0 0 
l 7 

0. 1 90092E+00 

1 7 

_ 0 . L9 4 8 9GE~ 0 1 
1 7 

0.285972E-03 

1 7 

0. I92S97E-04 

I 7 10 

0.1 19668EF01 


1 


7 


10 


101 



0.733<i 78L+00 


NOPM OF PRO JETC T E 0 G^ AD I CN'*’ 

LI ST -JF —ACTIVE -CONSTRAINTS l 7 

ITERATION 7 


'1 


NC»M OF PROJECTED GRADIENT 
LIST OF ACTIVE CONSTRAINTS 

iteration 8 


0. I 90O92E.TOO 

1 7 


NORM . OF_ PROJECT ED GRADIENT 
LIST OF ACTIVE CONSTRAINTS 
ITERATION 9 


0. 19 4 8 9 6E— 0 l 

1 7 


NORM OF PROJECTED GRADIENT 0.285972E-03 

LIST OF ACTIVE CCNSTRaYnTS " ~ I 7 

ITERATION 10 


NORM OF PROJECTED GRADIENT 



NORM OR PROJECTED GRAOICNT 
LIST OF ACTtVE CONSTRAINTS 
ITERATI ON L2 


N ORM OF P ROJECTED GRADIENT 

LIST OF ACTIVE CONSTRAINTS 
'TFatnnn L3 

NORM OF PROJECTED GRADIENT 
_ L 1ST OF AC fl V E ‘ CO N S T R A I N T S 
ITERATION 14 

NORM OF PROJECTEO GRADIENT 

LIST OF ACTIVE CONSTRAINTS 

ITERATION 15 

NORM OF PRO JECTECf ^RAIMENT 
LIST OF ACTIVE CONSTRAINTS 
iflRAT'i ON " T6 ' 


NORM -OF -PROJECTED GRADILNT 
LIST OF ACTIVE CONSTRAINTS 
.ITERATION . 17 

NORM OF PROJECTED GRADIENT 
LIST OF ACM V E’CC NS TR A I N T S 
ITERATION 18 

NCRM OF PROJECTED GRADIENT 
_L_I S T_0 F A C TIV E CCNST RAI MTS 


0. 192S97E-04 
_ l. 7 10 


0.1 19668E+0I 
1 7 10 

0.508975E+00 
1 7 10 


0 • I 5949 7£ + 00 
1 7 10 


0.244894 E— 0 1 
1 7 10 


0.746817E-03 
I 7 10 


0.202480E-C5 
I 7 10 4 


0.388537E+00 
I 7 4 


0. I990I9ET00 


O 

- 


7 


4 



ITERATI3N 19 


NORM OF PRO JECT £0 GRADIENT 0.755865E-01 
LIST OF ACTIVE CONSTRAINTS 1 7 4 
ITERATION 20 


tiCRM QF PRCJECT.£D GRADIENT 

LIST OF ACTIVE CONSTRAINTS 

UfiBATIatf 21_ 

NORM OF PROJECTED GRADIENT 
lTst of' ACTIVE CONSTRAINTS 
ITERATION NUMBER 
PRIMAL - VArTaBLES' EVALUATIONS 
NORM OF DUAL FUNCTION GRADIENT 
~~~OUAL "OBJECTIVE FUNCTION ~ 
FINAL WEIGHT 

LIST OF ACT XVE~CONSfRAl NTS 
LINEARIZED CONSTRAINTS 


CONSTRAINT CURRENT VALUE 

NC . 

0.999993L+00 
0.999993E+00 
0.999>99IE + 00 
0 • 99999 IE + 00 
0.36 fco 7 8E *00 
0.44 1806E+00 
0.999997E+00 
0. 974624E-01 
0.292556L+00 
0.2Bl930t *00 
0.40 92 73E TOO 
0. 512176E+C0 
0.352379E-01 
0.504031E+00 



0.121933E-OI 
1 7 4 


0.3301 7 1 E— 03 
1 7 4 

22 
47 

0.250481 £— 05 
473.837 1 0 9E + 0 l 
473.63671 9E+ 0 1 
17 4 


LIMITING VALUE 
LOWER UPPER 


OUAL VARIABLE 


-0. 999992E+30 
-0.9999926+30 
-0.999993E+30 
-0. 999993E+30 
— 0 . 99999 6E + 30 
-0.999998E+30 
— 0.99999 9E+ 30 
-0.999998E+30 
-0. 99999 7E+30 
-0.999999E+30 
— 0. 999997E+30 
-0.99999 7E + 30 
-0. 999999E+30 
-0.9S9999E+30 


0.999992E+00 
-0.999992E+00 
0.999993E+00 
— 0.999993E+00 
0.999996E+00 
0.99999 8 E+ 00 
0.9 9999 9E + 00 
0.99 999 6E + 0 0 
0.999997E+00 
0.9999 9 9E+00 
0.999997E+00 
0 .99999 7E + 0 0 
0.999999E+00 
0.99 99 9 9 £ + 0 0 


0.434644E+04 

0.0 

0.0 

0 .6 6S 807E+ 0 3 
0.0 
0.0 

0.I05776E+0A 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


H 

o 

CO 


END OVERLAY PRE0L2 




0 .c 


RESFCNSE FACTOR HEUUCEU TO 

TRUNCATION FACTORS MOoIFIEJ AS FOLLOWS 

JUS E-LA CEMENT ..CONST SAINTS 0. 1200 0 OE *0 0 

STRESS/STRAIN CONSTRAINT 0.120000E+00 

SELECTION FACTOR (F.S.D.l 0 . IOOOOOEj-OO 


UPJJATEi! SCALING FACTORS 

J.12J1ET01 0.4107E+00 


0.6977E-01 0* 3 2 40E + CO 0.1.1B9ET01 


UPOATEO HEIGHT COEFFICIENTS 

JX-a-OEE 1 £ ±03 0.2aiZ£±0i._ G.5O23E*02 0.2333E±Q3 

0.0 0.0 


0.QS59Et03 







STAGE M3. 


APPROXIMATL problem generator 


CURRENT MEMBER SIZE 

MEM UER~TYPE' 'NUMBER 1 

0 .246 IE +0 2 0 . 32 1 4E+ 0 1 0.139SEH01 0.6480EH01 0.237/E«-02 0.3031EFCI 0.1172E402 0. 1284E4-02 0.9164E»01 0.1162E+02 


CURRENT HEIGHT CATa 

member type number 


WEIGHT = 


0 .4 738J7E+04 


VARIABLE STRUCTURAL H EIGHT 0.473837E+04 
F I XEQ STRUCTURAL HEIGHT 0.0 

TOTAL STRUCTURAL WEIGHT 0.47383 7E+ 04 

NCN-STR UCTURAL WEIGHTMT 0,0 


TOTAL HEIGHT 


0.473B37E+04 


CON VERGENCE CHECK STAGE NO , = 2 O.T713EHOO 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 


0.U91E + 2? 
0.1 OOOOOEH31 


MUST BE LESS THAN O.lOOOOOE-03 
0.839292ET04 0 .473837EH04 


ENTER OVERLAY ANALYS 


AVAILABLE IN TEGER ARR AY = 250C OVERLAY ANALYS REQUIREMENTS 314 


POSTURE TABLE 


RETAINED TOTAL 


TYPE MEMBER 


NODE DIRECTION 


CONSTRAINT VALUES 


DISPLACEMENT CONSTRAINTS 
1 1 

2 . 2 

3 3 

4 4 

_SJB£LSS/S.TRAIN. CONSTRAINTS 


MOST CRITICAL * 


MOST CRITICAL 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 


-0.691996E-0 1 
2 4 

2 4 

2 S 

2 E 

0*38 70 08E+ 0 C 


; ST ANDS FQF_ the FCLLQhJNG 
NEGATI VE-LCHER BOUNO POSIT I VE*UPPER BOUNO 
FOR STRESS CONSTRAINT. (COOE+11 

1 s VCN MISES EQUIVALENT STRESS 

2 = LONGITUDINAL STRAIN 


-0.69I996E-01 

0.691 996E— Ol_ 

-6.61S72TE-01 

0.615721E-01 


0.6S7426E400 
0.6171 42E + 00 
0.387008E400 

&j,867999g400_ 

0.681680EHOO 
0.860421 £400 
0.569398E400 
.. _ 0.51|736g.t00 
0.867999E+00 
0.617 144E+00 


105 


l * » J *4 ‘ * f-uo l I 1 »L-Uhr'tl' l) l. U IN U 

FOR STRESS CONSTRAINT* (LODE+1) 

t - VCN MISES ECU I V AL fc NT STRESS 

i = VSi|j5IH?S& A fe.H?A ,N 

4 = SHEAR STRAIN 

5 s FIRST EQUATION OF STRESS INTERACTION 

6 = SECOND EQUATION OF STRESS INTERACTION 

7 = TSAI-AZZI CRITERION 

FOR FREQUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 


AVAILABLE REAL ARRAY = 7500 OVERLAY ANALYS R EQU IREK ENT = 304 


END OVERLAY ANALYS 


3 CALX EG FACTOR 0^.1 06 92 QE tO 1 


SCALED WEIGHT Q.506626E*Q4 


NEW LIST OF LINEARIZED CONSTRAINTS 

- 2 3 * 10 1617 20 21 24 


SIDE CONSTRAINTS 


RELATIVE MO VL LIMIT 

VARI ABLE 
NUMBER 

0. 

LOWER 

BOUND 

l Q00E-C2 

ACTUAL 

SIZE 

UPPER 

BOUND 

VARIABLE 

NUMBER 

LOWER 

BOUND 

ACTUAL 

SIZE 

UPPER 

BOUND 

- - l 

0. 0432c-*’ 01 

C.246 L E ♦ 02 

0 • 246 1 E + 0 5 

2 

0 .3 145E+0 l 

0 . 82 1 4E 401 

0*82 1 4E4 0 4 

3 

0. 1 00 OE«-00 

0.1 395E + 0 1 

0. 1395E+ 04 

4 

0 .8553E400 

0 • 648 OE 40 1 

0.6480E+04 

5 

0. 7‘56BE+0I 

0.237 7EFC2 

0. 2377E+05 

6 

0 • 1 00 OE + OO 

0.303 ie*ot 

0.303IE404 

. 7 

Q.l aODETOC 

Q.1 172E+02 

0.1 1 72Et 05 

d 

0 .IOOOETOO 

Q.I2U4E4Q2 

0.1284EtQ5 

9 

0. 1 QOOtTOO 

0.916 4 ET 0 1 

0.9 164E+04 

10 

O-IOOOE+OO 

0. 1162E+02 

0.1162E+03 


..MOST _Y1 PLATED SIDLC CJiS.T RA IHT . r DESIGN VARIABLE 2 CONSTRAINT VALUE Q_. 617.l£tfli>. 


ENTEP OVERLAY PRECL2 


AVAILABLE REAL ARRAY = 
-ARRAY- 


7b 0 0 OVERLAY PRE0U2 REQU I R EM ENT = 

2500 OVERLAY PREDU2 REQUIREMENT= 


275 

323 


l 


I TE&ATI ON 1 

NORM Or PRC JLC T ED GR A J I lNT 
LIST OF ACTIVE CONSTRAINTS 


MAXIMIZATION OF DUAL FUNCTION BY NEwTCN METHOD 3=s==s 

0.42tjB30E + 01 
1 4 5 


I TERATI ON 


2 


10$ 



===== MAXIMISATION Of DUAL fUNCTION BY NEwTCN METHOD 

XTERAXION X_ 

0.42cd30E+ J1 
1 4 5 


NORM OF PROJECTED GRADIENT 

LIST ..OF. ACT-IYt CONS I RAIMIS 

ITERATION 3 


0.39O776E+00 
l 4 5 


NORM OF PROJECTED .GRADIENT 
LIST OF ACTIVE CONSTRAINTS 
I.T£HALL2N_ . .5 . .. 

NORM OF PROJECTED GRADIENT 
Li st OF ACT I VE CON STR A I N t S 
ITERATION 6 


0.667348E-03 
1 4 5 

7 

26 

0.636S49E-05 
4 3 9*0 l 36 7 2E+ 9 1 
439.0 I40E2ET01 
1 4 5 


NORM OF PROJECTED GRADIENT 

ITERATION NUMBER 
.£aiMAL_.VLAa£Aai_ES.. EVALUAT IONS 
NORM OF DUAL FUNCTION GRADIENT 


FINAL WEIGHT 

.LIST OF ACTIVE. CONSTRAINTS 
LINEARIZED CONSTRAINTS 

=========== 3 = 333 ====== 


0. 181875E+00 
1 4 S 

0. I 3S036E-0I 
1 4 S 


'NORM' ~OF PROJECTED GRADIENT 
LIST OF ACTIVE CONSTRAINTS 
"ITERATION 4 


0. 4250990+01 
i 4 5 


NCRM Or PROJECTED GRADIENT 
L I ST OF ACT I ve"CON S TRA IN TS 
ITERATION 2 


CONSTRAINT 

NO 


current value 


LIMITING 

lower 


VALUE 

UPPER 


DUAL VARIABLE 


1 

2 

J 

4 

5 

6 
7 

a 


0. 9999960*00 
— 0 • 9 99996 0 1 0 J 
0. 999996L +00 
—0.999 9960+00 
0. lOOOOlE + OI 
J. 3094 1 9fc +C0 
0.396IS9E+00 
J. 5 2 40550 +00 


-O. 99999 SE + 30 
-0.999998E+30 
-0.999998E+30 
-0. 9999980+30 
-0. 1 OOOODE +31 
-0. 99999 6C +30 
-0. 999999E+30 
-0. 9 9999 dE + 30 


0.999998C+00 
-O. 9999960+00 
0.999998E+00 
— 0*99999eE+00 
0.1 000000+0 l 
0.99999ee+00 
0 .999999E+00 
0.999998E+00 


0.381 796E+04 
0.0 
0.0 

0. 4152566+03 
0 .987451E+03 
0.0 
0.0 
0.0 


9 


0.36161 4t— 0 1 


-0. 100000E+31 


0.1 0 00 0 00 + 0 1 


0.0 





L I NEA rtl ZL J CCNi A 1 NT 9 


CCNSTRA int 
NC 


CURRENT VALUE 


limiting value 
lower upper 


DUAL VARIABLE 


J 

4 

6 

6 

7 

a 

9 

1 C 


9 . 9999960+00 
— 0.999996E + JJ 
0. 999996L +00 
-0.999996E+00 
0. 1 OOOOIE + 01 
Q.3Q5419E+C0 
0.398l‘J9E + C0 
0. f< 405SE+C0 
0.26161 4E- J l 
0. 4 1 26 JOE + 00 


-0. 999998E+30 
-0.999998E+J0 
-0. 999998E +30 
-0.99999at+30 
-0.1 OOOOOE+Jl 
-0.999998E+J0 
-0. 999999E+30 
-0. 99999 at + 30 
-0. I 00 0 0 OE + 3 I 
-0. I 0000 OE +3 l 


0.99999EC+00 
-0.99999EC+00 
0. 999990E+OO 
-0.9999960+00 
0.1 OOOOOt+O 1 
0.999996E+00 
0.999999E+00 
0.99999QE+00 
0.1 0 00 0 OE +0 1 
O.tOOQOOE+Ol 


0.381 796E+04 
0.0 
0.0 

0.415256E + p3_ 
0 .987451E+03 
0.0 
0.0 
0.0 
0.0 
0.0 


END OVERLAT PREDU2 



RESPONSE FACTOR REDUCED TO 


TRUNCATJCN FACTORS MODIFIED AS FOLLOWS 


D I SBL.ACEMENT CONSJR A INTS 

STRESS/STRAIN CONSTRAINT 
SELECTION FACTOR (F.S.D.I 


0* 3583 I 7E + 00 
0.3583 l 7E ♦ 00 
0 .1 OOOOOE+OO 


BilAXEa_S.CAl_l^C.,EACrOR5-_ . 

0.1079EF01 0.4217E+00 0.5000E-02 0.7006E-01 0.ll33E*0l 


UPDATED WEIGHT COEFFICIENTS 

Q.Z76 dE*Q3 Q* 3036E ♦ 03 0.3600ET01 0.5044E+C2 0.81S9EF03 

0.0 0.0 



i 


‘ 


* 


0.5000E-C2 0.7272E400 0.6347EFC0 C.9908E-0I ' 0.5964E4-00 
0.36 0 OEF 01 Q.74 05EF03 0.6463E»03 Q.1Q09E*03 .... 0.6072EF03_ 




STAGE NO. 3 APPROXIMATE PROBLEM GENERATOR 


CURPEN T MEMBER SIZE 


MEMBER TYPE NUMBER 1 

3.2 i 5 8E +02 0.8434EF01 0.1000EF30 0.I4U1E+01 0.22<i6E*02 O.IOOCEFOO 0.1454E + 02 0.1269E+02 G.I982E+01 


’CURP.EMT * El GH T CAT A 


MEM B £ R TVPE NUMBER l WEIGHT = 0 . 4 3 4 38 1 £«■ 04 


VARIA6LE STRUCTURAL WEIGHT 0.404881E+04 
FIXED STRUCTURAL WEIGHT 0.0 

TOTAL STRUCTURAL WEIGHT 0 .4 0 4 88 1 EF 04 

NCti” STRUCTURAL WEIGHTHT 0,0 

TOTAL WEIGHT 0 .40468 1 E+ 04 

iaNViLECiLNCE CHECK STAGE Nl.= 8 0.E246c-04 0.9V22E-03 MLST BE LESS THAN O.lOOeOOe-83 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0.404501EF04 0.404S03E+04 0.404881EF04 


ENTER OVERLAV ANA L VS 


AVAI LAU LE INTEGER ARRAY= 2300 OVEPLAy ANAL V S REQUIREMENT 2 310 


POSTURE TABLE 


RETAINED TOTAL TYPE MEMBER NODE DIRECTION L.C. MODE CONSTRAINT VALUES 


DISPLACEMENT 

CONSTRAINTS 

MOST 

CRITICAL 2 

-0.371933E-04 




1 

1 



2 4 

1 

-l 

-0.37I933E-04 

2 

2 



2 4 

I 

1 

0.371 93 3 £—04 

3 

3 



2 S 

1 

-I 

-0.209808E-04 

4 

4 



2 5 

I 

1 

0.209808E-04 

S-THESS/STP.AIN 

CONSTRAINTS 

MOST 

CRITICAL 2 

— 0 .376562 E- 04 




5 

S 

1 

1 


1 

-l 

0.625415ET00 

6 

7 

1 

2 


I 

-I 

0.537578E+00 

7 

l 0 

1 

3 


1 

I 

-0.376562E-04 

a . 

.. L7 

1 

7 


1 

-1 

.Q*$L9U9£t.00_ 

9 

20 

1 

8 


t 

1 

0.S45171E400 

10 

24 


10 


I 

I 

Q.337578EFOO 


..MODE STANDS FOR t HE FOLLOWING 

NEGATI VE=LCWER BOUND FOS I TI VE=UPPEP BOUND 

FOR STRESS CONSTRAINT. (CODEFIJ 

l = VCN MISES EQUIVALENT STRESS 
_ _2 = LONGITUDINAL STRAIN 

3 2 TRANSVERSE STRAIN 

4 = SHEAR STRAIN 

5 2 FIRST EQUATION OF STRESS INTERACTION 

6 2 SECOND EQUATION OF STRESS INTERACTION 

7 2 TSAI-AZZI CRITERION 

FOR FRECUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 


AVAILABLE REAL ARRAY = 7500 OVERLAV ANALYS REQUIREMENT 2 300 


END OVERLAY ANALYS 

SSALLNG FACTOR 0. 1 00004E TO I SCALED WEIGHT 0.404896EF04 


NEW LIST OF L I NEAR I ZEO CONSTRAINTS 
i 2 3 4 10 


.1 193E+02 


SIUE CONSTRAINTS 


110 


M . IN 


4 = SFcAR STRAIN 

5 = FIRST EQLATIGN OF STRESS INTERACTION 

6 = SECOND EOLATION OF' STRESS INTERACTION 

7 = T S A I — AZZ I CRITERION 

FOR FRECUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 


AVAILABLE REAL ARRAY = 7500 OVERLAY ANALYS REQUIREMENT* 300 


END OVERLAY ANALYS 


£CA_LI.NG FACTOR 0.100004EF01 


SCALED WEIGHT 


0.404896E+04 


NEW LIST OF L I NEAR I ZEO CONSTRAINTS 
- g 3~ 4 10 " 


SIDE CONSTRAINTS 


RELATIVE MOVE LIMIT 


0.1 000E-02 


VARIABLE LOWER ’ ACTUAL UPPER 

NUMBER BOUND SIZE BCUND 


VARIABLE LOWER ACTUAL' UPPER 

NUMBER BCUND SIZE BOUND 


1 0 .80d3L+01 

3 0.1000E+00 

5 0. 100 OE -F 0 0 

— / P.5&40EFQ1 

9 0.10 OOE + OO 


0.2158E+02 
0. 1 00 OE* 00 
0.2266E+02 
0tl454EF0a 
0 . 1982EF0 1 


0 .2 158E+05 
0. 1000E+03 
0. 2266EF05 
0.1454E+05 
0. 1982E+04 


2 0 .3900E+0 1 0.8434EF01 0.S434E4-04 

4 0.100 OE + OO 0.1401EF01 0.1401EF04 

6 0.1 OOOE+OO 0.1000E+00 0.1000EF03 

d 0.5774EF01 0.1269E+02 0.1 269EW05 

10 0.5S1SE+01 " 0. 1193EF02 0.1 193E+0S 


MOST VIOLATED SIDE C C NSTR A I NT - DESIGN VARIABLE 6 CONSTRAINT VALUE 0.0 


ENTER OVERLAY PREDU2 


AVAILABLE REAL ARRAY = 7500 OVERLAY PRE0U2 REQUIREMENT* 155 

AVAILABLE IN T EG ER ARRAY* 2500 OVERLAY PREDU2 REQUIREMENT* 314 
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LTEfiATtnN 1. 

NORM OF PROJECTED GSAJltNT 

LIST of active constraints 

ITERATION NUMBER 
PRIMAL VARIABLES EVALUATIONS 
NORM OK OUAL FUNCTION GRADIENT 
DUAL OBJECTIVE FUNCTION ~ 

FINAL WEIGHT 

~ "LIST OF ACTIVE CONSTRAINTS 
LINEARIZEO CONSTRAINTS 


MAXIMIZATION OF DUAL FUNCTION UV NEKTON METHOD 


0.567979E-03 

1 3 6 

2 
2 

0.98a535E-06 
4O4.B9S5C8ET01 
40 4.B9 ‘jS5 7ETOI 
1 3 5 


CONSTRAINT 

. . MO 


CURRENT VALUE 


LIMITING VALUE 
LOWER UPPER 


DUAL VARIABLE 


1 

2 

J 

V 

5 


0. 999997E +00 
-0.999997E+00 
0.999998ET00 
-0. S9999bE 400 
J. lOOOOOE+Ol 


-0. 399997ET30 
-0.999997E+30 
-0.999998E+30 
-0. 999998ET30 
-0. I OOOOOE 4 3 1 


0 #9 999 9 ?£ 40 0 
-0.999997ET00 
0.99999EE+00 
-0.999998E+00 
0. 1 0 CO 0 QE + 01 


0. 302979E+04 

0.0 

0 • 354674E-F 03 

0.0 

0.S32447EF0 3 


END OVERLAY PRE0L2 


RESPONSE FACTOR REDUCED TO 0.0 


TRUNCAT I CN FACTORS MOOIFIED 


S TP ESS/STRAIN CONSTRAINT 
SELECTION FACTOR (F.S.D.) 


AS FOLLOWS 
.... 0..4299BOE400 
0.429980E+00 
0.1 000 0 OE TO 0 


UFJATED SCALING FACTORS 

0. 107 9E +01 0.421 7ET00 0.5000E-02 


0.7006E-01 0.1 1 33E TO l 0.5000E-02 0.7272E +00 0. 6347ET COQ.99O8E-010.5964ETO0 


UPDATED WEIGHT 


COEFFICIENTS 
_a.J0A&EAflJ ... 
0.0 


QJ&QQEtQI 


0.0 


0.5044£t02 0<81 59ET.Q3 0 » 36 3 QE T.OI 0.7405 ETQ3 0.6463ET03 0.1009ET03 0.6073ETQ3 


112 



0.0 




0.5000E-02 0.7272E<-00 0.6347E+C0 0.990BE-01 0.S964E400 

0,3&OOe»OJ 0, 74056403 0* 6463E+ Q3 0.1OQ9 E+03 0.6073E403 




9 


STAGE NO . 


APPROXIMATE 


PROBLEN GENERATOR 


CURRENT MEMBER SIZE 


MEMBER TYPE NUMBER 1 

0.2 1 58E +02 0.84J4E+01 O.IOOOEfOO 0.1401E*01 0.2266E+02 O.IOOCE+GO 0.1454E402 0.1269EF02 0. 1962E+01 0.U93E+02 


CURRENT HEIGHT CATA 


MLMUER Type NUMeEP 1 WEIGHT = 0 . 404 d96E* 0 4 


VARIABLE STRUCTURAL WE I GhT 0.404BS6ET04 
FIXED STRUCTURAL WEIGHT 0.0 

TOTAL STRUCTURAL HEIGHT 0.404896£*04 

NCN -STRUCTURAL WEIGHT.HT 0.0 


TOTAL WEIGHT 0.404896E+04 

iONV ERiLNCE _CF ECK STAGE ND.= 9 3.3497E-04 0.S246E-04 MOST BE LESS THAN O.lOOOOOE-03 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0.404903ET04 0.404881E*04 0.404896EF04 


ENTER OVERLAY ANALYS 


it* 


NQUE_ X 


NODAL DISPLACEMENTS 
Y Z NODE 



LOAD CONDITION l 


1 0.0 

3 -0. 7332 7E«-Q0 

5 0.31 441ET3C 


0.0 

— 0.2 900 OE TO 1 
-0. 1 OOOOE *-01 


0. 0 
0.0 
0.0 


2 

4 

6 


-O.JJ711ETOO -0.1IS58E+01 0.0 

C . 3 7 864 E TO 0 — Q.20000E + 01 0.0 

0.0 0.0 0.0 


:ss::::25ss=z:sss3is zssssasszsssssxassissaa 


- 4 DI 5PLACEVQN- C 3NS T P AI NT S 


TO 4 MOST CRITICAL CONSTRAINT* —0.143051 I E— 05 



NGC AL 


DISPLACEMENTS 


NODE 


X 


Y 


2 


NUDE 


L C AD CCNDITICN 1 


I 

3 

5 


0.0 

-0.75327E+00 

0.31441E+QC 

= = = = = = = =35-= = 


0.0 

- 0 . 290006+01 

-0.1 OOOOE+Ol 


0. 0 
0.0 
0.0 


2 

4 

6 


-0.33711E+00 — 0 . I 1 553E + 0 1 0.0 

C. 3 7 864E + 0 0 -0.20000E+01 0.0 

0.0 0.0 0.0 

3X = = S53S»35«S3S3SIZI»|»SSSailM«: 


4 X)I iELACEVEN’- C JNSTF AI NTS 

0. I 431E-05 -0.1431E-05 0.1132E-05 


1 

-0.U32E-05 


TO 4 


MOST CRITICAL CONSTRAINT* — Q .1430511 £-05 


MIXR . M _Li: S-CCMHlhED SX SY SXY . SX— THERM SY - THERM SXY-THERM 


1 

1 1 

-0* 93e425E*04 

-0 .9364E+04 

0.0 

0. 0 


I 

1 i 

-0. 93e425fc*04 

-0 .Q364E«-0<* 

0.0 

0.0 


_ . . . J. 

- 2_ I 

-Q.11560QE+05 

— 0.11 5 6E + 05 

0.0 

0*0 


1 

2 I 

-0.11 5600E+05 

-0 . 1 156E+05 

0.0 

0.0 


1 

3 1 

0.250000E+05 

0 .25006+05 

0.0 

0.0 


1 

2 1 

0.2500COE+05 

0.250 0E + 05 

0.0 

0.0 


1 

4 1 

_ fl*.17fl421£*0 4 

.0*178 4E + 04 

0*0 

0*4 


1 

4 1 

0. 17342IE+04 

0.1 734E+04 

0.0 

0.0 


I 

5 l 

0.3733S5E+04 

0.E734E+04 

0. 0 

0* 0 


1 

5 1 

0 .U73355E+04 

0 .3734E+04 

0.0 

0.0 



6. _ 1__._ 

Q.432727E+04 

0 . 432 7E + 04 

0.0 

O.Q 


l 

6 I 

0.432727E+04 

0.4327E+04 

0.0 

0.0 


1 

7 i 

-0. 9522C9E+04 

-0.9522E+04 

0.0 

0. 0 


1 

7 1 

-0.9522C9E+04 

-0 .9522E+04 

0.0 

0.0 



s i 

. Q*1 L37Q4E+Q5 

0. I U7E + Q5 . 

0.0 

0. 0 

1— ' 

1 

a l 

0. 1 1 37046+05 

0 . 1 137E+05 

0.0 

0.0 

M 

1 

9 1 

-Qm 1 78422E+04 

-0. 1784E+04 

0.0 

0.0 

in 

1 

9 1 

-0. 1 7 34 22E + 0 4 

-0 . 1734E+04 

0*0 

0.0 


10 1 

Q.llE600Et05 

0.11 56E+05 

0.0 

0. 0 


1 

1 0 1 

0.11 £6006+05 

0. 1 156E+05 

0. 0 

0. 0 



STRESS/STRAIN CONSTRAINTS 

k + Ql_ 0 «5376t+ 00 


TO 24 


0.1462E+Q1 0 • 2000E 1 0 1 -Q.4637E-Q6 


MOST CRITICAL CONSTRAINTS -0. 4637500E-06 


♦01 0.8269E+00 

0 . 

6191 E+ 00 0. 13316 + 01 

0.1455E+01 0.54! 

CONSTRAINTS 

OUT 

OF 

4 

CLTCFF 

PCINTs 

0.57002 OE + 00 

l 

2 

— 


3 

4 


CONSTRAINTS 

OUT 

OF 

0 

CLTCFF 

POINT* 

0.570020E+00 

CONSTRAINTS 

7 

OUT 
1 0 

OF 

_2Q 

CUTCF1 

20 

POINT* 

24 

Q.57Q020E+00 

CONSTRAINTS 

7 

JUT 

10 

OF 

4 

RETAINED 

20 

OLE T 0 
24 

VARIAbLE LINKING 

CONSTRAINTS 

OUT 

OF 

0 

CUTOFF 

POINT* 

0.5 70 0206 + 00 

CONSTRAINTS 

OUT 

□ F 

0 

CLTCFF 

PCINT= 

0. 5700206+00 


0. 1O71E+01_ 
0.9236E+00 


-Q«9236E+QQ a, 1349E+ 0 L_JU£50 7E+0 Ql_ 

0.10716+01 Q. 1462E+01 0.5376E+00 


availaole Integer ap°ay= 


2 30 0 


OVERLAY ANALYS REQUIREMENT: 


3oa 



POSTLRE TABLE 


RETAINED TOTAL TYPE MEM 3ER NODE DIRLCTICN L.C. MODE CONSTRAINT VALUES 




0 ISFLACLMtNT CONSTRAINTS 

MOST 

critical - 

—0. 143051 E— 

05 






1 

l 



2 

4 

1 

-1 

0. 143051E-05 



2 

2 



2 

4 

1 

I 

— 0. 143051 E— 05 



3 

3 



2 

c 

1 

- 1 

0. 1 13249E-05 



4 

4 



2 

c 

1 

1 

-O.II3249E-OS 



STRESS/STRAIN C3NSTRAI NT 3 

MOST 

critical = 

-0 •4taa750E 

-oc 






. 5 

7 

1 

2 



I 

-l 

0.537601 EHOO 



6 

1 0 

1 

3 



1 

1 

— 0 . 46 8 7 50 E— 06 



7 

20 

l 

a 



I 

1 

0.S45185EH00 



H 

24 

1 

1 0 



I 

1 

0.537601 E+00 

“ ■ 


MODE STANDS 

FOR THE FOLLOhING 









NEGATIVE 

=L C WLR BOUND 

PQS ITI VE=UFPER 

BOUND 






FOR STRESS CONSTRAINT. (CCDEM) 

1 - V C N MISES EQUIVALENT STRESS 

2 = LONGITUDINAL STRAIN 
J = TRANSVERSE STRAIN 

A s SHEAR STRAIN 

5 = FIRST EQUATION OF STRESS INTERACTION 

6 * SECCNO EOLATION OF STRESS INTERACTION 

7 a TSAI-AZZI CRITERION 

FOR FREQUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 


AVAILABLE REAL ARRAY = 7500 OVERLAY ANALYS REQUIREMENTS 298 


SELECTION “OF “l 1ST CROER APPROXIMATED STRESS CONSTRAINTS 


NEW NLMBER OF LINEARIZED CONSTRAINTS — NT CE ■ S 


ANA LYSIS TIME CATA 

ASSEMBLE MA S S /Sf I FFNE S S MATRIX 
ASSEMBLE LOAD VECTORS 
DECOMPOSE STIFFNESS MATRIX 

SOLUTION _.CF 01 SPJ-AO£.M£NTS 

FREQLENCY ANALYSIS 
FLUTTER ANALYSIS 
CONSTRAINT EVALUATION 

POST URE TAB LE SET . 

SELECTIVE GRADIENT EVALUATICN 


0.235825EH00 
0. A 1 1987E-03 
Q.619233E— 03 
0.6J9233E-03 
0.0 

0.10O812E-03 
0 • 33 737 2E— 01 
0.222015E-01 
0. 193954 E *00 


GRAND TOTAL CPU TIME 0 . 3 l 4 0 72 EHOO 


END OVERLAY ANALYS 

DIMIN ISHING RETURN .OF THREE CCNSECUTIVE STAGES 


DESIGN TIME STATISTICS 
TOTAL 

INITIAL PREPARATION 

UESIGN PHASE 

ANALYSIS TOTAL 
OPTIMIZER TOTAL 


4. 1873 
0.0397 
4.1477 
2. S3 IS 
0.2755 


ENO OVERLAY DESIGN 


MAIN PROGRAM TIME STATISTICS 
P RE— PROCESS CR 
DESIGN PHASE 


0.2516 

4.1915 


9TT! 





ZLTT ! 
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z 


< 

a 

%- 


o 

u 


>■ 

u 

z 

UJ 

3 

o 

UJ 

a 

u. 


o 
z z 

UJ< 

-J 

CO V) 
ocn 
a uj 
a a 

i- 

o« 

z 


> 

a 

n 


i 


*tn • 

OIL O 
Z< *K 

< UJO 

KIJ •- 

_iuj a 

a 

Q < < UJ 

JI 


h a o' 

Z 01 UJ 
»■* »-« X 
IQh 


CCTDB 



T 


l**aa* ****** ****** a***************** ******************************************************************************* ****************-- 


COMPLETE OPTIMIZATION BT DUAL 1 OPTIMIZER 
ANALYSIS PRINT OUT CONTROL = 2 


DESIGN IN LINKED SIZING VARIABLE SPACE 

NUMBER CF NODES 6 
NUMBER OF TOTAL ELEMENTS 28 
NUMBER OF LINKED VARIABLES 13 
NUMBER OF LOAD CONDITIONS 2 
NUMBER OF BOUNDARY NODES 3 
NUMBER OF ISOTROPIC MATERIALS t 
NUMBER OF ORTHOTROPIC MATERIALS 4 


TRUSS 


CST CST 

ISOTROPIC ORTHOTROPIC 


SSP 


PSP 


ELEMENTS 

LINKED VARIABLES 


16 

9 


TSP 


TBD 


TBD 


TBD 


TBD 


NODE NUMBER X 


Y 


Z 


1 

2 

3 

4 
_5 

6 


0.0 

0.0 

0.0 

400. 0000 
-400.0000 

800.0000 


960. 0000 
480. OCOO 

0.0 

480. 0000 

0.0 

0.0 


6.4680 

17.6900 

5.4920 

S.9660 

5.1430 

4.0000 


-tO- 


DI SPLACEMENT BOUNDARY CONDITIONS 


BOUNDARY CODES* 

NODE NO. X Y Z 


0.0 0.0 

o.o ....... o.n 

0.0 0.0 


4 — 1— PRESCR I BED. 0*FREE. INFIXED 


1 111 0.0 

.2 1 1 .1 0.0 

3 111 0.0 


FBESCRI BED DISPLACEMENT 
X Y Z 


ELEMENT 
. _ -MQ* ... ... 

N 1 

NODE 

.... N2 

NUMBERS 

N3 N4 

LINKED 

GROUP 

INITIAL 

element 

LQtf .BO 

SIZE 

UPP.BQ 

MATERIAL 

GROUP 

SIDE CONSTRAINT 
CODE* 







CONSTANT 

STRAIN 

TRIANGULAR - 

CRTFCTRCPIC 







1 

1 

2 

4 

1 

0.400003 

0.010000 

0.020000 

-l 

- 1 

999 

2 

1 

2 

4 

2 

0.350000 

O.CIOOOO 

0.020000 

-2 


999 

3 

1 

2 

4 

2 

0.350000 

0 .010000 

0.020000 

-3 

— | 

999 

4 

1 

2 

4 

3 

0. tooooo 

0. 010000 

0.020000 

-4 

• 1 

999 

. . . s~ 

2 

5 

4 

4 

1.250000 

0.010000 

0.020000 .. 

. . -1 

— 1 

. 99* 

6 

2 

5 

4 

5 

0.750000 

0.010000 

0.020000 

-2 

• % 

999 

7 

2 

5 

4 

5 

0.750000 

0.010000 

0.020000 

-3 

— 1 

999 

a 

2 

5 

4 

6 

0.250000 

0.010000 

0.020000 

-4 

• 1 

999 

9 . ...... 

. _ _2 

3 

5 

4 

1.250000 

0.019000 . 

0*020000 ... 

. . . -1 

-1 

999. ....... 

10 

2 

3 

5 

S 

0.750000 

0.010000 

0.020000 

-2 


999 

11 

2 

3 

5 

5 

0.750000 

0 .010000 

0.020000 

-3 

-1 

999 

12 

2 

3 

5 

6 

0.250000 

0.010000 

0.020000 

-4 

_ 1 

999 

13 

4 

5 

6 

7 

0.200000 

0.010000 

0.020000 

-I 

— i 

999 

14 

4 

5 

6 

e 

0.100000 

0 .010000 

0.020000 

-2 

* i 

999 

15 

4 

5 

6 

a 

0. 100000 

0. OlOOOO 

0.020000 

-3 


999 . 


4 

5 

6 

9 

0.100000 

0. OlOOOO 

0.020000 

-4 

-i 

999 * 

_ . . — . 


- . ^ „ 






■ 


■ “ 


0 





CCNSTANT S'SAIN t 

RI ANGULAR - 

CRTHITRCPIC 








1 1 

2 

4 

1 


0.40000 J 

0. 010000 

0.020000 

-1 

— 1 

999 

2 1 

3 

4 

2 


0*350000 

0 • Cl 0000 

0.020000 

-2 

— 1 

999 

3 1 


4 

2 


0*350000 

0 • 01 0000 

0.020000 

-3 

— 1 

999 

4 1 


4 

3 


0* 1 00000 

0* 01 0000 

0.020000 

-4 

— 1 

999 

5 2 

5 

4 

4 


1 .250000 

0.010000 

0.020000 

-l 

— 1 

990 

6 2 

5 

4 

5 


0.750000 

0.01 0000 

0 .020000 

-2 

— 1 

999 

7 2 

5 

4 

5 


0.750000 

0* 010000 

0.020000 

-3 

— 1 

999 

a 2 

5 

4 

6 


0.250000 

0*010000 

0.020000 

—4 

— l 

999 

_ 9 2. 

3 

5 

4 


1 *2 500 0 0 

0.010000 

0.020000 

-1 

— 1 

999 

10 2 

3 

5 

5 


0. 750000 

0*010000 

0*020000 

—2 

— 1 

999 

1 1 2 

3 

5 

5 


0*750000 

0 *010000 

0*020000 

—3 

— 1 

999 

12 2 

3 

5 

6 


0.2 50000 

0 • 01 0000 

0.020000 

-4 

— 1 

999 

13 4 

5 

6 . 

7 


0*200000 

0* 01 0000 

0.020000 

-1 

— 1 

999 

14 4 

5 

6 

e 


0.100000 

0 .01 oooo 

0*020000 

-2 

— 1 

999 

15 4 

5 

6 

a 


0* 100000 

0*010000 

0.020000 

•3 

— i 

999 

16 4 

5 

6 

9 


0*100000 

0. 01 0000 

0*020000 

-4 

— ‘1 

999 

SYMMETRIC SHEAR PANEL 






■ ■ 


* " 

l 4 

5 


i 


0*600000 

0 .020000 

0.020000 

1 

— 1 

0 

2 2 

4 


2 


0.600000 

0 .020000 

0.020000 

1 

— 1 

0 

3 . 3 

s 


3 


0.600000 

0.020000 

0. 020000 

1 

— 1 

0 .. ... 

4 5 

6 


3 


0.600000 

0 . 02 0 00 0 

0.020000 

1 

“1 

0 

5 1 

4 


4 


0*600000 

0*020 00 0 

0*020000 

1 

— 1 

0 

6 4 

6 


4 


0*600000 

0. 02 0 00 0 

0.020000 

1 

" 1 

0 

~ THERMAL SYMMETRIC 

PANEL 






■ ~ ■ • * 

_ - 

~ 

1 4 

5 


1 


0*350000 

0*020000 

0*020000 

1 

— 1 

0 

2 2 

4 


2 


0*350000 

0*020000 

0*020000 

1 

— 1 

0 

... 3 3 . 

3 


3 


0.350000 

0 .020000 

0.020000 

t 

— 1 

0 

4 5 

6 


3 


0*350000 

0.020000 

0*020000 

1 

— 1 

0 

5 l 

4 


4 


0*350000 

0.020000 

0.020000 

1 

—l 

0 

6 4 

6 


4 


0*350000 

0 .020 00 0 

0 .020000 

1 

“1 

0 

* -2 =F IXED AT 

INITIAL VALUE -1=L0WER 

BOUNOS 

ONLY “ 

— 

.... 

— — . 

. . . 

■ m " 

0 =NCN ! 

NEGATIVITY ONLY 

1 SUPPER 

BOUNDS 

CNLY 






2=BCTH 

UPPER AND LOWER 

bOUNDS 








. . SECOND 

JJUMHER_a DISCRETE VALUE GROUP 

NUMBER 





- _ . . 

..... - 


MATERIAL CONSTANTS - ISOTROPIC MATERIALS 

GR3UP 

.... n a*... 

YOUNG* S 
MODULUS 

PO IS SON * S SPECIFIC 

RATIO WEIGHT 

THERMAL COMPRESSIVE TENSILE 

EXPANSION. A. STRESS A. STRESS 


— — to-— 

i 

16400000*0 

0.3000 0.160000 

0.00000560 -125000.0 125000.0 


to 


MATERIAL CONSTANTS - ORTHOTROPIC MATERIALS 


GROUP 

YOUNG* S 

YOUNG* S 

SHEAR 

POISSON'S 

SPECIFIC 

THERMAL 

THERMAL 


DIRECTION COSINES 

NO. 

MO OULUS 

MODULUS 

MODULUS 

RATIO 

•(EIGHT 

expansi ON 

EXPANSION 

OF LONGITUDINAL 

AXIS 


(ELI 

JET) 

1 glT) 

1NULT1 

(GAMMA ) 

(ALPHAL) 

( ALPHAT ) 

X 

r 

z 

* 1 

21 000000. 0 

"T 7oooo oVo “ 

650000. 0~ 

0.210000 

0* 056000 

—0.000 00 02 1 

0* 0000 1600 

1.0000 

0*0 

0.0 

-2 

21000000. 0 

1 700000. 0 

650000.0 

0. 210000 

0.056000 

-0*00000021 

0.00001600 

0.7071 

0.7071 

0.0 

-3 

21 000000. 0 

1 700000* 0 

650000.0 

0.210000 

0.056000 

-0* 00000021 

0*00001600 

-0.7071 

0.7071 

0.0 

-4 

21 OOOOOO. 0 

1 7Q0Q_QQ*Q 

6SQQQQ.Q 

Q *210000 

0*056000 

-OjQOQ 00021 


0*0 

1*0000 


GROUP 

STRAI N 

strain 

STRAIN 

STRAIN 

STRAIN 

STRESS 

STRESS 

STRESS 

STRESS 

STRESS 

NO * 

TEN. LI MIT 

COMP. LIMIT 

... TEN. LI MIT 

COHP .LIMIT 

SHEAR LIMIT 

TEN. LONG. 

COMP .LONG. 

TEN .TRANS 

COMP. TRANS 

_ smear 


1EPSTL ) 

1 EPS CL 1 

(EPSTT) 

(EPSCT1 

(GAMMALT) 

(FTL> 

(FCLl 

<FTTT 

(FCT) 

IFLT) 

- 1 

o.ooasn 

-0.008571 

0.004706 

-0.017647 

0.018462 

0.0 

0.0 

0.0 

0.0 

0.0 

-2 

0.008571 

-0.008571 

0.004706 

-0.017647 

0.018462 

0*0 

0.0 

0.0 

0.0 

0.0 

-3 

0.008571 

-0.008571 

0.004706 

-0.017647 

0.018462 

0 .0 

0.0 

0.0 

0.0 

0.0 

-4 

0.008571 

-0.008571 

0.004706 

-0. 017647 

0*018462 

0 *0 

0.0 

0*0 

0.0 

0.0 


= = = === = = = = = = = = = = = = = = = == == = = ssss::ssssss33sssss: ==* *m = sssssaasaassaassssisxaasxxsstaasaisssflsisiisuUBMSUMaiasaiusilMlH 


LOAD CONDITIONS 


LUMPED LOAD AT NODES 


NODE NUMBER X V Z 


LOAD CONDITION 
4 

s 


94210.0000 

94210.0000 


0.0 

n. n 


C. 0 

0.0 




\ 





1 M»llh 

5lhMN 

strain 

STRESS 

STRESS 

STRESS 

STRESS 

STRESS 

NO. 

TEN. L! M! ’ 

C 3MP . L I MI T 

TEN. L I Ml T 

COMP. LIMIT 

SMEAR LIMIT 

TEN .LONG. 

COMP .LONG . 

TEN. TRANS 

COMP. TRANS 

SHEAR 


(fcPSTL > 

(EPSCU 

< EPSTT ) 

(EPSCT) 

idAMMALT* 

(FTL) 

IFCL) 

IFTT) 

< FCT ) 

<FLT> 

- 1 

0.0085/1 

-0.008571 

0.004 706 

-0.017647 

0.010462 

0.0 

0.0 

0.0 

0.0 

0.0 

-2 

0.00857 1 

-0.008571 

0 .004 706 

-0. 01 7647 

0.0 18462 

0.0 

0.0 

0.0 

0.0 

0.0 

-3 

0.00857 1 

-0. 0005 71 

0.004706 

-0.01 7647 

0.018462 

0 .0 

0.0 

0.0 

0.0 

0.0 

-4 

0.0085/1 

-0. 0085 7 l 

0.004706 

-0. 017647 

0.018462 

0.0 

0.0 

0.0 

0.0 

o.o 


it 

M 

M 

H 

II 

II 

II 

II 

II 

II 

II 

II 




ii 

it 

li 

ll 

il 

H 

H 

il 

ll 

II 

il 

il 

ii 

3=========== 

33=333333333 33X333= 3333= 

X 33 33 33 3 3 333= 

3X333333 IS 


LOAD CONDITIONS 


LUMPED LOAD AT NODES 


NODE NUMDER 


X 


y 


z 


LOAD CONDITION 1 

4 0.0 C.O 

5 0*0 C.O 

6 0.0 0.0 

LOAD CCNDITICN 2 

4 0.0 0.0 

5 0.0 0*0 

6 0.0 0.0 


94210.0000 

94210.0000 

94210.0000 

- 94210.0000 

- 94210.0000 

- 94240.0000 


PRESSURE LOAD 


NO PRESSURE LOAO SPECIFIED 


GRAVITY LOAO 


.LOAD CONDITION NQ. MAGNITUDE! G) DIRECTION COSINES _ 

0 NO GRAVITY LOAD 

_ . _ „ ro 

u> 

THERMAL LOAD 


THERMAL LCAD GROUP NUMBER 


.ELEMENT 


LOAD CONDITIONS 


6 

6 

6 

6 

6 

-6 


type 

NUMtJER 

i 

2 3 

4 5 

6 

7 

6 9 10 11 12 

13 

14 

IS 

16 17 18 

19 20 

3 

1 

1 

1 










- . _3- - 

2. 

... . _l ..... 

1 










3 

3 

1 

1 










3 

4 

1 

1 










3 

_ 3 . 

5 

6- 

A . 

1 

1 










3 

7 

1 

1 










3 

0 

1 

1 










3 

9 

1 

1 










3 

10 

1 

1 










3 

1 1 

1 

1 










3 

12 

1 

i 










3 

13 

1 

i 










3 

14 

1 

1 










3 

15 

1 

1 










3 

16 

1 

l 











2 2 

2 £ 

2 2 

2 2 

2 2 

2 2 


THERMAL LOAD MAGNI TLOE 


- 200.0000 - 100.0000 




.S3 3 3X3X3X33 33333 



THERMAL LOAD MAGNI TLOE 


3 

6 

6 

6 

to 

to 

6 


1 6 

1 

2 

3 

4 

5 

to 


1 

2 

2 

2 

2 

2 

2 


1 

2 
2 
2 
2 

2 




200*0000 -100.000Q 


NON— STRUCT URAL MASS 


NODE NUMBER WEIGHT OF LLMPED MASS 


4 2 I 000 *0 00 0 

5 l 5000*0000 

6 12000*0000 


CONSTRAINT DATA 

SIDE CONSTRAINTS 

S I L> £ CONSTRAINT SPECIFICATIONS ARE GIVEN IN THE ELEMENT DATA 
DISPLACEMENT CONSTRAINTS 
NO DISPLACEMENT CONSTRAINTS* 


SLOPE/RELATIVE DISPLACEMENT CONSTRAINTS 


NC SLOPE CONSTRAINTS . 


STRESS/STRAIN CONSTRAINTS 


INITIAL TRUNCATION FACTOR 0*1000 

MAXIMUM TRUNCATION FACTOR _Q*9O0O 

BASIS CUTOFF FACTOR 1*0000 

MULTIPLIER FOR TfiF UPDATING 1*5000 

MINIMUM NORMALIZATION FACTOR 

STRESS 1000.0000 

STRAIN 0*1 QOOQOE-02 

NO EULER BUCKLING CONSTRAINTS IMPOSED 

ELEMENT TYPE 1 

NO STRESS/STRAIN CONSTRAINTS SPECIFIED 
ELEMENT TYPE 2 

NO STRESS/STRAIN CONSTRAINTS SPECIFIED 
ELEMENT TYPE 3 


ALL ELEMENTS ARE SUBJECT TO STRAIN ENVELOPE CONSTRAINTS 
ELEMENT TYPE 4 

ALL ELEMENTS ARE CONSTRAINED BY UPPER BOUNDS ONLY 
element type 5 

NO STRESS/STRAIN CONSTRAINTS SPECIFIED 
LLE ML NT TYPE 6 

NO STRESS/STRAIN CONSTRAINTS SPECIFIED 



ELEMENT TYPE 2 

NO STRESS/STRAIN 

ELEMENT TYPE 3 

ALL ELEMENTS ARE 

ELEMENT TYPE 4 

ALL ELEMENTS ARE 

ELEMENT TYPE 5 

NO STRESS/STRAIN 

ELEMENT TYPE 6 

NO STRESS/STRAIN 


CONSTRAINTS SPECIFIED 

SUBJECT TO STRAIN ENVELOPE CONSTRAINTS 


CONSTRAINED 

CONSTRAINTS 

CONSTRAINTS 


BY UPPER BOUNDS ONLY 

SPECIFIED 

SPECIFIED 


FREQUENCY CONSTRAINTS 


B-SPACE FIRST ORDER EXPANSION 

INITIAL TRUNCATION FACTOR O.IOOO 

MAXIMUM TRUNCATION FACTOR 0.8000 

BASIS CUTOFF FACTOR I. 0000 

MULTIPLIER FOR TRF UPDATING 1.2000 

MINIMUM NORMALIZATION FACTOR 1.0000 

CONVERGENCE CRITERIA 0.100000E-03 

ACCELELATION OF GRAVITY 0.386070E+03 

.FREQUENCY NO. CODE* LOWER BOUND UPPER BOUND 

1 -l 40.00000 0.0 


*- I=LOWER 
l=UPPER 


BOUND ONLY, 0=NO CONSTRAINTS 

BOUND ONLY, 2»L0WER AND UPPER BOUNDS 


OPTIMIZER CONTROL PARAMETERS 


OUAL 1 

PRINT OUT CONTROL 
MAX. NO. CF STAGES 
MAX. NO. OF RESTARTS 

MAX. NO. OF ONE DIM. MIN. / RESTART 
COMBINATORIAL METHOD (DISCRETE CASEI 
DIMINISHING RETURN CRITERION AMONG STAGES 
DIMINISHING RETURN CRITERION AMONG O.D.M. 
MAX. STEP SIZE ALLOWED IN A SINGLE STAGE 
STEP SIZE MODIFICATION FACTOR 
STEP SIZE MINIMUM ALLOWABLE 
STEP SIZE IN DUAL SPACE 


2 

10 

10 

100 

J» 

0.5000E-02 
0. I000E-03 
0.5000E401 
0. 7SQ0E+QQ 
0.1200E+01 
0.0 



CONSTRAINT IDENTIFICATION COOES 

CONSTRAINT TYPE I 13 CONSTRAINTS IN THIS TYPE 


sjayjQQOi 

99930002 

99930004 

99930005 

99930006 

999300 

08 

99930013 

99930Q14 

99930016 

-40002 

—40003 

-40005 








CONSTRAINT TYFE 

2 

0 CONSTRAINTS IN 

THIS TYPE 






CON 5 TO A I NT f YPfcf 

3 

o constraints in 

THIS TYPE 






constraint type 

4 

172 CCNSTRAINTS IN 

THIS TYPE 






-I 10 JO 00 l 

l 103000 l 

-1203000 l 

1203000 l 

1 303000 l 

-11 0300 

02 

11030002 

-1203 0002 

12030002 

-1 1030 003 

1 1 030003 

-12030003 

12030003 

13030003 

-11 0300 

◦ 4 

11030004 

- I 2030004 

12030004 

-l 1030005 

11030005 

-12030005 

12C30005 

1 3030005 

-11 0300 

06 

1 1030006 

-12030006 

12030006 

-11 030007 

1 1030007 

-12030007 

12030007 

13030007 

-11 0200 

Od 

1 1030008 

— 12030008 

12030008 

-1 103QQ09 

1 1 030009 

-12030009 

12030009 

1 30 3 00 0 9 

-110300 

10 

1 1030010 

-12030010 

12030010 

-I 103001 l 

1 103001 1 

-1203001 l 

1 2030011 

1 303001 1 

-11 0300 

12 

1 1030012 

-12030012 

12030012 

-1 103001 3 

1 1030013 

-12030013 

12030013 

13030013 

-110200 

14 

11030014 

— 1 203001 A 

12030014 

-1 1030015 

11030015 

-12030015 

12030015 

1303001 5 

-11 0300 

16 

1 1030016 

-12030016 

12030016 

1 0040001 

10040002 

10040003 

10040004 

I 0040005 

100400 

06 

-21030001 

21030001 

-22030001 

2303000 1 

-21030002 

21030002 

-22030002 

22030002 

230300 

02 

-21030003 

21030003 

-22030003 

230 30003 

-21030004 

21030004 

-22030004 

22030004 

230300 

04 

—2 1 030005 

21 030005 

-22030005 

23030005 

-21030006 

2 i 030006 

-22030006 

22030006 

230300 

Qo 

-21030007 

21 030007 

-22030007 

23030007 

-21030005 

21030008 

-22030008 

2203000 a 

230200 

08 

-21030009 

21030009 

-22030009 

230 JO 009 

-21030010 

2103001 0 

-22030010 

2203001 0 

230300 

10 

-21030011 

21030011 

-22030011 

2303001 1 

-21030012 

21 030012 

-22030012 

22030012 

230300 

12 

-21030013 

21030013 

-22030013 

230 J0013 

-21 030014 

21030014 

-22030014 

22030014 

230300 

14 

-21030015 

21030015 

-22030015 

23030015 

-21030016 

21030016 

-22030016 

2203001 6 

230300 

16 

20040001 

2004 0002 

20040003 


2004000S 20040006 


-40001 


13030002 

13030004 

13030006 

I303000S 

13030010 

13030012 

13030014 

13030016 

22030001 

22030003 

22030005 

22030007 

22030009 

22030011 

22030013 

22030015 

20040004 



constraint type 

4 

172 CONSTRAINTS in 

this type 






-I 10 30 00 1 

1 103000 1 

-1203000 1 

l 203000 l 

13030001 

-11 030002 

1 1 030002 

-12030002 

12030002 

13030002 

- 1 1030 003 

1 1 030003 

-l 2030003 

1 2030003 

13030003 

-11 030004 

1 1 030004 

- 12030004 

12030004 

13030004 

- 1 1030005 

1 1030005 

- 1 2030005 

12030005 

l 3030005 

-ll 030006 

1 1 030006 

-12030006 

12030006 

13030006 

- 1 1 030 00 7 

1 1030007 

-12030007 

12030007 

13030007 

- 1 1 030006 

1 1030008 

-12030008 

12030008 

13030008 

-1 1030 00 9 

11030000 

-12030009 

12030009 

13030009 

-11030010 

11030010 

-12030010 

12030010 

13030010 

-1 103001 1 

I 103001 l 

-1203001 l 

1 203001 1 

1303001 1 

-11030012 

11030012 

-12030012 

12030012 

13030012 

-1 103001 3 

1 103001 3 

-12030013 

12030013 

1303001 3 

-11030014 

11030014 

-12030014 

12030014 

13030014 

- 1 10 3001 5 

1 1030015 

-12030015 

12030015 

1303001 5 

-11030016 

1 1030016 

-12030016 

12030016 

13030016 

1 0Q4QQ0 I 

1QQ40Q02 

10040003 

10040004 

10040005 

10040006 

-21030001 

21030001 

-22030001 

22030001 _ 

2303000 1 

-21030002 

21 030002 

-22030002 

22030002 

23030002 

-21030003 

21030003 

-22030003 

22030003 

23030 003 

-21030004 

21030004 

-22030004 

22030004 

23 0300 04 

-2 1030005 

21 030005 

-22030005 

22030005 

23030005 

-210 30006 

21 030006 

-22030006 

22030006 

23 030 0 Oo 

-21030007 

21 030007 

-22030007 

22030007 

. . 2JQ3Q007 

-21030000 

21 030006 

-2203000B 

22030000 

230300 OB 

-21030009 

21Q3Q0Q9 

-22030009 

22030009 

2 30 30 00 9 

-21030010 

2103001 0 

-22030010 

2203001 0 

23030010 

-21030011 

21030011 

-22030011 

22030011 

2303001 I 

-21030012 

21 030012 

-22030012 

22030012 

23030012 

-21030013 

21030013 

-22030013 

22030013 

23030013 

-21030014 

21030014 

-22030014 

22030014 

23030014 

-21030015 

21030015 

-22030015 

22030015 

2J0JQQ1S 

— 210 300 16 

21030016 

-22030016 

22030016 

23030016 

20040001 

20040002 

-2Q040QQ3 

2Q04QQB4 

20040005 

20040006 









CONSTRAINT TYPE 
-1 

5 

1 CONSTRAINTS IN 

this type 






CONSTRAINT TYPE 

6 

0 CONSTRAINTS IN 

THIS TYPE 






NUMBER 3F BEHAVIOUR CONSTRAINTS 173 








NUMBER OF DESIGN 

variables 

13 








NUMBER OF DISCRETE VARIABLES 

9 









END OVERLAY PREPOO 
ENTER OVERLAY DESIGN 


STAGE NO. 1 APPROXIMATE PROBLEM GENERATOR 


curre - !’’ member size 


MEMBER TYPE NUMBER 3 

0.4000E+00 0.3500E+00 

0 • 350 0 E tOO 

0. 10 OOE + OO 

0. 1 250E ♦ 0 l 

0. 7 5 0 CE +• C 0 

0* 7500E+00 0 • 25 OOE*- 00 

0.1250E401 

0.7500E+00 

0.7500E+00 0.2500E+0 0 
MEMBER TYPE NUMBER 4 

O.bOOOEtOO 0*60 OOE + 03 

■j • 

o* 2 oooe*oo 

0 • 6000E 4-0 0 

0* 1 0 OOE ♦00 
0 • 6 0 OOE 4* 00 

0. 1000E+00 
0*600 OE *00 

0. 10 00£ + 00 
0 • 60 G CE + GO 





126 



I 


STAGE SO. 1 APPROXIMATE PROBLEM GENERATOR 


CUR rtI'J T MtMt3ER SIZE 


MEMBER TYPE NUMBER 3 

0.4000E+00 O . 35006*00 
0. 7500E +00 0.25006*00 

member type number a 

0.60006*00 0.6000E*00 

MEMBER type NUMBER 6 

0 . 600 OE *-00 0.60006*00 


CUR REN T WEIGHT CA’A 


MEMBER TYPE NUMBER 3 

MEMBER TYPE NUMBER A WEIGHT = 0.1761696*04 


VARIABLE STRUCTURAL WEIGHT 0 .43 I 736E*05 
FIXED STRUCTURAL WEIGHT 0.0 

TOTAL STRUCTURAL WEIGHT 0.4317366*05 

NON— STRUC T UR AL WEIGFThT 0.480000E+05 


. . TOTAL WEIGHT 0.91I736E+05 . . 

CONVERGENCE CHECK STAGE NC.= 1 0.23166*26 0.1000E*01 MUST BE LESS THAN 0.500000E-02 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0.200000E*3l 0.100000E*31 0.431736E*05 


0.36006*00 

0.20006*00 

0. 10 OOE+OO 
0« 1 0 OOE + OO 

0. 12506*01 
0.1 OOOE+OO 

0 . 75 0 CE* 00 0.75006*00 0.25006*00 0.1250E*01 

0. 10 006*00 

0.7500E*00 

0.60006*00 

0*600 OE 4-00 

0.60006*00 

0 • 60 0 CE * OO 


0.60006*00 

0* 6000E+C0 

0.60006*00 

0. 60 0 OE* 00 



WEIGHT = 0.4141196*05 


ENTER OVERLAY ANALYS 


NODAL DISPLACEMENTS 


NODE 


X 


Y 


2 


NODE 


X 


z 


LOAD CONDITIO! 1 



0.0 

0.0 

— 0.234 34 E *00 


0.0 

0.0 

0.700 7 IE— 0 1 


0.0 

0.0 

0 . 916 3 3E *01 


2 

4 

6 


0.0 

-0.48375E-01 
-0. 110 I 2E*01 


0.0 

0. 57807E-01 
0.3685 0E*00 


0.0 

0.23022E*01 
0 • 69293 E *02 


« run r^KmTTTi’M o 


12 ? 



NODAL DISPLACEMENTS 


NODE . X 


Y 


Z 


NODE 


X 



LOAD CONDITION 1 

~I _ 0.0 0.0 0.0 2 0.0 0.0 0.0 

3 0.0 0.0 0.0 4 — 0.48 375 E— 01 0.57807E-0I 0.23022E401 

S -0.23434E400 0.70071E-01 0.91633E401 6 -0. I1012E+01 0.36850E+00 0.69293E+02 

Uc Ad'cONdTt ION 2 ~ 

1 0.0 0.0 0.0 2 0.0 0.0 0.0 

3 Q.Q 0.0 0.0 4 0.48375E-0I -0.57S07E-01 -Q.23022E401 

5 0. 23434E + 00 -0.70071E-01 -0.91633E+01 6 0.U012E+0I -0.36850E400 -0.69293E402 

SS233SSSS3S33S33SS3:33SSS SSSSSS353SSS SS333 31SS3 3SSS2S8SS 5SSSS S3S3X583 3S M 331 SSSX 9MM 3 33 S U SS«3 8138 SINSa SMUtSM 



NODAL DISPLACEMENTS 


NOOH 

X 

Y 

Z 

NODE 

X 

Y 

Z 

LCAL) CCNUITI 

1 0.0 

1 

0 .0 

0.0 

2 

0.0 

0.0 

0.0 

3 0.0 


0.0 

0.0 

4 

-0.3481 JE-0I 

-0.96274E-01 

-0.54249E-02 

5 -0.78961E-01 

0 . 1 094 7E + 00 

— 0.36802 E— 02 

6 

-0. I3376E400 

0.819S2E-01 

— 0.23 098 E— 02 


L CAD CCNDI'IJN 2 


l 



NODAL OI SPLACEMENTS 


NODE . X 


Y 


2 


NODE 


X 


Y 


Z 


L CAL) CONDITION 1 


1 

3 

5 


0.0 

0.0 

-0.78961E-01 


0.0 

0.0 

0.10947E+00 


0.0 

0.0 

-0.36802E-02 


2 

4 

6 


0.0 

-0.3481 JE-01 
-0. 13376E+00 


0.0 

— 0.96 274 E— 01 
0.S1952E-01 


0.0 

-0.54249E-02 

-0.23098E-02 


LOAD CCN0ITI3N 2 


1 

3 

5 


0.0 

0.0 

-0.78961E-0 


0.0 

0.0 

0.10947E*00 


0.0 

0.0 

-0.36S02E-02 


2 

4 

6 


0.0 

-0.348I3E-01 
-0. 13376E400 


0.0 

-0.96274E-01 

0.81952E-01 


0.0 

-0.54249E-02 
—0.23098 E— 02 


33 S3 = = = 33 333333B33333B3333Z33B33>333«33a>aaaa33333 


NEW AVAILABLE REAL ARRAY = 7473 


EIGEN VECTORS SCALED BY PAX. COMPONENTS 


VECTOR NO. ^ 

... —0.49 82 E— 03 

VECTOR NO. i 
-0.4632E-02 
VECTOR N0.= 

. -r-J3.1H5a£-01 


1 FREQUENCY = 0. 1 00863E+0 1 C/S 

0 .6857E- Q3 Q.2227E-01 -0.3075E-Q2 

2 FREQUENCY = 0.331536E + 01 C/S 

0.5255E-02 0.2572ET0O -0.1389E-01 

3 FREQUENCY = 0.608329E401 C/S 

Q.6Q63E-02 Q.lOOOEtOl Q.4682E-02 


0 • 97Q3E— 03 
0. 8978E-02 
0.3427E— 01 


0.111 400 — 0.1662E— Ql 0.5249E-02 0L. A QOUEtOl 

0.100 OE + 0 1 0.4398E— 01 0.1313E-01 -0.4307E400 

-0 .41 44ET00 -0 . 12S5E-Q1 -Q.5324E-Q1 -0.24J2E-Q1. 


*************** Bit********************** EIGEN ANALYSIS CONVERGED IN 2 ITERATIONS **************************************** 


EIGEN VALUES 

0.4016EF02 0.4339E+C3 0.1461E+04 ^ 

EIGEN VECTORS SCALED BY UMU 

VECTOR NO . = 1 FREQUENCY— 0 . 1 00863E *0 1 C/S 

— 0. 80 80E— 04 0.1 11 2E— 03 0.3612E-02 -0.4987E-03 0.1S87E-03 0.1808E-01 -0.2696E-02 0.B513E-03 0.1622E+00 


p 

M 

LC 

S — COMB I NED 

sx 

SY 

SXY 

SX-THERM 

SV-THERM 

SXY— THERM 

J 

... L 



-0 . 2894E— 03 

-0.6933E-08 

-0. 103SE-04 




3 

l 

1 


-0. 41 79E-03 

0.3200E— 02 

0. 2 29 7 E— 03 

—0. 1286E— 03 

0.3200E— 02 

0.2401E-03 

3 

2 

1 


-0 . 1396E— 03 

-0.1497E-03 

-0.2694E- 03 




3 

2 

l 


— 0 . 34 4 9 E— 03 

0.3127E-02 

— 0.3 75 9 E— 03 

—0 . 2053E— 03 

0.3277E-02 

— 0*8692 E— 04 

3 

3 



-0. 1498E-03 

-0. 1396E-03 

0. 2£94£— 03 




3 

3 

" l 


-0 . 1 1 49E— 03 

0.2897E-02 

0 • 3 76 0 E— 0 3 

0.3487E-04' 

~ 0.3037E— 02 

_ 0.8660E— 04 

3 

4 

l 


0.540 IE- 12 

-0.2894E-03 

0.9955E- 05 




3 

4 



-0.4200E— 04 

0 .2824E-02 

-0.2302E-03 

—0. 4200E— 04 

0.311 4E— 02 

— 0.2402E— 03 

3 

5 

1 


-0 .2894E-03 

—0.5001 E— 04 

-0. 9E03E-03 




3 

5 

1 


-0. 41 79E-03 

0.2721E-02 

-0. 8118E-03 

-0. 12B6E— 03 

0.2771E-02 

0* 1486E-03 

3 

6 

I 


0.3104E-03 

—0.64 9 BE— 03 

— 0 • 2 39 7 E— 03 




3 

6 

1 


-0.6340E-04 

0.2367E-02 

0.1 024E— 03 

-0.3738E-03 

0.3017E-02 

0.3421E-03 

3 

7 

1 


—0. 64 9 SE— 03 

0.31 0SE-03 

0. 2389E-03 




3 

7 

l 


— 0 . 875 1 E— 03 

0.31 79E-C2 

-0. 1 03 IE-03 

— 0. 2252E— 03 

0. 2 86 8 E— 02 

— 0*342 OE— 03 

3 

a 

1 


-0. 5006E-04 

— 0.2893E — 03 

0.9603 E— 03 




3 

a 

1 


-0.S207E-0J 

0.2824E-02 

0. B117E-03 

— 0.47 06 E— 03 

0.311 3E— 02 

—0 • 1486E— 03 

J 

9 

1 


-0.605aE-03 

0. 167OE-07 

— 0.7 57 l E— 03 




3 

9 

l 


-0. 8452E-C3 

0.3200E— 02 

-0. I 030E-02 

-0.2394E-03 

0.3200E-02 

-0.2733E-03 

3 

1 0 

l 


0. 7573E-04 

— 0 .60 l 5E — 03 

— 0. 6 056E— 03 




3 

1 0 



0 .71 74E-04 

0 .220 JE— 02 

-0.8029E-03 

-0.3988E-05 

0.2965E-02 

-0. 1973E-03 

3 

1 1 

1 


-0.6813E-03 

0. 755JE-04 

0. 606 1 E- 03 




3 

1 1 

1 


-0.9587E-03 

0.331 JE-02 

0. 8C36E- 03 

-0.2773E-03 

0.3238E— 02 

0.1975E-03 

3 

1 2 

1 


-0 . 1065E-08 

— 0.605BE— 03 

0. 757 1 E— 03 




3 

1 2 

1 


-0. 4200E-04 

0.239 7E — 02 

0. 1 030E-02 

—0. 4 2 00 E— 04 

0.3003E— 02 

0.2734E-03 

3 

l 3 

1 


-0.2597E-02 

-O.59Q5E-04 

-0. 1 432E- 02 




3 

1 J 

i 


-0.2776E-02 

0.2721E — 02 

-0.1 455E-02 

-0. 1790E-03 

0 .277 IE— 02 

-0.2319E-04 

3 

l 4 

1 


-0.6073E-03 

-0.2039E-02 

-0. 2547E- 02 




3 

1 4 

1 


-0.9205E— 03 

0.8662E-0J 

-0.2255E-02 

-0.3132E-03 

0.2906E— 02 

0.291 6E — 03 

3 

1 3 

1 


-0.2039E-02 

-0 .60 7JE-03 

0.2547 E— 02 




3 

l 'i 

1 


-0.237GE-02 

0.2 32 IE— 02 

0.2255E-02 

-0.3364E-03 

0.2929E-02 

— 0*291 6E— 03 

3 

1 o 

l 


-0 .5006E-04 

-0.2397E-02 

0. 1 432 E— 02 




3 

16 

I 


-0.5207E-03 

0.4664E-03 

0. 1455E-02 

-0.4706E-03 

0.3063E— 02 

0.231 9E— 04 


0.304057EF05 


1 0 
1 0 
1 1 
1 1 
1 2 
12 
13 
1 3 
1 4 
1 4 

15 
l 5 
1 6 

16 

1 
1 

2 
2 
3 

3 

4 

4 

5. 

5 

6 
6 

.1 

1 

2 

2 

. _3 ... 

3 

4 

4 

-.5 _ 

5 

6 
6 

Z___ 


._fl *.3040576+05 
0. 3070C8E+0S 
0*40566 08 + 05 
0.41 4035E+05 
- 0.* 135 4 4 18 1 05 
0* 1 l 82P96+05 
0 *5652 008+05 
0*5464 82 E ♦ 0 5 
0 • 26 Q525E+Q5 
0* 2746438+05 
0* 245984E+05 
0.2041 68E *-0 5 


3 7 

3 a 

3 8 

-3 A_. 

3 9 

3 1 O 

3 1 0 

-3- LJ 

3 1 l 

3 12 

3 12 

-3 -L3- 

3 13 


2 

2 

_2 

2 

2 

2 

.2 

2 

2 

2 


3 

3 

_3„ 

3 

3 

3 

6 

4 
6 
4 
6 
4 
6 
4 
6 
4 
6 


14 

14 

_ 1 - S _ 


1 5 
16 
16 

_1 . 

1 

2 

2 

J 

3 

4 

4 

5 

5 

6 
6 


2 

2 

2 

_2 

2 

2 

2 

_2 

2 

2 

2 

.2 

2 

2 

2 

_ 2 _. 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


0. 3040578+05 
0* 3078 12E+Q5 
0*4 Q5860E + 0 5 
0*4197958+05 
0 • 1354 4 IE 4-05 
0* 1 7909 2E + 0 5 
0.5852 00E +05 
0.6209948+05 
0* 260525E+05 
O.287068E+O5 
0.245984E+05 
0.2903478+05 


72 

0 • 
0* 
J • 
0. 
J. 
0* 
0. 
0 • 
0 • 
0* 
0. 
()• 
0* 
0 • 
0 . 


STRESS /STR 
9512E+00 0 

9366E+Q0 0 

95126+00 0 

89796+00 0 

901 4E ♦00 0 

88028+00 O 
67628+0 0 0 

7 22 8E + 00 0 

7544E *00 0 

9 86 4E f 00 0 

939QL 4-0 0 0 

93 9 9t 4-0 0 0 

9o£5E*00 C 
97308+00 0 

9/796+00 0 


AIN C3NSTRAI 
1 049E 0 1 0 

. 10138+01 0 

. 1 0498 + 01 0 

•11028+01 O 

* 1 099E + 0 1 0 

•11126+01 O 
•13246+01 0 

12778 + 01 0 

,66808+00 0 
992 3E + 00 0 

•99518+00 0 

9202E+00 0 

• 99 09E + 00 0 

.99076+00 0 

*9 75 1 E 4- 00 0 


-0* 0452E-CJ 
0. 75 73 E— 0 4 
0.71 74E-04 
-0. 601 38-03 
-0.9587E-03 
-0 . 1065E-06 
-0* 4200E-04 
-0.2597E-Q2 
-0.2 7766-02 
-0.6073E-03 
-0.9205E-03 
-Q. 20398-02 
-0.2376E-02 
-0 .5006E-04 
-0. 5207E-03 
-0 *41908+03 
0 « 5448E+03 
-0. 1983E+04 
0 .8328E+04 
- 0* 96086+04 
-O. 19168+04 
-0 • 3554E+ 05 
-0* 2843E+05 
-0 ■ 1978E+04 
0 .10276+05 
-0 • 2396E+05 
-0.2081E+05 
0 *2894E-Q3 
0. 16088-03 
0. 1396E-03 
-0 .65648-04 
Q-1498ErQ3 
0 .1847E-03 
-0.54018-12 
-0.42008-04 
0 .2094E-Q3 
0* 1608E— 03 
-0.3104E-03 
-0*68426-03 
Q.6499E-Q3 
0 .4247E-03 
0. 5006E-04 
-0.4206E-03 
0.6Q58E-Q3 
0* 3664E— 03 
-0*75738-04 
“0 • 797 2 E— 04 
0*681 3ErQ3 
0 * 4040E— 03 
0 • 1065E-08 
-0.4200E-04 
0.2S97E-02 
0.2418E-02 
0 • 6073E— 03 
0 • 294 1 E— 03 
0.2039E-02 
0. 1 703E—02 
O • 50 Q6E— 04 
-0* 4206E-03 
0*419 QE + Q3 
0* 1303E+O4 
0* 1983E+04 
0* 1229E+05 
0 . 96088+04 
O* 1730E+05 
0 • 3554E+ OS 
O* 4265E+05 
0. 1978E+04 
0.14238+C5 
0 • 23968+05 
0 • 27 1 OE+05 


0* J2U08-02 

— 0.601 58— 03 
O.220JE-O2 
0* 7553E-04 
0.331 JE-02 

-0.6058E-03 
0. 2397E— 02 
-0.50Q5E-04 
0.2/21E-02 
-0.2039E-02 
0 • 0662E— 0 J 
-0.6Q73E-03 
0*232 IE- 02 

— 0 *2597E— 02 
0* 4664E — 03 
0*0 

-0.3969E+04 

0.0 

0*85178 4-04 
0*0 

O.S016E+O4 

0*0 

0*5 723E 4-03 
0*0 

0*570484-04 

0*0 

— 0*26008 4-04 
0*69538-00 
0.32008-02 
0*1 4978— 03 
0*3427E— 02 
0*13968-03 
0*31 76E- 02 
0*28946—03 
0.3403E-02 
0*5001 E— 04 
0*2 821 E — 02 
0*64986-03 
0.3666E-02 

-Q.3105E-03 
0*25586-02 
0 *28936 — 03 
0.3403E-02 
-0* 1678E-07 
0*32006—02 
0.6815E-03 
0 *3646E— 02 
-0.75538-04 
0*31626-02 
0*6058E— 03 
0.3608E-02 
0*50058-04 
0*282 1 8—02 
0.2039E-02 
0*4945E-02 
Q *60738.-03 
0 *3536E— 02 
0.2597E-02 
0*566 0E-02 
0.0 

—0*39696+04 

0*0 

0*651 7E ♦ 04 
0*0 

0.5316E4-04 

0*0 

0*5723E 4-03 
0*0 

0*57046404 

0*0 

-0*26008404 


-U. 1 UJOt-D<! 
-0. 60SOE- 03 

— O .2394 fc— U 3 

u 

— U •£( JJC'UJ 

-o. aoa9e-03 
0.6 06 IE- 03 

-0. 3988E— 05 

0.296 5E— 02 

-0. 1V73E-03 

0. 8C36E-03 
0.7571 E- 03 

-0.2773E-03 

0.3238E-02 

0.1 975 E-03 

0. 1 030E-02 
-0. 1 432E— 02 

— 0.4200E— 04 

0.3003E— 02 

0.2734 E— 03 

-0. 1455E-02 
— 0. 2547E- 02 

-0. 1790E-03 

0.2771E-02 

-0.23I9E-04 

-0.2255E-02 
0 . 2547 E— 02 

-0.3132E-03 

0.2906E— 02 

0.2916E-03 

0.22S5E-02 

0.1432E-02 

—0. 3364E— 03 

0.2929E-02 

-0.2916E-03 

0. 1 455 E— 02 

-0.4706E-03 

0.3063E-02 

0.2319E— 04 

0. 1755E*05 



0. 1 755E*05 
0. 2340E»05 

0.9638E403 

-0.3969E+04 

o 

. 

o 

0.2240E«-05 

0.5512E+0* 

0.1 031E40S 

0.8517E404 

0.0 

0. 5E12E40* 
0.26U4E+05 

0.7691E+04 

O.S816E+04 

o 

. 

o 

0.26B4E405 

0.7109E404 

0.5723E+03 

o 

• 

o 

-0.150QE40S 




-0. I500E+05 
0.3227E+04 

0. 1225E*05 

0.5704E404 

0.0 

0.3 22 7E* 04 
0. 10JSE-04 

0.3148E404 

— 0.260 0E4 04 

o 

. 

o 

0. 2504E-03 
0. 2694E- 03 

-0. 1286E-03 

0.3200E— 02 

0.2401 E-03 

0.2029E-03 

-0.2053E-03 

0 . 3277E— 02 

-0.S652E-04 

-Q.2e94frr.QX_ 



— 0.2 02 0E— 03 
-0.9955E-05 

0.3487E-04 

0.3037E-02" 

0.666 0E -04 

-0.2502E-03 

0.9CQ3E-03 

-0.4200E-04 

0.3U4E-02 

-0.2402E-03 

0. 1 IQ9E-02 
0.2397E-03 

— 0 . 1286E— 03 

0.2771E-02 

0.14 86 E-03 

0.561 9E— 03 

n . oiane-nt 

-0.3738E-03 

0.3017E-02 

0.3421E-03 


-O 

-0 

-o 

Q 

0 

O 

0 

-Q 

-0 

-0 

-O 

O 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-o 

-0 

-o 

-0 

-0 

-o 

-0 

o 

0 

-o 

-0 


5E09E-03 

9603 E— 03 

l 109E-02 

757 lEr 03 

48376-03 

60568- 03 

4083E-03 

6 06 1 E— £L3_ 

4086E-03 

75718-03 

483 7 E- 03 

X432E-02 

1 4096—02 

2547E-02 

2838E-02 

2547Er-Q2_ 

2838E-02 

1432E-02 

I409E-02 

1 75564 05 
1 755E405 
2340 E4 05 
2340E4 05 
5S12E404 
55126404 
2684 E4 05 

2 C84 E4 05 
1 50064 05 
I500E405 
3227E404 
3227E404 


-0.2252E-03 0.2868E-02 -0.3420E-03 


-0.4706E-03 

0.31 13E—02 

-0.1 486 E-03 

— 0 . 2394E— 03 

0.3200E-02 

-0.2733E-03 

— 0.3986E— 05 

0.296SE-02 

-0.19736-03 

-0. 2773E-03 

0.3238E-02 

0 .19756— 03 

— 0.420 0E— 04 

0.3003E-02 

0.2734E-03 

-0. 1790E-03 

0.277 IE* 02 

—6.23196—04 

-0.3132E-03 

0.2906E— 02 

0.2916E— 03 

-0.3364E-03 

0.2929E-02 

-0.29166-03 

-0.4706E-03 

0.3063E-02 

0.23196-04 

0.9638E403 

-0.3969E404 

o.o 

0.1031E405 

0.8517E404 

O 

. 

o 

0.7691E404 

0.5816E404 

o 

. 

o 

0.7 109E404 

0.5723E403 

o 

• 

o 

0. 122SE40S 

0.5704E404 

0.0 

0.31 48E4 04 

— 0.260 0E+ 04 

o 

• 

o 


NTS 

• l 1018401 

•1 164E401 

• 1 15484 Cl 

•1 1808401 

• l 1018401 

• 1 I Udfc4 0 1 

•1 1548401 

• 1 13284 01 

• 9054E4 00 
l 0088401 
l 005E401 
10806401 
1 0498401 
1009E4 0 1 
1 0058401 


0*32006400 
0.38448400 
0.42178400 
0*32468400 
0 *32 0 064 00 
0.2959E400 
0.421 7E400 
0.50676400 
0*56286400 
0*1 194E401 
0* 11936401 
0*12086401 
0* U93E401 
0. 1207E40 l 
0*12046401 


TO 172 
0*98766400 
0*97968400 
0*95606400 
0*99446400 
0.9442E400 
0 *95656 400 
0 .9212E400 
0*87796400 
0*78036400 
O.2710E4OO 
0.27696400 
0*22096400 
0.27696400 
0.22526400 
0.23326400 


MOST CRIT 
0 • 9590 E 40 0 
0.99516400 
0*9926E400 
0*93936400 
0* l OO0E4O1 
0*9951 E400 
0 • 6926E 40 0 
0*93936400 
0*83676400 
0*90906400 
0* 9065 E 40 0 
0.9685E400 
0.9399E400 
0.97796400 
0*97306400 


ICAL CONSTRAINT 
0. 1040E401 0 

0*10056401 0 

0.1007E401 0 

0.1061E401 0 

0*99166400 O 
0*10056401 O 
0*11076401 0 

0*10618401 0 

0*10196401 0 

0*10226401 0 

0*10196401 0 

0*10506401 0 

0*10436401 O 
0*1 047E401 O 
0* 1282E4Q 1 0 


-0*20265356400 


•1 1 77E401 

.1 160E40I 

•1 134E401 

• 1 1 6064 01 

• 1 1 2964 o I 
.11366401 
•1049E401 
•1026E4Q1 
*9812 £400 
•978SE400 
•98I2E+00 

9505E400 
9572E+00 
9529E400 
•71 796400 


O *335564 00 
0* 3999E40 0 
0.497IE400 
0*39996400 
0.51496400 
0.49076400 
0.8159E400 
0.90096400 
0.1181 £401 
0*11806401 
0.11606401 
0.11456401 
0* 1181640 1 
0.11796401 
0.11606401 


0.97966400 
0. 98756*00 
0.99456400 
0.9560E400 
0.95656400 
0.9442E400 
0 .87796400 
0.9212E400 
0.32006400 
0.32506400 
0.40056400 
0.45656400 
0*32006400 
0.3280E400 
0.40056400 


i3d 



1 


1 72 
J. 
Cl. 

J • 

0. 
J . 
0 . 
0. 
0 . 
0. 
o. 
Q. 
()• 
0. 
0 • 
J • 
a • 
a . 
o. 



J 

2 

0. 1 7989 2E 40 5 

0. 1 7 3 OE 4 0 5 

0. 5816E4 04 

-0. 551 2E4 04 

0 • 769 1 b 4 0 A 

0 .5816E404 

4 

4 

2 

0. 585200b 405 

0 .3554E4 05 

0.0 

-0. 2684E4Q5 



ft 

4 

2 

0.O28994L4O5 

0. 4265E405 

0.5/23E 403 

-0. 2684E+ 05 

0.7109E404 

0.5723E403 

4 

5 

2 

0.260525E40S 

0 • 1 978E404 

0.0 

0. 1 500E4 05 



t> 

3 

2 

0.i>8?at>rtE + U5 

0 . 1 42 3b+ C5 

0.5 704b 404 

0. 1 500E405 

0. 1225E405 

0.5704E404 

4 

ft 

*> 

0.2459H4L405 

0.2396E4Q5 

0. 0 

-0. 3227E404 



6 

6 

2 

0.29034 7t 405 

0.2710E405 - 

0.260 Ob 404 

-0. 32276+04 

0.31 48E 4 04 

—0.26006404 


ST=L SS/S TRAIN 

constraints 

1 TO 

172 MOST CRITICAL 

CONSTRAINTS -0. 

2026535E400 


95 1 2t + 00 
936 6t 4-00 
9512btQQ 
0 9 7ML + 00 
90 1 4E.+0 0 
0ri82L*OO 
O762L+00 
72206+00 
7544E *00 
986 4E 4-00 
9890L tOO. 
93 9 9b 4 00 
9o£ 5E 40 0 
9738b 40 0 
9779b tQO 
y23 7b 40 0 
8463b 40 0 
76S 7b400 


0.1 0 4 9 E 4 0 1 
0. 131 JE+01 
0 • 1049E4Q 1 
0. 1 102b 4-01 
0 • 1 099E+ 01 
0. I 1 12E401 
0. I 324t>01 
0.12 77b 4- 01 
0 .6fc>BOE400 
0.9923E+00 
0.9951£tQQ 
0 .9202E4-00 
0.950 9E 400 
0.990 7E 4 0 0 
0 .99516+00 
0. 1 034E 4 01 
0 . 9509E 4-00 
0 • 7677E+Q0 


0 • 1 1 8 l b 4 0 1 
0.1 164E+01 
0.1 1S4E4 Cl 
0.1 1 80t 4- 0 1 

o • 1 1 ti l b+ o i 

O. 1 IU6E4 01 

0.1 lb4E+01 
0.1 132E+01 
O .9054E4 00 
0.1008E+01 
0.100 5E 4-01 
0 . I 080.E + 0 1 
0.1 049E4-01 
0. 1 U09E+ 0 1 
0 .1 005L+QI 
0.9657E+ 00 
0 • 1049E4-01 


0.3200L+C0 
0. J844E400 
0.421 7E+Q0 
0.3246E+00 
0.3200b* OO 
0.2959E+00 
0.421 7E + 0Q 
0.5067L+00 
0. 5628L+00 
0.1 194E401 
0. 1193E+01 
0. 1208L+01 
0.1 1 93E401 
0. 1207E+0 l 
0 . 1204E+0 1 
0. 1280L+01 
0. 1321E401 


0 .96 76E400 
0 . 9 796E 4-00 
0. 9S60E+00 
0.9944b 400 
0 .9442E400 
0. 9565E400 
0.921 2E400 
0 .87 79b 400 
0. 780JE+00 
0 .2718E400 
0.276 9b 4 00 
0 . 22 09b 4 00 
0*2 769b 40 0 
0 . 2252E4 00 
0.2332E+00 
-0. 5079E-01 
-0.2027E400 


0.9598E400 
0.9951 E+OO 
0 . 99 26 E 40 0 
0.939JE400 
0.10 08E 40 1 
0*9951 E 400 
0.8926E400 
0.939JE400 
0 .8367E 40 0 
0.98 90 E 40 0 
0.9865E400 
0.9685E400 
0.9399E400 
0.97 79 E 400 
0 .97 3d E 40 0 
0 .84 63 E 40 0 
0.9237E400 


1 005t4 0 1 0.1160E401 
1007E401 0.U34E40I 
1061E401 0.1 160E401 
9916E400 O.I129E401 
1005E401 0.11 36E401 
U07E401 0.1049E401 
1061E401 0.1026E401 
1019E401 0.9812E400 
1022E401 0.9785E+00 
1019E401 0.9812E+00 
1050E401 0.9505E400 
1 043E40 1 0.9572E400 
1047E401 0.9529E400 
I282E+01 0.71 796+00 
1199E401 0.8013E400 
7538E400 0 .66426400 


0 .0 
0.0 
0.0 
0*0 


0.3355E400 
0 • 3999E400 
0.4971E400 
0.3999E400 
0.5 149E400 
0.4907E400 
0 . 8 1 59E4 00 
0.9009E400 
0.1 181 £4Q| 
O.110OE4O1 
0.1 160E40 I 
0.1 145E401 
0. 1 I01E4O1 
0.1 179E4 01 
0.U6QE+01 
O. 1200E401 
0.8561E400 


0 -9796E400 
0. 90756*00 
0.9945E+00 
0 .95606400 
0.9565E400 
0.9442E400 
0.8779E400 
0.9212E400 
0.3200E+00 

0.3250E400 

0.40056400 
0.45656400 
0.32006400 
0.32006400 
0 *40056400 
0.24066400 
0.49666400 


l FREQUENCY CONSTRAINTS 
0 • 4070E-02 


173 TC 173 


HOST CRITICAL CONSTRAINT* 0.40698906-02 


0 

CONSTRAINTS 

OUT 

OF 

0 

CUTOFF 

POINT* 

0.0 






0 

CQNSTRAI NTS 

OUT 

OF 

0 

CUTOFF 

POINT* 

0.0 






.50. 

CQtiSTRAI NTS 

OUT 

OF 

172 

CUTQFF 

POINT* 

0.8797356400 







4 

9 



14 

19 

24 

29 

34 

39 

44 

49 


54 

59 



6 1 

64 

69 

70 

71 

74 

75 

61 


82 

84 



85 

86 

90 

95 

100 

1 05 

t to 

115 


L2Q 

125 



130 

135 

140 

145 

148 

150 

155 

156 


158 

160 



161 

165 

167 

168 

169 

170 

171 

172 

31 

CONSTRAINTS 

OUT 

OF 

50 

RETAINED 

DUE TO 

VARIABLE LINKING 









9 



19 

39 

44 

54 

6| 

64 

71 _ 

74 


75 

81 



82 

84 

85 

90 

95 

1 05 

115 

130 


145 

1 48 



1 50 

155 

150 

161 

165 

167 

168 

170 


172 












1 

CONSTRAINTS 

173 

OUT 

OF 

1 

CUTOFF 

POINT* 

0.900407E400 







CONSTHAl N T S 

OUT 

OF 

0 

CUTOFF 

POINT* 

0.900407E400 







AVAILABLE INTEGER ARRAY = 2500 OVERLAY ANALYS REQUIREMENT 3 812 


POSTURE TA9LE 


RETAINED 

TOTAL 

type 

MEMBER 

NODE DIRECTION 

L.C. 

MODE 

CONSTRAINT VALUES 

TRESS/STRAIN 

1 

CONSTRAINTS 

4 

MOST 

3 

CRITICAL 

1 

= — 0 .20 26 536+0 C 

1 

3 

0.31 9969E+00 

2 

9 

3 

2 


1 

3 

0 .3354076+00 

3 

19 

3 

4 



3 

0 *399656 E 4 00 

4 

39 

3 

8 


1 

3 

0.399S80E4 00 

5 

44 

3 

9 


l 

3 

0.3I9997E400 

6 

54 

3 

I 1 


1 

3 

0.295687E400 

7 

61 

3 

13 


i 

-2 

0.676I79E400 

8 

64 

3 

13 


1 

3 

0*421 723 £400 

9 

71 

3 

15 


1 

-2 

0. 722825E+00 

l 0 

74 

3 

15 


1 

3 

0 » 306687E400 

1 1 

75 

3 

15 


1 

4 

0 • 87 7 662 E 400 

12 

61 

4 

1 


1 

1 

0.754393E400 

13 

82 

4 

2 


1 

1 

0.668771E400 

14 

84 

4 

4 



1 

0.562014E+00 

1 5 

83 

4 

5 . 


l 

i 

0* 78 026 5E 400 

16 

90 

3 

1 


2 

3 

0*31 99 66 E 400 

1 7 

95 

3 

2 


2 

3 

0*271 844 £400 

1 8 

1 05 

3 

4 


2 

3 

0 .2768676+00 

19 

115 

3 

6 


2 

3 

0 . 22088 7E +00 

20 

130 

3 

9 


2 

3 

0.32 00 04 £+00 

21 

145 

3 

12 


2 

3 

0.233213E400 

22 

148 

3 

13 


2 

2 

0.71 7949E+00 

23 

150 

3 

13 


? 


A - .Art.COC a . Art 



MODE 


AVAILABLE REAL ARRAY 


<+ 

jy 

J 

8 

1 

3 

5 

44 

3 

9 

1 

3 

6 

54 

3 

1 1 

1 

3 

7 

61 

3 

13 

i 

-2 

8 

64 

3 

13 

1 

3 

9 

7 1 

3 

15 

1 

-2 

1 0 

74 

3 

15 

l 

3 

1 1 

/ 5 

3 

IS 

1 

4 

12 

81 

4 

1 

1 

i 

13 

82 

4 

2 

1 

1 

14 

84 

4 

4 

1 

1 

15 

85 

4 

5 

i 

i 

16 

90 

3 

1 

2 

3 

1 7 

95 

3 

2 

2 

3 

1 8 

1 05 

3 

4 

2 

3 

19 

11S 

3 

6 

2 

3 

20 

130 

3 

9 

2 

3 

21 

145 

3 

12 

2 

3 

22 

148 

3 

13 

2 

2 

23 ... 

.... ISO 

3 

13 

2 

3 

24 

155 

3 

14 

2 

3 

25 

158 

3 

15 

2 

2 

26 

16 1 

3 

15 

2 

4 

27 

1 65 

3 

16 

2 

3 

28 

167 

4 

1 

2 

i 

29 

168 

4 

2 

2 

i 

JO 

170 

4 

4 

2 

i 

ji 

172 

4 

6 

- 2 . 

. .. _ A 

CY CONSTRAINTS 
32 173 

MOST CRITICAL * 0.406990E-02 


-l 

stands for the following 

NEGAT IVE’LCwER BOUNO POS I TI VE=UPPER BOUND 

FOR STRESS CONSTRAINT. (CODE+1) 

-1_= YCN.MISE.S EQUI VALENT . STRESS 



2 = LONGITUDINAL STRAIN 

3 = TRANSVERSE STRAIN 

4 * SHEAR STRAIN 

5 ? FIRST EQUATION OF STRESS INTERACTION 

6 = SECOND EOUATION OF STRESS INTERACTION 

7 = TSAI-AZZI CRITERION 

FOR FREQUENCY CONSTRAINTS. ASSOCIATED MODE NUMBER 




7500 

OVERLAY ANALYS 

REQUIREMENTS 1000 




0*3998566+00 
0*3998806+00 
0.3 1999764-00 
0 • 295887E+Q0 
0*6761 79E+Q0 
0 *4217236+00 
0* 722825E+00 
0* 5066876+00 
0 • 877862 E tOO 
0* 754393E+00 
0 *66877 IE *00 
0*56281 46*00 
0.7a0285E+00 
0*3199666+00 
0*271 844 E+00 
0*27686764’ 00 
. 0*220aa7Et00 
0*3200046400 
0*2332136400 
0*71 79496400 
0*4004526+00. 
-0*5079206-01 
0*80 1 324E400 
0 • 846270E400 
-0* 2026 S3E tOO 
0*7537506400 
0*6641646400 
0*4968046400 


0*4069906-02 


ANALYSIS TIME DATA 

ASSEMBLE MASS/STIFFNESS MATRIX 

ASSEMBLE LOAD VEC TOR S 

DECOMPOSE STIFFNESS MATRIX 
SOLUTION CF DISPLACEMENTS 
FREQUENCY ANALYSIS 

-ELUT-TEfi -ANALYSIS 

CONSTRAINT EVALUATION 

POSTURE TABLE SET 

SELECTIVE GRAOIENT EVALUATION 


0*98541 3E— 01 
Q.519714E-QI 
0.239 563E—02 
0.271 606E-02 
0.32394 4E — 0 I 
0.0 

Q. 160202E400 
0.601 044E-0 I 
0.486 038E400 


GRAND TOTAL CPU TIME 


0.9 7 593 7 E 400 


ENO OVERLAY ANALYS 

SCALING FACTOR 0*1202656401 SCALED WEIGHT 0.S19228E405 

SIDE CONSTRAINTS 


RELATIVE MCVE LIMIT 0*2000E400 


VARIABLE 
BOMBER 

LOWER 

BOUND 

ACTUAL 

size 

UPPER 

BOUNO 

VARIABLE 

NUMBER 

LOWER 

BOUND 

ACTUAL 

SIZE 

UPPER 

BOUND 

01 S* 

1 

o.aoooE-oi 

0*40006+00 

0 • 2000E + 0 1 

DIS* 

2 

0*60006-01 

0*35006+00 

0* 1760E+0 1 

. DIS* 

_3 . 

0*20 006.^01 

0*10006+00 

0*5 100E + 0 0 

DIS* 

4 

0*25006+00 

0*12506+01 . 

Q.20Q0E+Q1- 

DI 5. 

5 

0. 1 5006+00 

0*75006+00 

0* 2000E+0 1 

DIS* 

6 

0*50006-01 

0*25006+00 

0.1260E+0I 

01 S. 

7 

0.4000E-01 

0 • 200 0E+ 0 0 

0* 1 010E + 0 1 

DIS* 

a 

0*20006-01 

0*10006+00 

0.5 1 OOE+OO 

DI S. 

9 

0. 20006-01 

0.1 OOOE+OO 

0*5 100E + 00 


10 

0*12006+00 

0*60006+00 

0.3000E401 


l_l 

Q . 1 20 OE 1 00 

0*600 OE + QO 

0* J000E+Q1 


12 

O*12Q0£+QO 

0*60006+00 

Q.30QQEtQl 


13 

0. 12006+00 

0 • 600 0E + 00 

0* 3000E + 0 1 







MOST — VIOLATED 5 inE CONSTRAINT 


DESIGN VARIABLE 


4 


rnwcTo* » .if 




t 


I 




di/iuco nc iUHf 0 • 51 922 8t +0 5 


SIDE CONSTRAINTS 


RELATIVE MCVE LIMIT 0.2000E+00 


VAR I ABLE 
NUMBER 

LOWER 

BOUND 

ACTUAL 

SIZE 

UPPER 

HOUND 

VAR I ABLE 
NUMBER 

LOWER 

BOUND 

ACTUAL 

SIZE 

upper 

BOUND 

BIS* 

1 

0.8000E-01 

0 *40006+00 

0 .20006+0 1 

BIS* 

2 

0 *60006—0 1 

0.3500E+00 

0*1 760E+0 I 

BIS. 

3 

0*20006-01 

0. 1 00 0E + 0 0 

0.51QQE+00 

DIS. 

4 

0 *2 5 GOE+ OQ 

0* 1250E+01 

0.2000E+01 

DI 5. 

5 

0. 1 500E+00 

0. 7500E+00 

0* 2000E + 0 l 

BIS* 

6 

0 .5000E-01 

0 • 250 0E + 0 0 

0* I260E+01 

01 s. 

7 

0.4000E-01 

0*20006+00 

0. 1 010E+0 I 

DIS. 

8 

0 .2000E-OI 

0. 1000E+00 

0*5 1 OOE+QO 

DI S. 

9 

0.2000E-01 

0. I 000E*00 

0*5 100E + 00 


10 

0*1 200E+00 

0* 60 0 0E +0 0 

0 • 3000E+0 1 


1 1 

0.1200E+Q0 

Q .600 OE + QO 

0 .JOOOE + U I 


12 

0 .12GQE+00 

0*600 OE tOO 

0.3 QQGE+Q1 


13 

0. 1200E+Q0 

0*60006+00 

0* 3000E+0 1 







MOST. VT OLATEU S IDE CONSTRAINT - DESIGN VARIABLE 4 


CONSTRAINT VALUE O.AOQOE+OO 


ENTER OVERLAY PREDL1 


AVAILABLE REAL ARRAY = 7500 OVLPLAY PPEOU 1 REQUIREMENTS 4796 

AVAILABLE ..IN TEQEJ3 ARRAYS 2500 OVERLAY PREDU1 REQUIREMENT? 770 


OJ 

OJ 


===== DUAL 


OPTIMIZER == * = 




.DUAL VARIABLES : 

O 67 9 4E 4-0 4 0*8d794F>04 

0*«a79*#E + 04 0.88794L4-04 

0.88 J94E+04 0.aa794E*04 

0.dfl794e>04 0.8B794L4-04 


0. 88794E+04 
0 *8B 794E >04 
0.88794E+04 
0-8Q794E4-04 


0* 88794E4-04 
0.88794E4-04 

o. ea794r:«-04 

0.88794E*04 


0 « 88 794E*04 
0 •B87Q4E+04 
O.B0794E+O 4 
0 • 88 754E ♦O 4 


0.88794E+04 

0.88794E+04 

0.88794E4-04 

0.88/94E+04 


0*807946+04 
0.8Q794E+04 
0* 88794E+04 
0*887946+04 


O • 88 7 94 E+ 04 
O *887946+04 
0.88794E+04 
0 *887946+04 


phase com dual obj. 


WEI GHT 


ZMCD 


nact nplan iconj 


NFWTflN Kin f c 


t nr c 


NALFA 



===== U U A L I 


QUA! VAR! AflLE S ; 

0 .88794E404 
0.88794E404. 
0.86794E404 
iUHB 29 A £ t 04 


0.88794E«-04 
0 • 88 7 94E404 
0.88794EF04 
. .ft.a3Z-94E + 04 


0.88794E404 
0 . 88794E *04 
0.88794E404 
Q > 88794E4Q4 


0.88794E404 
0* 88794E404 
0.88794E404 
0.88794EtQ4 


PHASE 0 DM 


OUAL OBJ. 


WEI GHT 


72QE + 0& 

-0. 1477E405 
-O. 1098ET05 
— 0*9 786E4 04 
. _-Q. 8fc54£tQ4 
-0.7762E404 
-0.6301E404 
-0.6285E+04 

-a .J560£<-04 

-0.301 6EA04 
-0 .2625E+04 
0.2447E404 

Q.J797fct04. 

0.960 IE+04 
0.1006E405 
0. 1289E.05 

0.1562E+05 
0. 1864EF05 
0.197 7E40S 

aji . s 77 et . QS . ... 

0. 1977E+05 
0.1977E405 
0.1993E405 
0 .2062E40S 
0.21 66E405 
0.2232E405 
O .2287E*05 

0.2:iL2Et 05 . .... 

0.2326E+05 
0 .2338E405 
0.2341E405 


Q.2313EJ05 

0.2277E+05 
0.2248E405 
0.2243E405 
0.221 1E4Q5 
0.2206E405 
0.21 83 E + 05 
0.2 183E4 0 5 
.. Q,2.134E«-05 
0.2133E+05 
0.2 132E+ 0 5 
0 .2074E405 
. 0.207 5E* 05 

0.2086EF05 
0.2086E405 
0.2085E405 

OjtgQ & SEi-QS 

0.2083EJ0S 
0.20B9E405 
0 • 2088 E+05 
0 . 2088E+Q5 
0.2086E+05 
0.2088E405 
0.2093E405 
0.2 10SE»05^ 
0. 2154 E* 05 
0.2167£«-05 
0.2240E405 
.. Q.23Q4E+Q5 
0.2293EJ05 
0.241 3E* 0 5 
0.2354E405 
0. 23S4E + 05 


0.347JE+01 
0 .302 J E 4 0 1 
0 . 2882E + 0 1 
0 .2696E + 01 
0 • 2536E + 0 1 
0 • 235SE+ 0 1 
0 .21 79E + 0 1 
0.2057E +01 
0 .203 1 E + 0 1 
0 • 1899E + 0 1 
0. 1 7 36E + 0 1 
0.1564E401 
0 .1 432E + 0 1 
0. 1412E + 01 
0. 1302E + 01 
0. 12 1 OE+O 1 
0, 1091E+0} 
6 .9739E + 00 
0 • 6832E+00 
0 . 80 8 5E + 0 0 
0.7424E+00 
0.6699E+00 
0 . 5886E+00 
0 .494 1 E 4 00 
0.38 7 1 E 400 
o . 25 aaE + oo ' 
0. 1773E+00 
0. 1331E400 
0.8177E-01 
0.8710E-0I 
0.3076E-01 
0 .2632E— 01 
. _Q m.Q. 


DUAL VARIABLES 

0.0 


0.0 0.0 

0.0 0.0 

0.0 0.0 

0 . 6834 IE 404 0.0 


0 .2341E405 
0.2341E*05 


_ WEIGHT 

0 .235 4E*- 05 
0 . 2354E + 05 


0.2532E-01 

0.0 


NUMBER OF RESTARTS 
TOTAL NUMBER OF O.C.M. 

NUMBER 0 F N ON— ZERO DUAL VARIABLES 
NUMBER OF DISCONTINUITY PLANES 


DUAL OBJECTIVE FUNC T ION 

NORM OF (PROJECTED) GRADIENT 

DUAL VARIABLES : 

0.0 O.C 

0.0 O.C 

0.0 ... O.C 

0.0 O.C 


0 .234 06 0E+ 05 


0.0 0.0 

0.0 0.0 

0.0 0.0 

0.68341EF04 0.0 


ANALYSIS OF DUAL SCLU"ION 


NUMBER ?F UlSC3NriNUI T Y PLANtS 
NUMBER OF POSSIBLE PRIMAL POINTS 


WEIGH" OF UPPER BOUND SOLUTION 


M I Z E R ===== 


0.88794E+04 0.88794E+04 0.88794E+04 0.88794E+04 
0.8B794E+04 0.88794E+04 0.88794E404 0.88794E404 
0.88794E+04 0.88794E+04 0. 88794EA04 0.88794E404 
0.8Q794E+04 0. 887946404 0.887 94640 4 0.8S794E404 



0.0 0.0 

0. 0 0.0 

0.0 0.0 

0.0 0.0 



= 100 


= 32 

= 13 

= O.IOCOOCE — 04 

= 0 .1 OOOOOE— 03 


0.0 0.0 

o . a . o.o 

0.0 0.0 

0.0 0.24892E405 


0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0 . 24692 E 405 



I 

NUMUth l'*- UISlON” INUI t Y PLANfcS 

DUAL OBJECTIVE FUNCTION 

NORM OF (PROJECTED) G D AD ILN t 

DUAL VARI ABLES : 

0.0 0.0 

0.0 0.0 

... 0.0 0 • 0 

0.0 0.0 

ANALYSIS IF DUAL SOLUTION 


0 .234060E405 

0.0 


0.0 

0.0 

0.0 

0.6634 IE+04 


NLD V 

EPSPH 

EPSODM 


0.0 

0.0 

0.0 

0.0 


13 

O.l 0C00GE-04 
0.1 000006-03 


0.0 
0. 0 
0. 0 
0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0*246926405 


NUMBER OF D ISCD NT I NU I T Y PLANES 
NUMBER OF POSSIBLE PRIMAL POINT S- 


WEIGHT OF UPPER BOUND SOLUTION 


2 

4 

0*2354356405 


INFEASIBLE 


MOST VULATtU CCNSTRAINT 


32 

VALUE = 

0*3860536-02 



PRIMAL VARIABLES 

0 • 2 2500E + 0 0 0.17M3E400 

0.20000E400 

0. 47200E400 

0.506C7E400 

0*200006400 

0 • 33500640 1 

0.18000E+0I 

. tuaaooaEtoa 

0^913396400 

0*907576400 

0 * 125Q4E40 1 

0 *63 7 60E 40 Q 



— • 

CCNSTRA INTS 

-0.3 1 989E400 

-0 .29248E400 

-0 .435186400 

-0 .435286400 

-0*31 999E40 0 

-0*259846400 

-0*892556400 

-0.565796400 

-0..7653 JE4Q0 

-0*333 71 E4QQ 

-0 . 10150E40I 

—0*695786400 

— 0 • 49 5 57E 40 Q 

—0*780276400 

-Q.46765£t00 

-0.31990E400 

-0. 14412E400 

-0.11 0 14E400 

—0 • 62869E— 0 1 

— 0*320016400 

— 0* 92 5 40 E—0 1 

-0*847916400 

-0*466586400 

-0.44490E-0I 

-0.85364E400 

-0.88975E400 

— 0*34 064 E—0 1 

-0 .6SS67E + 00 

-0.50304E400 

-0.67I07E400 

-Q.7816IE*00 

0.3S605E-02 

Mii_LCMX-nF- .U3BEH flCi-NQ SOLUTION 

0 *23382164 OS 


-- . . . _ 

- - - 

. . 

PRIMAL VARIABLES 

0.22500E400 0.17143E400 

0*200006400 

0*464006400 

0. 5066 7E 40 0 

0*200006400 

0*335006401 

O.ISOOOEtOI 

... . . .... - Q.2B00QE + QQ 

0* 91 J39E40 0 

0* 9075 7E4QQ 

0* 12 5Q4E + 0 1 

Q.6376QE4Q0 




constraints 

-0.3 1988E400 

-0.2881 7E400 

—0 • 4339 BE 40 0 

-0 *43409E400 

-0 *319996400 

-0*256916400 

-0.89527E400 

-0*578686400 

-0.760 43E.tQQ . 

-Q.333Q8Et00 

-0.1Q196Et01 

-Q.69l77EtQ0 

-0.A9348E+00 

-0.76374E+Q0 

-Oa*6699E400 

..-Q»3IS9#EtOJL 

-0 . 1 3932E4 0 0 

-0. 10704E+00 

-0 .63235E-01 

— 0*320016400 

-0.91 815E-01 

-0.814866400 

-0*476996400 

— 0 • 16 078 E—0 1 

-0.83804E400 

-0.871 75E400 

0.26137E-01 

— 0*691 65E 40 0 

-0.50 l 09E400 

-0*660326400 

-0.77469E400 

-0.62685E-02 


ttEJ&HX-OF B INA1 Ot SIGN 


0 • 235435E4 QS 


MOST VIOLATED CONSTRAINT 32 

- NLIMHER QF— gRXMAl VARIABLES FROM LOWER BOUND SOLUTION 


INFEASIBLE 
VALUE = 

0 


0.3 6 6053E-02 


PRIMAL VARIABLES 

0.2 2500E+ 0 0 0.171436400 

... Q.3QQQ QEtQO Q.91339E4Q0 


0 • 20 00 OE 4- 0 0 
0 . 9075 7E tOO 


CONSTRAINTS 
-0.3 1989E400 

tQ *7651 3£t Q Q 

-0.14412E400 
-0 .8536 4E400 


— 0 .29248E40 0 -0.43518E400 

.-Q*J3371EtQQ -0.10150E4QI 

— 0. 1 i 0 14E4Q0 -0.62869E-01 

— 0.88975E4QO -0.34064E-01 


0 .47200E *00 
Q* 125Q4E401 


-0.43528E400 
— Q.69578E+QQ 
-0.32001E400 
-0.69567E400 


0.50667E tOO 
0.637EQE+QQ 


— 0.31 999E40 0 
-0.49557E4QQ 
-0.925396-01 
— 0.50 3 04 E 40 0 


0 . 200 0 0E + 00 


— 0 • 259 8 4E 4 0 0 

— Q •_ 78 0 2 7E4 00 

— 0 . 84 79 IE 400 
-0 .671 07E400 


0. 33500E40 1 


— 0. 8925 SE 400 
-Q*46765£tOQ 
-0.46658E400 
— 0 . 78 16 IE 400 


0*180006401 


-0 *5657964 00 
-0*319906400 
-0 *444906— 0 I 
0*386056-02 


END OVERLAY PREDU1 


RESPONSE FACTOR REDUCED TO 0.0 


TRUNCATION FACTORS MODIFIED 
STRESS/STHA IN CONSTRAINT 
FREQUENCY CONSTRAINTS 


AS FOLLOWS 

0. 1 500 CUE400 
◦ • 1 20000E4 00 


UPDATED SCALING FAC T DP S 

0.225 OE 400 0.1 /1 4t 4 00 

0.9076E+00 0.1250E401 


0.2000L400 0.4720E400 0.5067E+00 0.2000E400 0.3350E401 

0.6376E400 


0.1800E40I 0*30006400 0*91346400 


UPDATED WEIGHT 
0 .4d4 2E *03 
0 .4122E 403 


COEFFICIENTS 
0.6456E40J 0.1076E403 

0.4/4bE40J 0.4283E403 


0. 63 46£ 4 04 

0.0 


0.81 75E 404 0.53786403 0*36026404 0*19356404 

0.4800E4Q5 


0*16X36403 0*23386403 


get 



RESPONSE FACTOR REDUCED TO 0.0 
TRUNCATION FACTORS MODIFIED AS FOLLOWS 

SJfLE-SS/ STRAIN CONSTRAINT 0.15QOOOEFOO . . .. - 

FREQUENCY CONSTRAINTS 0.120000E+00 

UPDATED SCALING FACTORS 

_ 0.225 QEtOO 0.171 4E + QQ Q.20QOE+QO 0.4720EFDO 0.5Q67ET00 O.2QOQET0O 0.3350EF01 Q* 1 EOOEFQl Q^3000EFQQ 0.9134EF0O 

0.9076E+00 0.1 2S0 E + 0 1 0.6376E«-00 

UPOATEO WEIGHT COEFFICIENTS 

0.48420*03 8.6456EFQ3 

0.4122E*03 0.4748ET03 


Q.1Q76ET03 0. 6346ETQ A 

0.42B3E+03 0.0 


0.81 75EF04 0 • 53 7EE* 03 

0.4800E+05 


0.3802EF04 O. 1 9 35 E* 04 0.1613EFQ3 



136 



i 


I 


STAGE NO. 6 APPROXIMATE PROBLEM GENERATOR 


CURRENT MEMBER SIZE 


MEMBER TYPE NUMBER 3 

O.IOOOE-OI O.IOOOE-OI 
0.4600E400 O.IOOOE-OI 

MEMB.ER ..TYPE .NUMBER 4. 

0.6823E+00 0.7736E400 

MEMBER TYPE NUMBER 6 

0.6823E+00 0.7736E400 


O.IOOOE-OI 
0 •5800E+00 

0.6715E+00 

0.6715E+00 


0. 1000E-01 
0.1 1 00E+00 

0.6715E+00 

0. 6715E+00 


0. 5600E+00 
0.1100E+00 

0. 570 OE F 00 

0.5700EF00 


0 . 46 0 OE 4 00 
0. 100CE-01 

0.570 CE4 00 

0 . 57 00E4 00 


0.4600E400 


O.IOOOE-OI 


0.5600E400 


0.4600E400 


CURRENT WEIGHT DATA 

ME MBER TY PE NUMBER 3 . WEIGHT * 0.205963E405 

MEMBER TYPE NUMBER 4 WEIGHT = 0.193988E404 


VARIABLE STRUCTURAL WEIGHT 0 .22536 IE + 05 

TOTAL STRUCTURAL WEIGHT - * 0.22 5 36 FeVoS 

NON-STRUCTURAL WEIGHThT 0.480000E405 


TOTAL WEIGHT 0.70S361E405 

CONVERGENCE CHECK STAGE NG.= 6 0.1077E-02 0.8Q10E-02 MUST BE LESS THAN 0.500000E-02 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0.2274 1 IE405 0.225604E405 0 .225361 E405 


ENTER OVERLAY ANALYS 


NEW AVAILABLE REAL ARRAY = 7473 


AVAILABLE INTEGER ARRAY - 2500 OVERLAY ANALYS REQUIREMENT a 772 


POSTURE TABLE 


RETAINEO 

TOTAL 

Type 

MEMBER 

NODE DIRECTION 

1 

! 

J! 
• | 
J 1 

l 

MODE 

CONSTRAINT values 

stress/strain 

CONSTRAINTS 

MOST CRITICAL 

= 0.742754E-01 




1 

54 

3 

1 1 


i 

3 

0. 2647 76 E 400 

2 

95 

J 

2 


2 

3 

0. I9 1498E + 00 

3 

105 

3 

4 


2 

3 

0.1591 I5E400 

4 

115 

3 

6 


2 

3 

0. 101879E+00 

5 

145 

3 

12 


2 

3 

0. 14 1490E+00 

6 

155 

J 

14 


2 

3 

0.742754E-01 

7 

165 

3 

16 


2 

3 

0.960189E-01 

FREQUENCY CONSTRAINTS 

MOST 

CRITICAL = - 

0.498237E— 02 




8 

1 73 





-1 

-0.498237E-02 


Mf|r>F RT/VNDS FOR THE F Cl I OWING 


4ei 


MOST CRI TICAL 


0.49U237E-02 


f-“tauENi_r cunoTha in : s 
H I 73 


-I -0.49B237E-02 


MODE STANDS FOR THE FOLLOWING 

NEGATIVE -LOWER BOUND POS I T I V E = UPP ER BOUND 

FOR STRESS CONSTRAINT, (CCDEFI) 

1 = VCN MISES EQUIVALENT STRESS 

2 * LCNGI T UDINAL STRAIN 

3 = transverse strain 

4 = S FEAR STRAIN 

5 = FIRST EQUATION OF STRESS INTERACTION 
6s SECOND EQUATION OF STRESS INTERACTION 
7 = T S A I — A ZZ I CRITERION 

FOR FREQUENCY CONSTRAINTS, ASSOCIATED MODE NUMBER 


AVAILABLE REAL ARRAY = 7500 OVERLAY ANALYS REQUIREMENT* 976 


END OVERLAY ANALYS 

_ SCALl.NG FACTOR Q.lOOSOlLtOl SCALED WEIGHT 0.226490EF05 

SIOE CONSTRAINTS 


RELATIVE MOVE LIMIT 0.8333EF00 




- VARIABLE. . 

LOWER 

ACTUAL 

UPPER 

VARIABLE 

LOWER 

ACTUAL 

UPPER . 



NUMBER 

BOUND 

SIZE 

BOUND 

NUMBER 

BOUND 

SIZE 

BOUNO 


JUS, 1 . ..3* 1 QQQE— Q l ..Q.IOOOE-Ol 0.20Q0E-01 DIS. 2 0 , 1 OOOE-Q1 _ 0. I OOOErQI 0.2000E-Q1 


01 5* 

3 

0. iOOOE-Ol 

0. IOOOE-Ol 

0.2000E-01 

DIS. 

4 

0 .4600E+00 

0 • 560 OE 400 

0.6800E400 

DIS. 

5 

0.3800E+00 

0 .460 0 1 *00 

0 .5600E+00 

DIS. 

6 

O.tOOOE-Ol 

0.1000E— 01 

0.2000E-01 

01 S. 

? 

0. 4 dO OE-*- 00 

0 .580QE+00 

0.7000E+00 

DIS. 

6 

0.9000E— 01 

0.11 OOE 400 

0. 1 400E+00 

JUS*,. 

__ 9 

.. 0. 1UQQE-Q1 

0. IOOOE-Ol 

Q « 2 QQ OE— 0 1 


10 

0.5686Et00 

0 .6B23EFQC 

O.BiBBEtQQ 


1 1 

0.6447E+00 

0.7736E+00 

0. 9283E+00 


12 

0.5596E+00 

0.67 1 5£ 40 0 

0.8038E400 


13 

0.4750E+00 

0.5700E+00 

0.6840E+00 







MOST VIOLATEO SIOE CONSTRAINT""- DESIGN VARIABLE " 9 CONSTRAINT VALUE 0.0 


ENTER OVERLAY PREDU1 


AVAJLLAHLE — RE At ASSAY = „ 75Q0 . OVERLAY PREDU 1 REQUIREMENT* 3872 

AVAILABLE INTEGER ARRAY* 2500 OVERLAY PREDU 1 REQUIREMENT* 746 


= =* = = OUALI OPTIMIZER ***** 


dual variables : 




0.0 0.0 

J .0 

0.0 

0.0 

o 

. 

o 

o 

. 

o 

0. 

291 18E40S 


PHASE 

C D A 

DUAL 3HJ. 

I GHT 

ZMCD 

NACT 

NPLAN 

ICON J NALFA 

NEWTON 

NO I S 

IDIS 

| 

1 

0.224ot*0b 

0 «22b OE Ob 

0. !06dE-0i 

1 

0 

0 0 

0 

2 

S 

1 

2 

0<*2246E*05 

0. 225 CE4- 05 

0.0 

I 

1 

0 0 

0 

2 

s 


DUAL VAR I ABLES : J t 
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D U A L 1 OPTIMIZER 


dual, variables : 

0.0 0.0 0.0 0.0 

PHASE COM DUAL OBJ. WEIGHT ZMCD 

1 1 0.2246E+05 0.22S0E+0S 0.106BE-01 

1 2 0.2246E+05 0.225CE+05 0.0 

. DUAL -VARI ABLES : ... 

0.0 0.0 0.0 0.0 

PHASE 00M DUAL OBJ. WEIGHT ZMOD 

2 l 0.224OE+05 0.2250E+05 0.0 

NUMBER OF RESTARTS 2 MAXPH = 

TOTAL NUMBER OF O.D.M. 3 MAXOOM = 

NUMdER OF NCN-ZERO DUAL VARIABLES 1 NTCE = 

NUMBER OF DISCONTINUITY PLANES 1 NLDV = 

DUAL OBJECTIVE FUNCTION 0.224646E+05 EPSPH = 

NORM OF (PROJECTED) GRADIENT 0.0 EPSODM = 

DUAL VARIABLES : 

- a.a o.o o.a o.o 

ANALYSIS OF DUAL SOLUTION 

Inumber of oiscontTnui ty pl anes i 

NUMBER OF POSSIBLE PRIMAL POINTS 2 

WEIGHT OF UPPER BOUND SOLUTION 0.224957E+05 FEASIBLE 

PRIMAL VARIABLES 

Q . 1 QQQQEtQ 1. . O.IOQCQE+Ol O.iaOOQE+ai O.IOOQOE+Ql 

O.IOOOOE+Ol 0 • l 0366E+01 0.1 0490E +0 1 0.99547E+00 

CONSTRAINTS 

- -Q - .2.69AXEJ-0-0 -0.19952E+00 -0.16876E+00 -0 . 10666E+00 

WEIGHT OF LOWER BOUND SOLUTION 0.222B06E+05 

PRIMAL VARIABLES. 

0.1 0 0 0 0 L +0 l O.IOOOOE+Ol O.IOOOOE+Ol 0 . X 00 OOE + 0 1 

0.1 OOOOE+Ol 0. 1 0366E+01 J.10490E+01 0.99547E+00 

CCNSTKAINTS 

-0.26770E+00 -0. 199206+00 -0.168226+00 -0.1036BE+00 

WEIGHT OF FINAL OESIGN 0.224957E+05 FEASIBLE 

NUMBER OF PRIMAL VARIABLES FROM LOWER BOUND SOLUTION 0 

PRIMAL VARIABLES 

0.1 OOOOE+Ol O.IOOOOE+Ol 0.1 OOOOE +0 1 O.IOOOOE+Ol 

0.1 000 OE + 01 0.10366E+01 0.10490E+01 0.9S547E+00 

CONSTRAINTS 

—0.2694 1 E + 0 0 -0.19952E+00 -O.X GB76E+00 - 0 • 1 OOOuE 1 00 


0.0 

NACT 

0.0 

nplan 

ICON J 

0. 0 
nalfa 

0. 

NEWTON 

291 18E+0S 
NO I S 

IDIi 

1 

0 

0 

0 

0 

2 

~5 

1 

1 

0 

0 

0 

2 

5 

0.0 

NACT 

0.0 

NPLAN 

ICONJ 

0.0 

nalfa 

0. 

NEWTON 

29830E+05 

NOIS 

IDIS 

1 

1 

0 

0 

0 

'2 

s’ 


1 0 
100 
6 
13 

0.1 OOOOOE— 04 

0*10 CO OOE — 03 

0.0 0.0 0 . 0 0 . 2983 QE +05 


0.XQ6E2E + 01 Q.1QQ00EFQ1 0. R27S 9E+00 

0.104 62 E + 0 1 


CO 

>£> 


-Q. 146E3E+0Q -0.42562E-01 -0.A8527E-QX -O.X0435E-Q2 


0.10435E+01 O.IOOOOE+Ol 0.82759E+00 0.818 16E+00 

O.X0462E+01 

— 0. 14372E+00 -0.40227E-01 -0.48694E-01 0.61680E-02 


0.10652E+01 0.1 OOOOE+Ol 0.62759E+00 0.81816E+00 

0 • X0462E +0 1 

-0.X4653E+00 - 0 . 42S62E-01 -Q.48527E-0I -0.X0435E-02 


END OVERLAY PREDLl 





0.4000E-01 _.Q.2400Et0l . .Q*9QQQE*.OQ_jO*10JOO£tA<L....a.«li79e«-.Ol 

0.I076E4C3 O.250OE4O4 0.96 77E+03 0.3376E+02 _Q»30I7E403 


f 
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t I V 


STAGE NO. 7 APPROXIMATE PROBLEM GEN ERA T OR 


CURRENT MEMBER SIZE 


MEMBER TYPE NUMBER 3 

O.IOOOE-Ol O.IOOOE-Ol O.IOOOE-Ol O.IOOOE-Ol 0.5600E+00 0.490CEWC0 0.4900E400 O.IOOOE-Ol 0.5600E4-00 0.4900E+00 

0.490 OE +0 0 O.IOOOE-Ol 0.4800E+00 G.9000E-01 0.9000E-01 0.100CE-01 

MEMBER TYPE NUMBER 4 

0.7C73E+00 0 . ti 1 1 5E+ 00 0.6b8bE+00 0.6685EH00 0.5963E+00 0.5963EW-00 

MEMBER T yPE NUMBER 6 

0.707 3E +00 0.81 1 5E TO 0 0.6685E+00 0.6665E+00 0.5963E+00 0.5963EFQ0 


CURRENT HEIGHT DATA 


MEMBER TYPE NUMBER 3 WEIGhT = 0.204890E+05 

MEMBER TYPE NUMBER 4 HEIGHT = 0.200677E*-04 


VARIABLE STRUCTURAL WEIGHT 0.224957E+05 
FlJtEJ STRUCTURAL WEIGHT 0.0 

TOTAL STRUCTURAL WEIGHT 0.224957EW05 

NON- STRUCTURAL WEIGHThT 0.480000EF05 


... ... . TOTAL WEIGHT 0.704957E+Q5 

CONVERGENCE CHECK STAGE N0.= 7 0.1796E-02 0.I077E-02 MUST BE LESS THAN 0.500000E-02 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARE 0 . 2256 04E «■ 05 0 .2 2536 1 EF05 0.224957EF05 


ENTER OVERLAY ANALYS 


NOCAL DISPLACEMENTS 

NODE X Y i NODE 

LOAD CCNCIT'N l 

1 0.0 0.0 0.0 2 0.0 0.0 0.0 

3 0.0 0.0 0.0 4 -0.1431 OEWOO 0.42576E-01 0.38610E401 

-> -0.40O6SE + 00 0.1 9068E + 00 0. 15345E + 02 6 - 0. 76908EF0 0 0.52852EF00 0.68358E+02 

L C AD CCNC I T I CN 2 

1 0.0 0.0 0.0 

J J.U 0.0 0.0 


2 0.0 0.0 0.0 
4 0. 1431 QE+OO -0.42576E-01 - 0 . 386 1 OE ♦ 0 1 
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NOCAL DISPLACEMENTS 


NODE 

X 

y 

i 

NODE 

X 

Y 

z 

L C AO CCNC I T I :.N 1 







1 

0.0 

0.0 

0. 0 

2 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

4 

-0. 1431 OE+OO 

0.42576E-01 

0.3B610E401 

‘i 

-0.40U65E+00 

0 .1 9066E+00 

0. 1534SE+02 

6 

-0.76908E+00 

0.52852E+00 

0.68358E+02 

L C AD CCNDI T 

ICN 2 







1 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

J 

J.O 

0.0 

0. 0 

4 

0.1431 OE*00 

-0.42576E-01 

-0.38610E+01 

5 

j.4oao‘SE*-oo 

-0.1 9JoaE»oo 

-0. 1534SE*02 

6 

0. 76908E*00 

-0.52B52E*00 

— 0.68358 E 4 02 

H 

II 

II 

II 

II 

II 

It 

II 

II 

II 


=;rss=33=ss3=s 

:5=SSS333SS233:333 

II 

II 

II 

II 

II 

II 

u 

N 

II 

II 

II 

w 

II 

S3 X33 1SIISSSI3 

sasxxiissisasa siaiui sa n a ■■ 





NOCAL 

DISPLACEMENTS 




N00E 

X 

Y 

2 

NODE 

• X 

Y 

z 

LOAD CCNDl T I3N 

i 







1 J.O 


0.0 

0. 0 

2 

0.0 

0.0 

0.0 

3 0.0 


0.0 

0.0 

4 

0.29981E-01 

-0.21788E400 

—0 . 55807 E— 02 

5 - J.97t>8Jd-0i 

0 .1 2948E+00 

-0. 34034E-02 

6 

— 0. 9Q203E— 0 1 

0.2621 IE— 0 1 

— 0.2336 OE— 02 

LC AO CCND ITICN 

2 







l 0. 0 


0.0 

0. 0 

2 

0.0 

0.0 

0.0 

1 

l 




■ 


n rc6rt>»e_A»» 
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NODAL DISPLACEMENTS 


NODE X Y t NODE X Y Z 


LOAD C CND I T I ON 1 


1 

3 

5 


0.0 

0.0 

— 0.975 83 E— 01 


0.0 

0.0 

0.129486*00 


0.0 

0.0 

-0.34034E-02 


2 

4 

6 


0.0 

0.2 9 98 IE— 01 
— 0.98203E— 01 


0.0 

—0.21 788E400 
0.2621 IE-01 


0.0 

— 0.SS807E— 02 
— 0. 23360E— 02 


LOAD CONDITION 2 


1 

J 

5 




0. 0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

O.Q 

0 • 0 

0. 0 

4 

0.2998 1 E— 0 1 

-0.21788E400 

— 0.55807E— 02 

0 .97583E-01 

0.1 2S46E+00 

— 0.34 034 E— 02 

6 

— 0.98203E— 0 1 

0.2621 IE— 0 1 

-0. 23360 E-02 


rSS3SSSSSSSX53SSSa3<SZ23r33SS3S>X83SIISHS33aS3a33Z8Sa 333 3 33 331333333 3333333 33 ..1X3.3. .1333.1.333331* 


NEW AVAILABLE REAL ARRAY = 74 73 


EIGEN VECTORS SCALEO BY MAX. COMPONENTS 

VECTOR NO. = 1 FREQUENCY = 0.1005026*01 C/S 

^a.l381£r^Q2„ 0*581 9E-Q3 Q. 41 756-01 -0.5727E-02 O.2836E-02 0.20856*00 -0.1181E-01 0.80356-02 0. 10006*01 

VECTOR N0.= 2 FR EQUENCY = 0.395797E+01 C/S 

— 0 . 1471 E— 01 0.2455E-02 0.4279E*00 -0.1342E-01 0.8765E-02 0.10006*01 0.4S40E-01 0.1912E-02 -0.69796400 
VECTOR NO . — 3 FREQUENCY = 0.543532E*01 C/S 

35106- Q l -0.2519E-Q3 0.1 QQOEtQl Q.5961E-C2 -0.34406-01 -O.S3O7E*O0 -Q.1697E-Q1 -0.6252E-01 0.3 832E-01 


***************4**********4************* EIGEN ANALYSIS CONVERGED IN 2 ITERATIONS 4*4**4*4* *••******•*•***•♦»•*•••******•• 


EIGEN VALUES 

0.3988E*02 0.6185E*03 0.1166E+04 


EIGEN VECTORS SCALEO BY U MU 

VECTOR NO. l FREQUENCY= 0.1005026*01 C/S 


—0. 24 1 8E— 03 

0.8897E— 04 0.6384E 

S— COMBINED 

-02 -0.8756E-03 
SX 

0. 4336 E- 03 
SY 

0.3188E-01 

-0.1805E-02 0 

.12296-02 

}.1S29E*00 

MTYP 

M 

LC 

SXY 

SX- THERM 

SY-THERM 

SXV-THERN 

3 

l. _. 



— 0.64 0 IE— 03 

0.8129E— 07 

0. 1 183 Er 03 




3 

1 

1 


— 0 • 606 8 E— 03 

0.3200E-02 

0.6624E— 03 

0. 333 OE— 04 

0.3200E— 02 

0.5441 E— 03 

3 

2 

1 


-0 . 37 94 E— 03 

-0.2606E-03 

— 0. 640 1 E— 03 




3 

2 

1 


-0.6556E-03 

0.3249E— 02 

-0. 56476—03 

— 0.2762E— 03 

0.3S10E— 02 

0 • 7 538 E— 04 

3 

3_ 

1. 


-0 .2607E— 03 

— 0.3793E— 03 

Qa 6401 E- 03 




3 

3 

1 


0.714 7E— 05 

0. 25866-02 

0 . 5649E— 03 

0.2679E-03 

0.2966E— 02 

-0. 75216-04 

3 

4 



-0.6464E-11 

— 0. 6400E— 03 

-0. 1 191 E— 03 




3 

4 

1 


— 0.4200E— 04 

0.2636E-02 

— 0.6631 E— 03 

— 0. 4200E— 04 

0.3276E— 02 

-0.54406-03 

3 

5 



-0. 6401E-03 

— 0.3495E— 03 

-0.1377E-02 




3 

5 



— 0 . 6 06 8 E— 03 

0.2127E-02 

-0. 1 098E-02 

0 • 3330E— 04 

0.24766-02 

0.27876—03 

3 

6 

1 


0. 1935E-03 

-0.1 183E-02 

-0.291 1 E— 03 




3 

6 

1 


— 0.311 9E— 03 

0,1 8326—02 

0. 5 080 E— 03 

— O.5054E— 03 

0.30156-02 

0.79916-03 

3 

7. 

1 


-0. 1183E-02 

0.1937E-03 

0. 2900E— 03 




3 

7 



-0.1410E-02 

0.2930E— 02 

—0. 5 08 9 E— 03 

-0.2267E-03 

0.2736E-02 

— 0.7989E— 03 

3 

<3 

1 


-0.34V6E-03 

— 0. 640 IE — 03 

0. 13776-02 




3 

a 

1 


-0 . 11 1 5E-02 

0.2635E-02 

0. 1 098E-02 

-0.7657E-03 

0.32756-02 

— 0.2788E— 03 

J 

) 

I 


-0. 1055E-02 

0.32166-07 

-0. 1 451 6-02 




3 

9 

1 


-0. 1 34 IE— 02 

0.3200E-02 

-0. 1 774E-02 

-0.2859E-03 

0.3200E-02 

-0 .32346—03 

3 

l 0 

1 


0. I982E-03 

-0. 1253E-02 

-0.1055E-02 




3 

10 



0. 196 OE— 03 

0.1663E-02 

-0. 1 298 E— 02 

-0 . 2245E-05 

0.2916E— 02 

— 0.2438E— 03 

3 

1 1 

1 


— 0 . 12S3E— 02 

0.1979E-03 

0. 1 0566— 02 




3 

1 1 



-0. 1579E-02 

o*j43ae-o2 

0. 1300E-02 

-0.3256E- 03 

0.32406—02 

0.24416-03 

3 

1 2 

I 


—0 . 204 IE— 08 

-0.1055E-02 

0. 145 IE- 02 




3 

1 2 

1 


—0 .42006—04 

0. 19016-02 

0. 1 775E-02 

-0.42006-04 

0.29566-02 

0.3234E— 03 

3 

1 3 

1 


-0 .1280E-02 

— 0.3 496 E- 03 

-0. 169JE-02 




3 

1 3 

1 


-0. 1 J23E-02 

0.2 1276—02 

-0. 1 70 1 E— 02 

— 0.4356E— 04 

0.24766-02 

—0 • 7605E— 05 

3 

1 4 



0.3203E-04 

—0.1 66 1 E— 02 

— 0. 9202E— 03 




J 

1 4 

1 


— 0 • 36QQE— 0 J 

0. 1 1 726-02 

— 0.2 08 1 E— 03 

-0.40086-03 

0.28346-02 

0. 722 IE— 03 

3 

1 ii 

1 


-0. 166 IE— 02 

0. 3203E — 04 

0.93026-03 




3 

1 L S 

1 


— 0.2070E— 02 

0.2873E-02 

0.208 IE- 03 

—0. 4084E— 03 

0 .284 IE— 02 

— 0 . 7221 E— 03 

3 

1 6 

t 


—0. 349 6 E— 03 

-0.1280E-02 

0. 1 693 E— 02 




3 

I 6 

I 


-0. 1 1 1 5E-02 

0. 19196-02 

0. 1 70 1 E— 02 

-0.76576-03 

0.3198E-02 

0.7596E— 05 

4 - 

1 

1 

0. 3238926+05 

-0.50606*04 

0*0 

0. 18476*05 




6 

I 

1 

0 • 3 3 OO n * 


“ 






143 


J 

8 

1 

J 

9 

1 

3 

9 

l 

3 

1 0 

1 

3 

1 0 

1 

3 

1 L 

1 

3 

1 t 

1 

3 

1 2 

I 

3 

l 2 

1 

... 

-La 


3 

13 

1 

3 

1 4 

1 

3 

I 4 

l 

3 

15 


3 

l 5 

1 

3 

16 

1 

3 

1 6 

1 

A 

l 


6 

1 

1 

4 

2 

l 

6 

2 

l 

. 4 

a 

1 

6 

3 

1 

4 

4 

1 

6 

4 

1 

A 

s 

a 

6 

5 

1 

4 

6 

1 

6 

6 

1 

■j 

X 

_2_ 

3 

1 

2 

3 

2 

2 

3 

2 

2 

_3_ 

3 

P 

3 

3 

2 

3 

4 

2 

3 

4 

2 

. _3 

5 

2 

3 

5 

2 

3 

6 

2 

3 

6 

2 

a 

. 7 

2 

3 

7 

2 

3 

8 

2 

3 

a 

2 

3 

Q 

2 

3 

9 

2 

3 

10 

2 

3 

10 

2 

3 

1 1 

g 

3 

1 1 

2 

3 

1 2 

2 

3 

12 

2 

3 

1 3 _ 


3 

13 

2 

3 

14 

2 

3 

14 

2 

-3 

15 

2 

3 

1 5 

2 

3 

16 

2 

J 

16 

2 

_4._ 

L 

.2 

6 

1 

2 

4 

2 

2 

6 

2 

2 

4 

3 

2 

6 

3 

2 

4 

4 

2 

6 

4 

2 

4 

5 

2 

6 

5 

2 

4 

6 

2 

b 

6 

2 



-0 • 1 1 1 bL-02 

0.26j3t-u2 

0 * 1 oyat - 02 

-0.7657E-0J 


-0. 1055E-02 

0. J216E-07 

— 0.1 45 1 6— 02 



-0. 134 IE— 02 

0. 3200E-02 

-0. 1 7746-02 

-0.2859E-03 


0. 1982E-03 

-0. i 25 BE— 02 

-0. 10556-02 



0. 19606-03 

0.1 663E-02 

-0. 12986-02 

-0.2245E-05 


-0 • 1253E-02 

0.1979E-C3 

0. 10566- 02 


-0.15796-02 

0 .34 3 BE- 02 

0.13006-02 

-0.3256E-03 


-0. 20416-08 

-0.1 05SE-02 

0. 14516-02 



- 0 • 4200E— 04 

0. 1901E-02 

0. 17756-02 

-0.4200E-04 


-0 .1280E-02 

-0.3496E-03 

-0. 1C93E-02 


-0. 1323E-02 

0.2 127E-02 

-0. 1 7016-02 

-0.4356E-04 


0.3203E-04 

— 0.1661 E— 02 

-0.93026-03 



— 0 • 368 8 E— 03 

0.11 72E-02 

-0.20816-03 

-0.4008E-03 


-0. 166 IE— 02 

0 . 3203E— 04 

0.93026-03 



— 0 • 2070E— 02 

0.2 8 73 E— 02 

0.20816-03 

-0.4084E-03 


—0. 34 9 6 E— 03 

—0 • l 28 OE — 02 

0. 16936-02 


0^3238926 + 05 

-0. 11 15E-02 

0. 1919E-02 

0.17016-02 

-0.7657E-03 

— 0.5060E+04 

0.0 

0.18476+05 

0.3300656+05 

-0.93556+04 

-0.5367E*04 

0. 18476+05 

“0 . 42946*04 

0.3980126+05 

-0.58676+04 

0.0 

0.22736+05 

0. 402 6 73E* 05 

0.73286+04 

0.9274E*04 

0.22736*05 

0. 1320E*0S 

0.1676836*05 

-0. 16756+05 

0. 0 

-0.39016* 03 

0.1382196+05 

-0.976 16+04 

0.6034E404 

-0.39016+03 

0.6993E*04 

0. 4325476+05 
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